This white paper submitted for the 2020 Canadian Long-Range Planning process (LRP2020) presents the prospects for Canada and the Square Kilometre Array (SKA) from 2020-2030, focussing on the first phase of the project (SKA1) scheduled to begin construction early in the next decade. SKA1 will make transformational advances in our understanding of the Universe across a wide range of fields, and Canadians are poised to play leadership roles in several. Canadian key SKA technologies will ensure a good return on capital investment in addition to strong scientific returns, positioning Canadian astronomy for future opportunities well beyond 2030. We therefore advocate for Canada's continued scientific and technological engagement in the SKA from 2020-2030 through participation in the construction and operations phases of SKA1.

We contrast predictions for the high-redshift galaxy population and reionization history between cold dark matter (CDM) and an alternative self-interacting dark matter model based on the recently developed ETHOS framework that alleviates the small-scale CDM challenges within the Local Group. We perform the highest resolution hydrodynamical cosmological simulations (a 36~Mpc$^3$ volume with gas cell mass of $\sim10^5\mathrm{M}_{\odot}$ and minimum gas softening of $\sim180$~pc) within ETHOS to date -- plus a CDM counterpart -- to quantify the abundance of galaxies at high redshift and their impact on reionization. We find that ETHOS predicts galaxies with higher ultraviolet (UV) luminosities than their CDM counterparts and a faster build-up of the faint end of the UV luminosity function. These effects, however, make the optical depth to reionization less sensitive to the power spectrum cut-off: the ETHOS model differs from the CDM $\tau$ value by only 10 per cent and is consistent with Planck limits if the effective escape fraction of UV photons is 0.1-0.5. We conclude that current observations of high-redshift luminosity functions cannot differentiate between ETHOS and CDM models, but deep JWST surveys of strongly-lensed, inherently faint galaxies have the potential to test non-CDM models that offer attractive solutions to CDM's Local Group problems.

The Canadian Hydrogen Intensity Mapping Experiment (CHIME) is a novel transit radio telescope operating across the 400-800-MHz band. CHIME is comprised of four 20-m x 100-m semi-cylindrical paraboloid reflectors, each of which has 256 dual-polarization feeds suspended along its axis, giving it a >200 square degree field-of-view. This, combined with wide bandwidth, high sensitivity, and a powerful correlator makes CHIME an excellent instrument for the detection of Fast Radio Bursts (FRBs). The CHIME Fast Radio Burst Project (CHIME/FRB) will search beam-formed, high time-and frequency-resolution data in real time for FRBs in the CHIME field-of-view. Here we describe the CHIME/FRB backend, including the real-time FRB search and detection software pipeline as well as the planned offline analyses. We estimate a CHIME/FRB detection rate of 2-42 FRBs/sky/day normalizing to the rate estimated at 1.4-GHz by Vander Wiel et al. (2016). Likely science outcomes of CHIME/FRB are also discussed. CHIME/FRB is currently operational in a commissioning phase, with science operations expected to commence in the latter half of 2018.

We formulate an effective theory of structure formation (ETHOS) that enables cosmological structure formation to be computed in almost any microphysical model of dark matter physics. This framework maps the detailed microphysical theories of particle dark matter interactions into the physical effective parameters that shape the linear matter power spectrum and the self-interaction transfer cross section of non-relativistic dark matter. These are the input to structure formation simulations, which follow the evolution of the cosmological and galactic dark matter distributions. Models with similar effective parameters in ETHOS but with different dark particle physics would nevertheless result in similar dark matter distributions. We present a general method to map an ultraviolet complete or effective field theory of low energy dark matter physics into parameters that affect the linear matter power spectrum and carry out this mapping for several representative particle models. We further propose a simple but useful choice for characterizing the dark matter self-interaction transfer cross section that parametrizes self-scattering in structure formation simulations. Taken together, these effective parameters in ETHOS allow the classification of dark matter theories according to their structure formation properties rather than their intrinsic particle properties, paving the way for future simulations to span the space of viable dark matter physics relevant for structure formation.

Our understanding of the Universe is known to be incomplete and new gauge forces beyond those of the Standard Model might be crucial to describing its observed properties. A minimal and well-motivated possibility is a pure Yang-Mills non-Abelian dark gauge force with no direct connection to the Standard Model. We determine here the relic abundances of the glueball bound states that arise in such theories and investigate their cosmological effects. Glueballs are first formed in a confining phase transition, and their relic densities are set by a network of annihilation and transfer reactions. The lightest glueball has no lighter states to annihilate into, and its yield is set mainly by 3 to 2 number-changing processes which persistently release energy into the glueball gas during freeze-out. The abundances of the heavier glueballs are dominated by 2 to 2 transfer reactions, and tend to be much smaller than the lightest state. We also investigate potential connectors between the dark force and the Standard Model that allow some or all of the dark glueballs to decay. If the connection is weak, the lightest glueball can be very long-lived or stable and is a viable dark matter candidate. For stronger connections, the lightest glueball will decay quickly but other heavier glueball states can remain stable and contribute to the dark matter density.

We present results from a new incoherent-beam Fast Radio Burst (FRB) search on the Canadian Hydrogen Intensity Mapping Experiment (CHIME) Pathfinder. Its large instantaneous field of view (FoV) and relative thermal insensitivity allow us to probe the ultra-bright tail of the FRB distribution, and to test a recent claim that this distribution's slope, $\alpha\equiv-\frac{\partial \log N}{\partial \log S}$, is quite small. A 256-input incoherent beamformer was deployed on the CHIME Pathfinder for this purpose. If the FRB distribution were described by a single power-law with $\alpha=0.7$, we would expect an FRB detection every few days, making this the fastest survey on sky at present. We collected 1268 hours of data, amounting to one of the largest exposures of any FRB survey, with over 2.4\,$\times$\,10$^5$\,deg$^2$\,hrs. Having seen no bursts, we have constrained the rate of extremely bright events to $<\!13$\,sky$^{-1}$\,day$^{-1}$ above $\sim$\,220$\sqrt{(\tau/\rm ms)}$ Jy\,ms for $\tau$ between 1.3 and 100\,ms, at 400--800\,MHz. The non-detection also allows us to rule out $\alpha\lesssim0.9$ with 95$\%$ confidence, after marginalizing over uncertainties in the GBT rate at 700--900\,MHz, though we show that for a cosmological population and a large dynamic range in flux density, $\alpha$ is brightness-dependent. Since FRBs now extend to large enough distances that non-Euclidean effects are significant, there is still expected to be a dearth of faint events and relative excess of bright events. Nevertheless we have constrained the allowed number of ultra-intense FRBs. While this does not have significant implications for deeper, large-FoV surveys like full CHIME and APERTIF, it does have important consequences for other wide-field, small dish experiments.

Cosmological observations are becoming increasingly sensitive to the effects of light particles in the form of dark radiation (DR) at the time of recombination. The conventional observable of effective neutrino number, $N_{\rm eff}$, is insufficient for probing generic, interacting models of DR. In this work, we perform likelihood analyses which allow both free-streaming effective neutrinos (parametrized by $N_{\rm eff}$) and interacting effective neutrinos (parametrized by $N_{\rm fld}$). We motivate an alternative parametrization of DR in terms of $N_{\rm tot}$ (total effective number of neutrinos) and $f_{\rm fs}$ (the fraction of effective neutrinos which are free-streaming), which is less degenerate than using $N_{\rm eff}$ and $N_{\rm fld}$. Using the Planck 2015 likelihoods in conjunction with measurements of baryon acoustic oscillations (BAO), we find constraints on the total amount of beyond the Standard Model effective neutrinos (both free-streaming and interacting) of $\Delta N_{\rm tot} < 0.39$ at 2$\sigma$. In addition, we consider the possibility that this scenario alleviates the tensions between early-time and late-time cosmological observations, in particular the measurements of $\sigma_8$ (the amplitude of matter power fluctuations at 8$h^{-1}$ Mpc), finding a mild preference for interactions among light species. We further forecast the sensitivities of a variety of future experiments, including Advanced ACTPol (a representative CMB Stage-III experiment), CMB Stage-IV, and the Euclid satellite. This study is relevant for probing non-standard neutrino physics as well as a wide variety of new particle physics models beyond the Standard Model that involve dark radiation.

The recent detection of gravitational waves indicates that stellar-mass black hole binaries are likely to be a key population of sources for forthcoming observations. With future upgrades, ground-based detectors could detect merging black hole binaries out to cosmological distances. Gravitational wave bursts from high redshifts ($z \gtrsim 1$) can be strongly magnified by gravitational lensing due to intervening galaxies along the line of sight. In the absence of electromagnetic counterparts, the mergers' intrinsic mass scale and redshift are degenerate with the unknown magnification factor $\mu$. Hence, strongly magnified low-mass mergers from high redshifts appear as higher-mass mergers from lower redshifts. We assess the impact of this degeneracy on the mass-redshift distribution of observable events for generic models of binary black hole formation from normal stellar evolution, Pop III star remnants, or a primordial black hole population. We find that strong magnification ($\mu \gtrsim 3$) generally creates a heavy tail of apparently massive mergers in the event distribution from a given detector. For LIGO and its future upgrades, this tail may dominate the population of intrinsically massive, but unlensed mergers in binary black hole formation models involving normal stellar evolution or primordial black holes. Modeling the statistics of lensing magnification can help account for this magnification bias when testing astrophysical scenarios of black hole binary formation and evolution.

In the standard structure formation scenario based on the cold dark matter paradigm, galactic halos are predicted to contain a large population of dark matter subhalos. While the most massive members of the subhalo population can appear as luminous satellites and be detected in optical surveys, establishing the existence of the low mass and mostly dark subhalos has proven to be a daunting task. Galaxy-scale strong gravitational lenses have been successfully used to study mass substructures lying close to lensed images of bright background sources. However, in typical galaxy-scale lenses, the strong lensing region only covers a small projected area of the lens's dark matter halo, implying that the vast majority of subhalos cannot be directly detected in lensing observations. In this paper, we point out that this large population of dark satellites can collectively affect gravitational lensing observables, hence possibly allowing their statistical detection. Focusing on the region of the galactic halo outside the strong lensing area, we compute from first principles the statistical properties of perturbations to the gravitational time delay and position of lensed images in the presence of a mass substructure population. We find that in the standard cosmological scenario, the statistics of these lensing observables are well approximated by Gaussian distributions. The formalism developed as part of this calculation is very general and can be applied to any halo geometry and choice of subhalo mass function. Our results significantly reduce the computational cost of including a large substructure population in lens models and enable the use of Bayesian inference techniques to detect and characterize the distributed satellite population of distant lens galaxies.

The Canadian Hydrogen Intensity Mapping Experiment (CHIME) Pathfinder radio telescope is currently surveying the northern hemisphere between 400 and 800 MHz. By mapping the large scale structure of neutral hydrogen through its redshifted 21 cm line emission between $z \sim 0.8-2.5$ CHIME will contribute to our understanding of Dark Energy. Bright astrophysical foregrounds must be separated from the neutral hydrogen signal, a task which requires precise characterization of the polarized telescope beams. Using the DRAO John A. Galt 26 m telescope, we have developed a holography instrument and technique for mapping the CHIME Pathfinder beams. We report the status of the instrument and initial results of this effort.

We present the first simulations within an effective theory of structure formation (ETHOS), which includes the effect of interactions between dark matter and dark radiation on the linear initial power spectrum and dark matter self-interactions during non-linear structure formation. We simulate a Milky Way-like halo in four different dark matter models and the cold dark matter case. Our highest resolution simulation has a particle mass of $2.8\times 10^4\,{\rm M}_\odot$ and a softening length of $72.4\,{\rm pc}$. We demonstrate that all alternative models have only a negligible impact on large scale structure formation. On galactic scales, however, the models significantly affect the structure and abundance of subhaloes due to the combined effects of small scale primordial damping in the power spectrum and late time self-interactions. We derive an analytic mapping from the primordial damping scale in the power spectrum to the cutoff scale in the halo mass function and the kinetic decoupling temperature. We demonstrate that certain models within this extended effective framework that can alleviate the too-big-to-fail and missing satellite problems simultaneously, and possibly the core-cusp problem. The primordial power spectrum cutoff of our models naturally creates a diversity in the circular velocity profiles, which is larger than that found for cold dark matter simulations. We show that the parameter space of models can be constrained by contrasting model predictions to astrophysical observations. For example, some models may be challenged by the missing satellite problem if baryonic processes were to be included and even over-solve the too-big-to-fail problem; thus ruling them out.

In a popular class of models, dark matter comprises an asymmetric population of composite particles with short range interactions arising from a confined nonabelian gauge group. We show that coupling this sector to a well-motivated light mediator particle yields efficient darkleosynthesis, a dark-sector version of big-bang nucleosynthesis (BBN), in generic regions of parameter space. Dark matter self-interaction bounds typically require the confinement scale to be above \Lambda_{QCD}, which generically yields large (>>MeV/dark-nucleon) binding energies. These bounds further suggest the mediator is relatively weakly coupled, so repulsive forces between dark-sector nuclei are much weaker than coulomb repulsion between standard-model nuclei, which results in an exponential barrier-tunneling enhancement over standard BBN. Thus, dark nuclei are easier to make and harder to break than visible species with comparable mass numbers. This process can efficiently yield a dominant population of states with masses significantly greater than the confinement scale and, in contrast to dark matter that is a fundamental particle, may allow the dominant form of dark matter to have high spin > 3/2.

We propose that "standard pings", brief broadband radio impulses, can be used to study the three-dimensional clustering of matter in the Universe even in the absence of redshift information. The dispersion of radio waves as they travel through the intervening plasma can, like redshift, be used as a cosmological distance measure. Because of inhomogeneities in the electron density along the line of sight, dispersion is an imperfect proxy for radial distance and we show that this leads to calculable dispersion-space distortions in the apparent clustering of sources. Fast radio bursts (FRBs) are a new class of radio transients that are the prototypical standard ping and, due to their high observed dispersion, have been interpreted as originating at cosmological distances. The rate of fast radio bursts has been estimated to be several thousand over the whole sky per day and, if cosmological, the sources of these events should trace the large-scale structure of the Universe. We calculate the dispersion-space power spectra for a simple model where electrons and FRBs are biased tracers of the large-scale structure of the Universe and we show that the clustering signal could be measured using as few as 10 000 events. Such a survey is in line with what may be achieved with upcoming wide-field radio telescopes.

We present the design and implementation of a custom GPU-based compute cluster that provides the correlation X-engine of the CHIME Pathfinder radio telescope. It is among the largest such systems in operation, correlating 32,896 baselines (256 inputs) over 400MHz of radio bandwidth. Making heavy use of consumer-grade parts and a custom software stack, the system was developed at a small fraction of the cost of comparable installations. Unlike existing GPU backends, this system is built around OpenCL kernels running on consumer-level AMD GPUs, taking advantage of low-cost hardware and leveraging packed integer operations to double algorithmic efficiency. The system achieves the required 105TOPS in a 10kW power envelope, making it among the most power-efficient X-engines in use today.

In the standard model neutrinos are assumed to have streamed across the Universe since they last scattered at the weak decoupling epoch when the temperature of the standard-model plasma was ~MeV. The shear stress of free-streaming neutrinos imprints itself gravitationally on the Cosmic Microwave Background (CMB) and makes the CMB a sensitive probe of neutrino scattering. Yet, the presence of nonstandard physics in the neutrino sector may alter this standard chronology and delay neutrino free-streaming until a much later epoch. We use observations of the CMB to constrain the strength of neutrino self-interactions G_eff and put limits on new physics in the neutrino sector from the early Universe. Recent measurements of the CMB at large multipoles made by the Planck satellite and high-l experiments are critical for probing this physics. Within the context of conventional LambdaCDM parameters cosmological data are compatible with G_eff < 1/(56 MeV)^2 and neutrino free-streaming might be delayed until their temperature has cooled to as low as ~25 eV. Intriguingly, we also find an alternative cosmology compatible with cosmological data in which neutrinos scatter off each other until z~10^4 with a preferred interaction strength in a narrow region around $G_{\rm eff} \simeq 1/({\rm 10 \, MeV})^{2} \simeq 8.6\times10^8 G_{\rm F}$, where $G_{\rm F}$ is the Fermi constant. This distinct self-interacting neutrino cosmology is characterized by somewhat lower values of both the scalar spectral index and the amplitude of primordial fluctuations. While we phrase our discussion here in terms of a specific scenario in which a late onset of neutrino free-streaming could occur, our constraints on the neutrino visibility function are very general.

During and after recombination, in addition to Thomson scattering with free electrons, photons also coupled to neutral hydrogen and helium atoms through Rayleigh scattering. This coupling influences both CMB anisotropies and the distribution of matter in the Universe. The frequency-dependence of the Rayleigh cross section breaks the thermal nature of CMB temperature and polarization anisotropies and effectively doubles the number of variables needed to describe CMB intensity and polarization statistics, while the additional atomic coupling changes the matter distribution and the lensing of the CMB. We introduce a new method to capture the effects of Rayleigh scattering on cosmological power spectra. Rayleigh scattering modifies CMB temperature and polarization anisotropies at the $\sim\!1 \%$ level at $353 {\rm GHz}$ (scaling $\propto \nu^4$), and modifies matter correlations by as much as $\sim\!0.3\%$. We show the Rayleigh signal, especially the cross-spectra between the thermal (Rayleigh) E-polarization and Rayleigh (thermal) intensity signal, may be detectable with future CMB missions even in the presence of foregrounds, and how this new information might help to better constrain the cosmological parameters.

A sufficiently rigid relativistic elastic solid can be stable for negative pressure values and thus is capable of driving a stage of accelerated expansion. If a relativistic elastic solid drove an inflationary stage in the early Universe, quantum mechanically excited perturbations would arise in the medium. We quantize the linear scalar and tensor perturbations and investigate the observational consequences of having such an inflationary period. We find that slowly varying sound speeds of the perturbations and a slowing varying equation of state of the solid can produce a slightly red-tilted scalar power spectrum that agrees with current observational data. Even in the absence of nonadiabatic pressures, perturbations evolve on superhorizon scales, due to the shear stresses within the solid. As such, the spectra of perturbations are in general sensitive to the details of the end of inflation and we characterize this dependence. Interestingly, we uncover here accelerating solutions for elastic solids with (1 + P/\rho) significantly greater than 0 that nevertheless have nearly scale-invariant scalar and tensor spectra. Beyond theoretical interest, this may allow for the possibility of viable inflationary phenomenology relatively far from the de Sitter regime.

The Canadian Hydrogen Intensity Mapping Experiment (CHIME) is a transit interferometer currently being built at the Dominion Radio Astrophysical Observatory (DRAO) in Penticton, BC, Canada. We will use CHIME to map neutral hydrogen in the frequency range 400 -- 800\,MHz over half of the sky, producing a measurement of baryon acoustic oscillations (BAO) at redshifts between 0.8 -- 2.5 to probe dark energy. We have deployed a pathfinder version of CHIME that will yield constraints on the BAO power spectrum and provide a test-bed for our calibration scheme. I will discuss the CHIME calibration requirements and describe instrumentation we are developing to meet these requirements.

A pathfinder version of CHIME (the Canadian Hydrogen Intensity Mapping Experiment) is currently being commissioned at the Dominion Radio Astrophysical Observatory (DRAO) in Penticton, BC. The instrument is a hybrid cylindrical interferometer designed to measure the large scale neutral hydrogen power spectrum across the redshift range 0.8 to 2.5. The power spectrum will be used to measure the baryon acoustic oscillation (BAO) scale across this poorly probed redshift range where dark energy becomes a significant contributor to the evolution of the Universe. The instrument revives the cylinder design in radio astronomy with a wide field survey as a primary goal. Modern low-noise amplifiers and digital processing remove the necessity for the analog beamforming that characterized previous designs. The Pathfinder consists of two cylinders 37\,m long by 20\,m wide oriented north-south for a total collecting area of 1,500 square meters. The cylinders are stationary with no moving parts, and form a transit instrument with an instantaneous field of view of $\sim$100\,degrees by 1-2\,degrees. Each CHIME Pathfinder cylinder has a feedline with 64 dual polarization feeds placed every $\sim$30\,cm which Nyquist sample the north-south sky over much of the frequency band. The signals from each dual-polarization feed are independently amplified, filtered to 400-800\,MHz, and directly sampled at 800\,MSps using 8 bits. The correlator is an FX design, where the Fourier transform channelization is performed in FPGAs, which are interfaced to a set of GPUs that compute the correlation matrix. The CHIME Pathfinder is a 1/10th scale prototype version of CHIME and is designed to detect the BAO feature and constrain the distance-redshift relation.

We demonstrate that self-interacting dark matter models with interactions mediated by light particles can have significant deviations in the matter power-spectrum and detailed structure of galactic halos when compared to a standard cold dark matter scenario. While these deviations can take the form of suppression of small scale structure that are in some ways similar to that of warm dark matter, the self-interacting models have a much wider range of possible phenomenology. A long-range force in the dark matter can introduce multiple scales to the initial power spectrum, in the form of dark acoustic oscillations and an exponential cut-off in the power spectrum. Using simulations we show that the impact of these scales can remain observationally relevant up to the present day. Furthermore, the self-interaction can continue to modify the small-scale structure of the dark matter halos, reducing their central densities and creating a dark matter core. The resulting phenomenology is unique to this type of models.

In this paper we continue to develop the m-mode formalism, a technique for efficient and optimal analysis of wide-field transit radio telescopes, targeted at 21 cm cosmology. We extend this formalism to give an accurate treatment of the polarised sky, fully accounting for the effects of polarisation leakage and cross-polarisation. We use the geometry of the measured set of visibilities to project down to pure temperature modes on the sky, serving as a significant compression, and an effective first filter of polarised contaminants. We use the m-mode formalism with the Karhunen-Loeve transform to give a highly efficient method for foreground cleaning, and demonstrate its success in cleaning realistic polarised skies observed with an instrument suffering from substantial off axis polarisation leakage. We develop an optimal quadratic estimator in the m-mode formalism, which can be efficiently calculated using a Monte-Carlo technique. This is used to assess the implications of foreground removal for power spectrum constraints where we find that our method can clean foregrounds well below the foreground wedge, rendering only scales $k_\parallel < 0.02 h \,\mathrm{Mpc}^{-1}$ inaccessible. As this approach assumes perfect knowledge of the telescope, we perform a conservative test of how essential this is by simulating and analysing datasets with deviations about our assumed telescope. Assuming no other techniques to mitigate bias are applied, we recover unbiased power spectra when the per-feed beam width to be measured to 0.1%, and amplifier gains to be known to 1% within each minute. Finally, as an example application, we extend our forecasts to a wideband 400-800 MHz cosmological observation and consider the implications for probing dark energy, finding a medium-sized cylinder telescope improves the DETF Figure of Merit by around 70% over Planck and Stage II experiments alone.

We investigate the effects of warm dark matter (WDM) on the cosmic 21-cm signal. If dark matter exists as WDM instead of cold dark matter (CDM), its non-negligible velocities can inhibit the formation of low-mass halos that normally form first in CDM models, therefore delaying star-formation. The absence of early sources delays the build-up of UV and X-ray backgrounds that affect the 21-cm radiation signal produced by neutral hydrogen. With use of the 21CMFAST, code, we demonstrate that the pre-reionization 21-cm signal can be changed significantly in WDM models with a free-streaming length equivalent to that of a thermal relic with mass mx of up to ~10-20 keV. In such a WDM cosmology, the 21-cm signal traces the growth of more massive halos, resulting in a delay of the 21-cm absorption signature and followed by accelerated X-ray heating. CDM models where astrophysical sources have a suppressed photon-production efficiency can delay the 21-cm signal as well, although its subsequent evolution is not as rapid as compared to WDM. This motivates using the gradient of the global 21-cm signal to differentiate between some CDM and WDM models. Finally, we show that the degeneracy between the astrophysics and mx can be broken with the 21-cm power spectrum, as WDM models should have a bias-induced excess of power on large scales. This boost in power should be detectable with current interferometers for models with mx < 3 keV, while next generation instruments will easily be able to measure this difference for all relevant WDM models.

If all or a fraction of the dark matter (DM) were coupled to a bath of dark radiation (DR) in the early Universe we expect the combined DM-DR system to give rise to acoustic oscillations of the dark matter until it decouples from the DR. Much like the standard baryon acoustic oscillations, these dark acoustic oscillations (DAO) imprint a characteristic scale, the sound horizon of dark matter, on the matter power spectrum. We compute in detail how the microphysics of the DM-DR interaction affects the clustering of matter in the Universe and show that the DAO physics also gives rise to unique signatures in the temperature and polarization spectra of the cosmic microwave background (CMB). We use cosmological data from the CMB, baryon acoustic oscillations (BAO), and large-scale structure to constrain the possible fraction of interacting DM as well as the strength of its interaction with DR. Like nearly all knowledge we have gleaned about dark matter since inferring its existence this constraint rests on the betrayal by gravity of the location of otherwise invisible dark matter. Although our results can be straightforwardly applied to a broad class of models that couple dark matter particles to various light relativistic species, in order to make quantitative predictions, we model the interacting component as dark atoms coupled to a bath of dark photons. We find that linear cosmological data and CMB lensing put strong constraints on existence of DAO features in the CMB and the large-scale structure of the Universe. Interestingly, we find that at most ~5% of all DM can be very strongly interacting with DR. We show that our results are surprisingly constraining for the recently proposed Double-disk DM model, a novel example of how large-scale precision cosmological data can be used to constrain galactic physics and sub-galactic structure.

The GMRT Epoch of Reionization (EoR) experiment is an ongoing effort to measure the power spectrum from neutral hydrogen at high redshift. We have previously reported an upper limit of (70 mK)^2 at wavenumbers of k=0.65 h/Mpc using a basic piecewise-linear foreground subtraction. In this paper we explore the use of a singular value decomposition to remove foregrounds with fewer assumptions about the foreground structure. Using this method we also quantify, for the first time, the signal loss due to the foreground filter and present new power spectra adjusted for this loss, providing a revised measurement of a 2-sigma upper limit at (248 mK)^2 for k=0.50 h/Mpc. While this revised limit is larger than previously reported, we believe it to be more robust and still represents the best current constraints on reionization at z=8.6.

In this paper we describe the spherical harmonic transit telescope, a novel formalism for the analysis of transit radio telescopes. This all-sky approach bypasses the curved sky complications of traditional interferometry and so is particularly well suited to the analysis of wide-field radio interferometers. It enables compact and computationally efficient representations of the data and its statistics that allow new ways of approaching important problems like map-making and foreground removal. In particular, we show how it enables the use of the Karhunen-Loeve transform as a highly effective foreground filter, suppressing realistic foreground residuals for our fiducial example by at least a factor twenty below the 21cm signal even in highly contaminated regions of the sky. This is despite the presence of the angle-frequency mode mixing inherent in real-world instruments with frequency-dependent beams. We show, using Fisher forecasting, that foreground cleaning has little effect on power spectrum constraints compared to hypothetical foreground-free measurements. Beyond providing a natural real-world data analysis framework for 21cm telescopes now under construction and future experiments, this formalism allows accurate power spectrum forecasts to be made that include the interplay of design constraints and realistic experimental systematics with twenty-first century 21cm science.

The cosmological origin of both dark and baryonic matter can be explained through a unified mechanism called hylogenesis where baryon and antibaryon number are divided between the visible sector and a GeV-scale hidden sector, while the Universe remains net baryon symmetric. The "missing" antibaryons, in the form of exotic hidden states, are the dark matter. We study model-building, cosmological, and phenomenological aspects of this scenario within the framework of supersymmetry, which naturally stabilizes the light hidden sector and electroweak mass scales. Inelastic dark matter scattering on visible matter destroys nucleons, and nucleon decay searches offer a novel avenue for the direct detection of the hidden antibaryonic dark matter sea.

While, to ensure successful cosmology, dark matter (DM) must kinematically decouple from the standard model plasma very early in the history of the Universe, it can remain coupled to a bath of "dark radiation" until a relatively late epoch. One minimal theory that realizes such a scenario is the Atomic Dark Matter model, in which two fermions oppositely charged under a new U(1) dark force are initially coupled to a thermal bath of "dark photons" but eventually recombine into neutral atom-like bound states and begin forming gravitationally-bound structures. As dark atoms have (dark) atom-sized geometric cross sections, this model also provides an example of self-interacting DM with a velocity-dependent cross section. Delayed kinetic decoupling in this scenario predicts novel DM properties on small scales but retains the success of cold DM on larger scales. We calculate the atomic physics necessary to capture the thermal history of this dark sector and show significant improvements over the standard atomic hydrogen calculation are needed. We solve the Boltzmann equations that govern the evolution of cosmological fluctuations in this model and find in detail the impact of the atomic DM scenario on the matter power spectrum and the cosmic microwave background (CMB). This scenario imprints a new length scale, the Dark-Acoustic-Oscillation (DAO) scale, on the matter density field. This DAO scale shapes the small-scale matter power spectrum and determines the minimal DM halo mass at late times which may be many orders of magnitude larger than in a typical WIMP scenario. This model necessarily includes an extra dark radiation component, which may be favoured by current CMB experiments, and we quantify CMB signatures that distinguish an atomic DM scenario from a standard $\Lambda$CDM model containing extra free-streaming particles. [Abridged]

We explore the model-independent constraints from cosmology on a dark-matter particle with no prominent standard model interactions that interacts and thermalizes with other particles in a hidden sector. Without specifying detailed hidden-sector particle physics, we characterize the relevant physics by the annihilation cross section, mass, and temperature ratio of the hidden to visible sectors. While encompassing the standard cold WIMP scenario, we do not require the freeze-out process to be nonrelativistic. Rather, freeze-out may also occur when dark matter particles are semirelativistic or relativistic. We solve the Boltzmann equation to find the conditions that hidden-sector dark matter accounts for the observed dark-matter density, satisfies the Tremaine-Gunn bound on dark-matter phase space density, and has a free-streaming length consistent with cosmological constraints on the matter power spectrum. We show that for masses <1.5 keV no region of parameter space satisfies all these constraints. This is a gravitationally-mediated lower bound on the dark-matter mass for any model in which the primary component of dark matter once had efficient interactions -- even if it has never been in equilibrium with the standard model.

The thermal and expansion history of the Universe before big bang nucleosynthesis is unknown. We investigate the evolution of cosmological perturbations through the transition from an early matter era to radiation domination. We treat reheating as the perturbative decay of an oscillating scalar field into relativistic plasma and cold dark matter. After reheating, we find that subhorizon perturbations in the decay-produced dark matter density are significantly enhanced, while subhorizon radiation perturbations are instead suppressed. If dark matter originates in the radiation bath after reheating, this suppression may be the primary cutoff in the matter power spectrum. Conversely, for dark matter produced nonthermally from scalar decay, enhanced perturbations can drive structure formation during the cosmic dark ages and dramatically increase the abundance of compact substructures. For low reheat temperatures, we find that as much as 50% of all dark matter is in microhalos with M > 0.1 Earth masses at z=100, compared to a fraction of 1e-10 in the standard case. In this scenario, ultradense substructures may constitute a large fraction of dark matter in galaxies today.

Current theories of the origin of the Universe, including string theory, predict the existence of a multiverse containing many bubble universes. These bubble universes will generically collide, and collisions with ours produce cosmic wakes that enter our Hubble volume, appear as unusually symmetric disks in the cosmic microwave background (CMB) and disturb large scale structure (LSS). There is preliminary observational evidence consistent with one or more of these disturbances on our sky. However, other sources can produce similar features in the CMB temperature map and so additional signals are needed to verify their extra-universal origin. Here we find, for the first time, the detailed three-dimensional shape and CMB temperature and polarization signals of the cosmic wake of a bubble collision in the early universe consistent with current observations. The predicted polarization pattern has distinctive features that when correlated with the corresponding temperature pattern are a unique and striking signal of a bubble collision. These features represent the first verifiable prediction of the multiverse paradigm and might be detected by current experiments such as Planck and future CMB polarization missions. A detection of a bubble collision would confirm the existence of the Multiverse, provide compelling evidence for the string theory landscape, and sharpen our picture of the Universe and its origins.

We investigate new and unusual signals that arise in theories where dark matter is asymmetric and carries a net antibaryon number, as may occur when the dark matter abundance is linked to the baryon abundance. Antibaryonic dark matter can cause {\it induced nucleon decay} by annihilating visible baryons through inelastic scattering. These processes lead to an effective nucleon lifetime of 10^{29}-10^{32} years in terrestrial nucleon decay experiments, if baryon number transfer between visible and dark sectors arises through new physics at the weak scale. The possibility of induced nucleon decay motivates a novel approach for direct detection of cosmic dark matter in nucleon decay experiments. Monojet searches (and related signatures) at hadron colliders also provide a complementary probe of weak-scale dark-matter--induced baryon number violation. Finally, we discuss the effects of baryon-destroying dark matter on stellar systems and show that it can be consistent with existing observations.

Prior to recombination photons, electrons, and atomic nuclei rapidly scattered and behaved, almost, like a single tightly-coupled photon-baryon plasma. We investigate here the accuracy of the tight-coupling approximation commonly used to numerically evolve the baryon and photon perturbation equations at early times. By solving the exact perturbations equations with a stiff solver starting deep in the radiation-dominated epoch we find the level of inaccuracy introduced by resorting to the standard first-order tight-coupling approximation. We develop a new second-order approximation in the inverse Thomson opacity expansion and show that it closely tracks the full solution, at essentially no extra numerical cost. We find the bias on estimates of cosmological parameters introduced by the first-order approximation is, for most parameters, negligible. Finally, we show that our second-order approximation can be used to reduce the time needed to compute cosmic microwave background angular spectra by as much as ~17%.

We predict the polarization of cosmic microwave background (CMB) photons that results from a cosmic bubble collision. The polarization is purely E-mode, symmetric around the axis pointing towards the collision bubble, and has several salient features in its radial dependence that can help distinguish it from a more conventional explanation for unusually cold or hot features in the CMB sky. The anomalous "cold spot" detected by the Wilkinson Microwave Anisotropy Probe (WMAP) satellite is a candidate for a feature produced by such a collision, and the Planck satellite and other proposed surveys will measure the polarization on it in the near future. The detection of such a collision would provide compelling evidence for the string theory landscape.

We present a new upper limit to the 21cm power spectrum during the Epoch of Reionization (EoR) which constrains reionization models with an unheated IGM. The GMRT-EoR experiment is an ongoing effort to make a statistical detection of the power spectrum of 21cm neutral hydrogen emission at redshift z~9. Data from this redshift constrain models of the EoR, the end of the Dark Ages arising from the formation of the first bright UV sources, probably stars or mini-quasars. We present results from approximately 50 hours of observations at the Giant Metrewave Radio Telescope in India from December 2007. We describe radio frequency interference (RFI) localisation schemes which allow bright sources on the ground to be identified and physically removed in addition to automated flagging. Singular-value decomposition is used to remove remaining broadband RFI by identifying ground sources with large eigenvalues. Foregrounds are modelled using a piecewise linear filter and the power spectrum is measured using cross-correlations of foreground subtracted images.

We present a novel mechanism for generating both the baryon and dark matter densities of the Universe. A new Dirac fermion X carrying a conserved baryon number charge couples to the Standard Model quarks as well as a GeV-scale hidden sector. CP-violating decays of X, produced non-thermally in low-temperature reheating, sequester antibaryon number in the hidden sector, thereby leaving a baryon excess in the visible sector. The antibaryonic hidden states are stable dark matter. A spectacular signature of this mechanism is the baryon-destroying inelastic scattering of dark matter that can annihilate baryons at appreciable rates relevant for nucleon decay searches.

We demonstrate a new model which uses an ADD type braneworld scenario to produce a multi-state theory of dark matter. Compactification of the extra dimensions onto a sphere leads to the association of a single complex scalar in the bulk with multiple Kaluza-Klein towers in an effective four-dimensional theory. A mutually interacting multi-state theory of dark matter arises naturally within which the dark matter states are identified with the lightest Kaluza-Klein particles of fixed magnetic quantum number. These states are protected from decay by a combination of a global U(1) symmetry and the continuous rotational symmetry about the polar axis of the spherical geometry. We briefly discuss the relic abundance calculation and investigate the spin-independent elastic scattering off nucleons of the lightest and next-to-lightest dark matter states.

Recent claims in the literature have suggested that the {\it WMAP} quadrupole is not primordial in origin, and arises from an aliasing of the much larger dipole field because of incorrect satellite pointing. We attempt to reproduce this result and delineate the key physics leading to the effect. We find that, even if real, the induced quadrupole would be smaller than claimed. We discuss reasons why the {\it WMAP} data are unlikely to suffer from this particular systematic effect, including the implications for observations of point sources. Given this evidence against the reality of the effect, the similarity between the pointing-offset-induced signal and the actual quadrupole then appears to be quite puzzling. However, we find that the effect arises from a convolution between the gradient of the dipole field and anisotropic coverage of the scan direction at each pixel. There is something of a directional conspiracy here -- the dipole signal lies close to the Ecliptic Plane, and its direction, together with the {\it WMAP} scan strategy, results in a strong coupling to the $Y_{2,\,-1}$ component in Ecliptic co-ordinates. The dominant strength of this component in the measured quadrupole suggests that one should exercise increased caution in interpreting its estimated amplitude. The {\it Planck} satellite has a different scan strategy which does not so directly couple the dipole and quadrupole in this way and will soon provide an independent measurement.

We address the question of whether the upcoming generation of dark matter search experiments and colliders will be able to discover if the dark matter in the Universe has two components of weakly interacting massive particles (WIMPs). We outline a model-independent approach, and we study the specific cases of (1) direct detection with low-background 1 ton noble-gas detectors and (2) a 0.5 TeV center of mass energy electron-positron linear collider. We also analyze the case of indirect detection via two gamma-ray lines, which would provide a verification of such a discovery, although multiple gamma-ray lines can in principle originate from the annihilation of a single dark matter particle. For each search "channel", we outline a few assumptions to relate the very small set of parameters we consider (defining the masses of the two WIMPs and their relative abundance in the overall dark matter density) with the relevant detection rates. We then draw general conclusions on which corners of a generic dual-component dark matter scenario can be explored with current and next generation experiments. We find that in all channels the ideal setup is one where the relative mass splitting between the two WIMP species is of order 1, and where the two dark matter components contribute in a ratio close to 1:1 to the overall dark matter content of the Universe. Interestingly, in the case of direct detection, future experiments might detect multiple states even if only ~ 10% of the energy-density of dark matter in the Universe is in the subdominant species.

We show that hidden hot dark matter, hidden-sector dark matter with interactions that decouple when it is relativistic, is a viable dark matter candidate provided it has never been in thermal equilibrium with the particles of the standard model. This hidden hot dark matter may reheat to a lower temperature and number density than the visible Universe and thus account, simply with its thermal abundance, for all the dark matter in the Universe while evading the typical constraints on hot dark matter arising from structure formation. We find masses ranging from ~3 keV to ~10 TeV. While never in equilibrium with the standard model, this class of models may have unique observational signatures in the matter power spectrum or via extra-weak interactions with standard model particles.

These lectures are intended to provide a brief pedagogical review of dark matter for the newcomer to the subject. We begin with a discussion of the astrophysical evidence for dark matter. The standard weakly-interacting massive particle (WIMP) scenario--the motivation, particle models, and detection techniques--is then reviewed. We provide a brief sampling of some recent variations to the standard WIMP scenario as well as some alternatives (axions and sterile neutrinos). Exercises are provided for the reader.

We compute the probability distribution for bubble collisions in an inflating false vacuum which decays by bubble nucleation. Our analysis generalizes previous work of Guth, Garriga, and Vilenkin to the case of general cosmological evolution inside the bubble, and takes into account the dynamics of the domain walls that form between the colliding bubbles. We find that incorporating these effects changes the results dramatically: the total expected number of bubble collisions in the past lightcone of a typical observer is N ~ \gamma V_f / V_i, where \gamma is the fastest decay rate of the false vacuum, V_f is its vacuum energy, and V_i is the vacuum energy during inflation inside the bubble. This number can be large in realistic models without tuning. In addition, we calculate the angular position and size distribution of the collisions on the cosmic microwave background sky, and demonstrate that the number of bubbles of observable angular size is N_{LS} \sim \sqrt{\Omega_k} N, where \Omega_k is the curvature contribution to the total density at the time of observation. The distribution is almost exactly isotropic.

It has recently been argued that the PAMELA, ATIC and PPB-BETS data showing an anomalous excess of high-energy cosmic ray positrons and electrons might be explained by dark matter annihilating in the Galactic halo with a cross section resonantly enhanced compared to its value in the primeval plasma. We find that with a very large annihilation cross section the flash of energetic photons and electron-positron pairs expected from dark-matter annihilation in the first protohalos that form at redshift z~40 is likely substantial and observable. As a consequence, bounds on the allowed energy injection into the primordial gas and the energy density of the diffuse gamma-ray background give rise to limits on the low-velocity dark matter cross section that can be difficult to reconcile with this interpretation of the PAMELA, ATIC and PPB-BETS results.

The GMRT reionization effort aims to map out the large scale structure of the Universe during the epoch of reionization (EoR). Removal of polarized Galactic emission is a difficult part of any 21 cm EoR program, and we present new upper limits to diffuse polarized foregrounds at 150 MHz. We find no high significance evidence of polarized emission in our observed field at mid galactic latitude (J2000 08h26m+26). We find an upper limit on the 2-dimensional angular power spectrum of diffuse polarized foregrounds of [l^2 C_l/(2 PI)]^{1/2}< 3K in frequency bins of width 1 MHz at 300<l<1000. The 3-dimensional power spectrum of polarized emission, which is most directly relevant to EoR observations, is [k^3 P_p(k)/(2 PI^2)]^{1/2}< 2K at k_perp > 0.03 h/Mpc, k < 0.1 h/Mpc. This can be compared to the expected EoR signal in total intensity of [k^3 P(k)/ (2 PI^2) ]^{1/2} ~ 10 mK. We find polarized structure is substantially weaker than suggested by extrapolation from higher frequency observations, so the new low upper limits reported here reduce the anticipated impact of these foregrounds on EoR experiments. We discuss Faraday beam and depth depolarization models and compare predictions of these models to our data. We report on a new technique for polarization calibration using pulsars, as well as a new technique to remove broadband radio frequency interference. Our data indicate that, on the edges of the main beam at GMRT, polarization squint creates ~ 3% leakage of unpolarized power into polarized maps at zero rotation measure. Ionospheric rotation was largely stable during these solar minimum night time observations.

How did the universe evolve? The fine angular scale (l>1000) temperature and polarization anisotropies in the CMB are a Rosetta stone for understanding the evolution of the universe. Through detailed measurements one may address everything from the physics of the birth of the universe to the history of star formation and the process by which galaxies formed. One may in addition track the evolution of the dark energy and discover the net neutrino mass. We are at the dawn of a new era in which hundreds of square degrees of sky can be mapped with arcminute resolution and sensitivities measured in microKelvin. Acquiring these data requires the use of special purpose telescopes such as the Atacama Cosmology Telescope (ACT), located in Chile, and the South Pole Telescope (SPT). These new telescopes are outfitted with a new generation of custom mm-wave kilo-pixel arrays. Additional instruments are in the planning stages.

Using the 21 cm line, observed all-sky and across the redshift range from 0 to 5, the large scale structure of the Universe can be mapped in three dimensions. This can be accomplished by studying specific intensity with resolution ~ 10 Mpc, rather than via the usual galaxy redshift survey. The data set can be analyzed to determine Baryon Acoustic Oscillation wavelengths, in order to address the question: 'What is the nature of Dark Energy?' In addition, the study of Large Scale Structure across this range addresses the questions: 'How does Gravity effect very large objects?' and 'What is the composition our Universe?' The same data set can be used to search for and catalog time variable and transient radio sources.

Fluctuations in high-redshift cosmic 21-cm radiation provide a new window for observing unconventional effects of high-energy physics in the primordial spectrum of density perturbations. In scenarios for which the initial state prior to inflation is modified at short distances, or for which deviations from scale invariance arise during the course of inflation, the cosmic 21-cm power spectrum can in principle provide more precise measurements of exotic effects on fundamentally different scales than corresponding observations of cosmic microwave background anisotropies.

We carry out a model independent study of resonant photon scattering off dark matter (DM) particles. The DM particle chi_1 can feature an electric or magnetic transition dipole moment which couples it with photons and a heavier neutral particle chi_2. Resonant photon scattering then takes place at a special energy set by the masses of chi_1 and chi_2, with the width of the resonance set by the size of the transition dipole moment. We compute the constraints on the parameter space of the model from stellar energy losses, data from SN 1987A, the Lyman-alpha forest, Big Bang nucleosynthesis, electro-weak precision measurements and accelerator searches. We show that the velocity broadening of the resonance plays an essential role for the possibility of the detection of a spectral feature originating from resonant photon-DM scattering. Depending upon the particle setup and the DM surface mass density, the favored range of DM particle masses lies between tens of keV and a few MeV, while the resonant photon absorption energy is predicted to be between tens of keV and few GeV.

The 21-cm hyperfine line of atomic hydrogen (HI) is a promising probe of the cosmic dark ages. In past treatments of 21-cm radiation it was assumed the hyperfine level populations of HI could be characterized by a velocity-independent ``spin temperature'' T_s determined by a competition between 21-cm radiative transitions, spin-changing collisions, and (at lower redshifts) Lyman-alpha scattering. However we show here that, if the collisional time is comparable to the radiative time, the spin temperature will depend on atomic velocity, T_s=T_s(v), and one must replace the usual hyperfine level rate equations with a Boltzmann equation describing the spin and velocity dependence of the HI distribution function. We construct here the Boltzmann equation relevant to the cosmic dark ages and solve it using a basis-function method. Accounting for the actual spin-resolved atomic velocity distribution results in up to a 2 per cent suppression of the 21-cm emissivity, and a redshift and angular-projection dependent suppression or enhancement of the linear power spectrum of 21-cm fluctuations of up to 5 per cent. The effect on the 21-cm line profile is more dramatic -- its full-width at half maximum (FWHM) can be enhanced by up to 60 per cent relative to the velocity-independent calculation. We discuss the implications for 21-cm tomography of the dark ages.

We investigate signatures of neutrino scattering in the Cosmic Microwave Background (CMB) and matter power spectra, and the extent to which present cosmological data can distinguish between a free streaming or tightly coupled fluid of neutrinos. If neutrinos have strong non-standard interactions, for example, through the coupling of neutrinos to a light boson, they may be kept in equilibrium until late times. We show how the power spectra for these models differ from more conventional neutrino scenarios, and use CMB and large scale structure data to constrain these models. CMB polarization data improves the constraints on the number of massless neutrinos, while the Lyman--$\alpha$ power spectrum improves the limits on the neutrino mass. Neutrino mass limits depend strongly on whether some or all of the neutrino species interact and annihilate. The present data can accommodate a number of tightly-coupled relativistic degrees of freedom, and none of the interacting-neutrino scenarios considered are ruled out by current data -- although considerations regarding the age of the Universe disfavor a model with three annihilating neutrinos with very large neutrino masses.

We calculate the kinetic-decoupling temperature for weakly interacting massive particles (WIMPs) in supersymmetric (SUSY) and universal-extra-dimension (UED) models that can account for the cold-dark-matter abundance determined from cosmic microwave background measurements. Depending on the parameters of the particle-physics model, a wide variety of decoupling temperatures is possible, ranging from several MeV to a few GeV. These decoupling temperatures imply a range of masses for the smallest protohalos much larger than previously thought -- ranging from 10^{-6} earth masses to 10^{2} earth masses. We expect the range of protohalos masses derived here to be characteristic of most particle-physics models that can thermally accommodate the required relic abundance of WIMP dark matter, even beyond SUSY and UED.

We discuss how measurements of fluctuations in the absorption of cosmic microwave background (CMB) photons by neutral gas during the cosmic dark ages, at redshifts z ~ 7--200, could reveal the primordial deuterium abundance of the Universe. The strength of the cross-correlation of brightness-temperature fluctuations due to resonant absorption of CMB photons in the 21-cm line of neutral hydrogen with those due to resonant absorption of CMB photons in the 92-cm line of neutral deuterium is proportional to the fossil deuterium to hydrogen ratio [D/H] fixed during big bang nucleosynthesis (BBN). Although technically challenging, this measurement could provide the cleanest possible determination of [D/H], free from contamination by structure formation processes at lower redshifts, and has the potential to improve BBN constraints to the baryon density of the Universe \Omega_{b} h^2. We also present our results for the thermal spin-change cross-section for deuterium-hydrogen scattering, which may be useful in a more general context than we describe here.

The curl (B) modes of cosmic microwave background (CMB) polarization anisotropies are a unique probe of the primordial background of inflationary gravitational waves (IGWs). Unfortunately, the B-mode polarization anisotropies generated by gravitational waves at recombination are confused with those generated by the mixing of gradient-mode (E-mode) and B-mode polarization anisotropies as CMB photons propagate through the Universe and are gravitationally lensed. We describe here a method for delensing CMB polarization anisotropies using observations of anisotropies in the cosmic 21-cm radiation emitted or absorbed by neutral hydrogen atoms at redshifts 10 to 200. While the detection of cosmic 21-cm anisotropies at high resolution is challenging, a combined study with a relatively low-resolution (but high-sensitivity) CMB polarization experiment could probe inflationary energy scales well below the Grand Unified Theory (GUT) scale of 10^{16} GeV -- constraining models with energy scales below 10^{15} GeV (the detectable limit derived from CMB observations alone). The ultimate theoretical limit to the detectable inflationary energy scale via this method may be as low as 3 \times 10^{14} GeV.

We review the physical and cosmological consequences of two possible electromagnetic couplings to the dark sector: (i) a neutral lightest dark-matter particle (LDP) with nonzero electric and/or magnetic dipole moments and (ii) a charged next-to-lightest dark-matter particle (NLDP) which decays to a neutral LDP. For scenario (i) we find that a relatively light particle with mass between a few MeV and a few GeV and an electric or magnetic dipole as large as $\sim 3\times 10^{-16}e$ cm (roughly $1.6\times10^{-5} \mu_B$) satisfies experimental and observational bounds. In scenario (ii), we show that charged-particles decaying in the early Universe result in a suppression of the small-scale matter power spectrum on scales that enter the horizon prior to decay. This leads to either a cutoff in the matter power spectrum, or if the charged fraction is less than unity, an effect in the power spectrum that might resemble a running (scale-dependent) spectral index in small-scale data.

We show that charged-particles decaying in the early Universe can induce a scale-dependent or `running' spectral index in the small-scale linear and nonlinear matter power spectrum and discuss examples of this effect in minimal supersymmetric models in which the lightest neutralino is a viable cold-dark-matter candidate. We find configurations in which the neutralino relic density is set by coannihilations with a long-lived stau, and the late decay of staus partially suppresses the linear matter power spectrum. Nonlinear evolution on small scales then causes the modified linear power spectrum to evolve to a nonlinear power spectrum similar (but different in detail) to models parametrized by a constant running $\alpha_{s}=d n_{s}/d{\rm ln} k$ by redshifts of 2 to 4. Thus, Lyman-$\alpha$ forest observations, which probe the matter power spectrum at these redshifts, might not discriminate between the two effects. However, a measurement of the angular power spectrum of primordial 21-cm radiation from redshift $z \approx 30$--$200$ might distinguish between this charged-decay model and a primordial running spectral index. The direct production of a long-lived charged particle at future colliders is a dramatic prediction of this model.

We review several aspects of the Sunyaev-Zel'dovich (SZ) effect associated with the large scale baryon distribution and its characteristic signatures in the statistics of cosmic microwave background (CMB) anisotropies. We discuss (1) the contributions to the angular power spectrum from the thermal and kinetic SZ effects, (2) the effect of SZ non-Gaussianities on cosmological parameter estimation, and (3) the SZ-induced CMB polarization towards galaxy clusters. We also discuss the extent to which recent measurements of small angular scale CMB anisotropy can be accounted for by SZ clusters and the potential use of SZ polarization towards large samples of galaxy clusters as a probe of dark energy parameters.

We consider the consequences of a neutral dark-matter particle with a nonzero electric and/or magnetic dipole moment. Theoretical constraints, as well as constraints from direct searches, precision tests of the standard model, the cosmic microwave background and matter power spectra, and cosmic gamma rays, are included. We find that a relatively light particle with mass between an MeV and a few GeV and an electric or magnetic dipole as large as $\sim 3\times 10^{-16}e$ cm (roughly $1.6\times10^{-5} \mu_B$) satisfies all experimental and observational bounds. Some of the remaining parameter space may be probed with forthcoming more sensitive direct searches and with the Gamma-Ray Large Area Space Telescope.

A Bayesian likelihood analysis of the stellar density profile shapes of several dark matter dominated Galactic dwarf spheroidal galaxies (dSphs) has been completed. This analysis indicates that the inner power law, \rho \propto r^{-\gamma}, of each of these profiles satisfies \gamma < 1.21 at greater than 99% confidence with a slight preference for \gamma \sim 0.7 . Speculating that this bound on the stellar profiles may be related to the inner power law of the underlying dark matter distribution reveals a statistical discrepancy with the values favoured by recent N-body simulations of cold dark matter (CDM), \gamma_{CDM} \sim 1.4-1.5 .

We study the suppression of the small-scale power spectrum due to the decay of charged matter to dark matter prior to recombination. Prior to decay, the charged particles couple to the photon-baryon fluid and participate in its acoustic oscillations. However, after decaying to neutral dark matter the photon-baryon fluid is coupled only gravitationally to the newly-created dark matter. This generically leads to suppression of power on length scales that enter the horizon prior to decay. For decay times of approximately 3.5 years this leads to suppression of power on subgalactic scales, bringing the observed number of Galactic substructures in line with observation. Decay times of a few years are possible if the dark matter is purely gravitationally interacting, such as the gravitino in supersymmetric models or a massive Kaluza-Klein graviton in models with universal extra dimensions.

We study the effects on cosmic microwave background (CMB) temperature and polarization anisotropies of spatial fluctuations of the fine-structure parameter between causally disconnected regions of the Universe at the time of recombination. Analogous to weak gravitational lensing, in addition to modifying the mean power spectra and inducing a curl component (B mode) to the polarization, spatial fluctuations of the fine-structure parameter induce higher-order (non-Gaussian) temperature and polarization correlations in the CMB. We calculate these effects for the general case of arbitrary correlation between temperature fluctuations and fine-structure parameter fluctuations, and show the results for a model where these two types of fluctuations are uncorrelated. The formalism we present here may also be applied to other modifications of recombination physics that do not significantly alter the evolution of the dominant density perturbations. We discuss the constraints on the effective Lagrangian for variable fine-structure parameter necessary to realize this scenario.

We show how estimates of parameters characterizing inflation-based theories of structure formation localized over the past year when large scale structure (LSS) information from galaxy and cluster surveys was combined with the rapidly developing cosmic microwave background (CMB) data, especially from the recent Boomerang and Maxima balloon experiments. All current CMB data plus a relatively weak prior probability on the Hubble constant, age and LSS points to little mean curvature (Omega_{tot} = 1.08\pm 0.06) and nearly scale invariant initial fluctuations (n_s =1.03\pm 0.08), both predictions of (non-baroque) inflation theory. We emphasize the role that degeneracy among parameters in the L_{pk} = 212\pm 7 position of the (first acoustic) peak plays in defining the $\Omega_{tot}$ range upon marginalization over other variables. Though the CDM density is in the expected range (\Omega_{cdm}h^2=0.17\pm 0.02), the baryon density Omega_bh^2=0.030\pm 0.005 is somewhat above the independent 0.019\pm 0.002 nucleosynthesis estimate. CMB+LSS gives independent evidence for dark energy (Omega_\Lambda=0.66\pm 0.06) at the same level as from supernova (SN1) observations, with a phenomenological quintessence equation of state limited by SN1+CMB+LSS to w_Q<-0.7 cf. the w_Q=-1 cosmological constant case.

The past, present and future of cosmic microwave background (CMB) anisotropy research is discussed, with emphasis on the Boomerang and Maxima balloon experiments. These data are combined with large scale structure (LSS) information and high redshift supernova (SN1) observations to explore the inflation-based cosmic structure formation paradigm. Here we primarily focus on a simplified inflation parameter set, {omega_b,omega_{cdm},Omega_{tot}, Omega_Q,w_Q, n_s,tau_C, sigma_8}. After marginalizing over the other cosmic and experimental variables, we find the current CMB+LSS+SN1 data gives Omega_{tot}=1.04\pm 0.05, consistent with (non-baroque) inflation theory. Restricting to Omega_{tot}=1, we find a nearly scale invariant spectrum, n_s =1.03 \pm 0.07. The CDM density, omega_{cdm}=0.17\pm 0.02, is in the expected range, but the baryon density, omega_b=0.030\pm 0.004, is slightly larger than the current nucleosynthesis estimate. Substantial dark energy is inferred, Omega_Q\approx 0.68\pm 0.05, and CMB+LSS Omega_Q values are compatible with the independent SN1 estimates. The dark energy equation of state, parameterized by a quintessence-field pressure-to-density ratio w_Q, is not well determined by CMB+LSS (w_Q<-0.3 at 95%CL), but when combined with SN1 the resulting w_Q<-0.7 limit is quite consistent with the w_Q=-1 cosmological constant case. Though forecasts of statistical errors on parameters for current and future experiments are rosy, rooting out systematic errors will define the true progress.

The confrontation of the Cosmic Microwave Background (CMB) theoretical angular power spectrum with available data often requires the calculation of large numbers of power spectra. The standard practice is to use a fast code to compute the CMB power spectra over some large parameter space, in order to estimate likelihoods and constrain these parameters. But as the dimensionality of the space under study increases, then even with relatively fast anisotropy codes, the computation can become prohibitive. This paper describes the employment of a "morphing" strategy to interpolate new power spectra based on previously calculated ones. We simply present the basic idea here, and illustrate with a few examples; optimization of interpolation schemes will depend on the specific application. In addition to facilitating the exploration of large parameter spaces, this morphing technique may be helpful for Fisher matrix calculations involving derivatives.