Understanding the non-linear dynamics of satellite halos (a.k.a. ``sub-halos'') is important for predicting the abundance and distribution of dark matter substructures and satellite galaxies, and for distinguishing among microphysical dark matter models using observations. Typically, modeling these dynamics requires large N-body simulations with high resolution. Semi-analytic models can provide a more efficient way to describe the key physical processes such as dynamical friction, tidal mass loss, and tidal heating, with only a few free parameters. In this work, we present a fast Monte Carlo Markov Chain fitting approach to explore the parameter space of such a sub-halo non-linear evolution model. We use the dynamical models described in an earlier work and calibrate the models to two sets of high-resolution cold dark matter N-body simulations, ELVIS and Caterpillar. Compared to previous calibrations that used manual parameter tuning, our approach provides a more robust way to determine the best-fit parameters and their posterior probabilities. We find that jointly fitting for the sub-halo mass and maximum velocity functions can break the degeneracy between tidal stripping and tidal heating parameters, as well as providing better constraints on the strength of dynamical friction. We show that our semi-analytic simulation can accurately reproduce N-body simulations statistics, and that the calibration results for the two sets of N-body simulations agree at 95% confidence level. Dynamical models calibrated in this work will be important for future dark matter substructure studies.

Circular polarization is a relatively unexplored realm of CMB physics. Given the substantial community effort towards building next generation CMB polarization experiments, including those which will be sensitive to circular polarization, it behooves theorists to understand the possible sources and relevant physics of circular polarization, as encoded in the Stokes V parameter. In this work we develop and derive the requisite formalism, namely the Boltzmann hierarchy for V-mode scalar, vector, and tensor, anisotropies. We derive the V-mode anisotropies induced by a general source term, and demonstrate how existing proposals for the generation of V can be incorporated as source terms in the Boltzmann hierarchy. A subset of these effects may be correlated with 21cm intensity; we provide a worked example and derive an estimator to extract this information from observations. We conclude by computing the CMB TV cross-correlation generated by axions, and find a relation between TV and VV spectra in axion models.

The EXperiment for Cryogenic Large-Aperture Intensity Mapping (EXCLAIM) is a cryogenic balloon-borne instrument that will survey galaxy and star formation history over cosmological time scales. Rather than identifying individual objects, EXCLAIM will be a pathfinder to demonstrate an intensity mapping approach, which measures the cumulative redshifted line emission. EXCLAIM will operate at 420-540 GHz with a spectral resolution R=512 to measure the integrated CO and [CII] in redshift windows spanning 0 < z < 3.5. CO and [CII] line emissions are key tracers of the gas phases in the interstellar medium involved in star-formation processes. EXCLAIM will shed light on questions such as why the star formation rate declines at z < 2, despite continued clustering of the dark matter. The instrument will employ an array of six superconducting integrated grating-analog spectrometers (micro-spec) coupled to microwave kinetic inductance detectors (MKIDs). Here we present an overview of the EXCLAIM instrument design and status.

In the cold dark matter (CDM) picture of structure formation, galaxy mass distributions are predicted to have a considerable amount of structure on small scales. Strong gravitational lensing has proven to be a useful tool for studying this small-scale structure. Much of the attention has been given to detecting individual dark matter subhalos through lens modeling, but recent work has suggested that the full population of subhalos could be probed using a power spectrum analysis. In this paper we quantify the power spectrum of small-scale structure in simulated galaxies, with the goal of understanding theoretical predictions and setting the stage for using measurements of the power spectrum to test dark matter models. We use a sample of simulated galaxies generated from the \texttt{Galacticus} semi-analytic model to determine the power spectrum distribution first in the CDM paradigm and then in a warm dark matter scenario. We find that a measurement of the slope and amplitude of the power spectrum on galaxy strong lensing scales ($k\sim 1$ kpc$^{-1}$) could be used to distinguish between CDM and alternate dark matter models, especially if the most massive subhalos can be directly detected via gravitational imaging.

MegaMapper is a proposed ground-based experiment to measure Inflation parameters and Dark Energy from galaxy redshifts at 2<z<5. A 6.5-m Magellan telescope will be coupled with DESI spectrographs to achieve multiplexing of 20,000. MegaMapper would be located at Las Campanas Observatory to fully access LSST imaging for target selection.

In response to the need for the Astro2020 Decadal Survey to explicitly engage early career astronomers, the National Academies of Sciences, Engineering, and Medicine hosted the Early Career Astronomer and Astrophysicist Focus Session (ECFS) on October 8-9, 2018 under the auspices of Committee of Astronomy and Astrophysics. The meeting was attended by fifty six pre-tenure faculty, research scientists, postdoctoral scholars, and senior graduate students, as well as eight former decadal survey committee members, who acted as facilitators. The event was designed to educate early career astronomers about the decadal survey process, to solicit their feedback on the role that early career astronomers should play in Astro2020, and to provide a forum for the discussion of a wide range of topics regarding the astrophysics career path. This white paper presents highlights and themes that emerged during two days of discussion. In Section 1, we discuss concerns that emerged regarding the coming decade and the astrophysics career path, as well as specific recommendations from participants regarding how to address them. We have organized these concerns and suggestions into five broad themes. These include (sequentially): (1) adequately training astronomers in the statistical and computational techniques necessary in an era of "big data", (2) responses to the growth of collaborations and telescopes, (3) concerns about the adequacy of graduate and postdoctoral training, (4) the need for improvements in equity and inclusion in astronomy, and (5) smoothing and facilitating transitions between early career stages. Section 2 is focused on ideas regarding the decadal survey itself, including: incorporating early career voices, ensuring diverse input from a variety of stakeholders, and successfully and broadly disseminating the results of the survey.

Weak gravitational lensing of the 21 cm radiation is expected to be an important cosmological probe for post-reionization physics. We investigate the reconstruction of the matter density perturbations using a quadratic minimum variance estimator. The next generation of line intensity mapping (LIM) surveys such as HIRAX and CHIME will cover a larger sky fraction, which requires one to account for the curvature in the sky. Thus, we extend the plane-parallel flat-sky formalism for lensing reconstruction to account for a full-sky survey using the Spherical Fourier-Bessel (SFB) expansion. Using the HIRAX 21 cm survey as a basis, we make predictions for lensing-reconstruction noise in our formalism and compare our results with the predictions from the plane-parallel formalism. We find agreement with the plane-parallel noise power spectrum at small scales and a significant deviation at scales $L\lesssim \ell_{\rm res}-k_{\rm eq}R$ where $R$ is the radius of the shell volume, $k_{\rm eq}$ is the wavenumber for matter-radiation equality, and $\ell_{\rm res}$ is the angular resolution scale. Furthermore, we derive the SFB flat-sky reconstruction noise and compare it with the full-sky SFB case as well as the plane-parallel case, finding minor deviations from the full-sky noise due to sphericity. We also determine that, in the absence of non-Gaussian statistics of the intensity field but accounting for foregrounds, the signal-to-noise ratio (SNR) for $C_\ell^{\phi\phi}$ using our SFB estimator increases by 107%. This shows that accounting for the curved sky in LIM weak lensing will be crucial for large-scale cosmology.

CII is one of the brightest emission lines from star-forming galaxies and is an excellent tracer for star formation. Recent work measured the CII emission line amplitude for redshifts $2<z<3.2$ by cross-correlating Planck High Frequency Instrument emission maps with tracers of overdensity from the Baryon Oscillation Spectroscopic Sky Survey, finding $I_{\mathrm{CII}}=6.6^{+5.0}_{-4.8}\times 10^4$ Jy/sr at $95\%$ confidence. In this paper, we present a refinement of this earlier work by improving the mask weighting in each of the Planck bands and the precision in the covariance matrix. We report a detection of excess emission in the 545 GHz Planck band separate from the cosmic infrared background (CIB) present in the 353-857 GHz Planck bands. This excess is consistent with redshifted CII emission, in which case we report $b_{\mathrm{CII}}I_{\mathrm{CII}} =2.0^{+1.2}_{-1.1}\times 10^5$ Jy/sr at $95\%$ confidence, which strongly favors many collisional excitation models of CII emission. Our detection shows strong evidence for a model with a non-zero CII parameter, though line intensity mapping observations at high spectral resolution will be needed to confirm this result.

Line-Intensity Mapping is an emerging technique which promises new insights into the evolution of the Universe, from star formation at low redshifts to the epoch of reionization and cosmic dawn. It measures the integrated emission of atomic and molecular spectral lines from galaxies and the intergalactic medium over a broad range of frequencies, using instruments with aperture requirements that are greatly relaxed relative to surveys for single objects. A coordinated, comprehensive, multi-line intensity-mapping experimental effort can efficiently probe over 80% of the volume of the observable Universe - a feat beyond the reach of other methods. Line-intensity mapping will uniquely address a wide array of pressing mysteries in galaxy evolution, cosmology, and fundamental physics. Among them are the cosmic history of star formation and galaxy evolution, the compositions of the interstellar and intergalactic media, the physical processes that take place during the epoch of reionization, cosmological inflation, the validity of Einstein's gravity theory on the largest scales, the nature of dark energy and the origin of dark matter.

Intensity mapping has attracted significant interest as an approach to measure the properties of the interstellar medium in typical galaxies at high redshift. Intensity mapping measures the statistics of surface brightness as a function of frequency, making it sensitive not only to all line emission of interest but also radiation from all other sources. Significant effort has gone into developing approaches that reject foreground contamination. Additionally, the target galaxies have multiple sources of emission that can complicate the interpretation of the line brightness. We describe the problem of jointly estimating correlated continuum emission and cleaning uncorrelated continuum emission, such as from the Milky Way. We apply these considerations to a cross-correlation of Planck data with BOSS quasars for a determination of CII for 2 < z < 3.2. Intensity mapping surveys with few bands have unique challenges for treating foregrounds and avoiding bias from correlated continuum emission. We show how a future intensity mapping survey with many bands can separate line from continuum emission in cross-correlation.

We present a search for CII emission over cosmological scales at high-redshifts. The CII line is a prime candidate to be a tracer of star formation over large-scale structure since it is one of the brightest emission lines from galaxies. Redshifted CII emission appears in the submillimeter regime, meaning it could potentially be present in the higher frequency intensity data from the Planck satellite used to measure the cosmic infrared background (CIB). We search for CII emission over redshifts z=2-3.2 in the Planck 545 GHz intensity map by cross-correlating the 3 highest frequency Planck maps with spectroscopic quasars and CMASS galaxies from the Sloan Digital Sky Survey III (SDSS-III), which we then use to jointly fit for CII intensity, CIB parameters, and thermal Sunyaev-Zeldovich (SZ) emission. We report a measurement of an anomalous emission $\mathrm{I_\nu}=6.6^{+5.0}_{-4.8}\times10^4$ Jy/sr at 95% confidence, which could be explained by CII emission, favoring collisional excitation models of CII emission that tend to be more optimistic than models based on CII luminosity scaling relations from local measurements; however, a comparison of Bayesian information criteria reveal that this model and the CIB & SZ only model are equally plausible. Thus, more sensitive measurements will be needed to confirm the existence of large-scale CII emission at high redshifts. Finally, we forecast that intensity maps from Planck cross-correlated with quasars from the Dark Energy Spectroscopic Instrument (DESI) would increase our sensitivity to CII emission by a factor of 5, while the proposed Primordial Inflation Explorer (PIXIE) could increase the sensitivity further.

We assess the effect of magnification bias on the $E_G$ statistic for probing gravity. $E_G$, a statistic constructed from power spectrum estimates of both weak lensing and redshift space distortions (RSD), directly tests general relativity (GR) while in principle being independent of clustering bias. This property has motivated its recent use in multiple tests of GR. Recent work has suggested that the magnification bias of galaxies due to foreground matter perturbations breaks the bias-independence of the $E_G$ statistic. The magnitude of this effect is very sensitive to the clustering and magnification biases of the galaxy sample. We show that for realistic values of the clustering and magnification biases, the effect for magnification bias is small relative to statistical errors for most spectroscopic galaxy surveys but large for photometric galaxy surveys. For the cases with significant magnification bias, we propose a method to calibrate magnification bias using measurements of the lensing auto-power spectrum. We test this calibration method using simulations, finding that our calibration method can calibrate $E_G$ from 2-4 times the simulation error to well within the errors, although its accuracy is sensitive to the precision of the measured redshift distribution and magnification bias, but not the clustering bias. This work gives strong evidence that this method will work increasingly well in future CMB lensing surveys.

Following the first two annual intensity mapping workshops at Stanford in March 2016 and Johns Hopkins in June 2017, we report on the recent advances in theory, instrumentation and observation that were presented in these meetings and some of the opportunities and challenges that were identified looking forward. With preliminary detections of CO, [CII], Lya and low-redshift 21cm, and a host of experiments set to go online in the next few years, the field is rapidly progressing on all fronts, with great anticipation for a flood of new exciting results. This current snapshot provides an efficient reference for experts in related fields and a useful resource for nonspecialists. We begin by introducing the concept of line-intensity mapping and then discuss the broad array of science goals that will be enabled, ranging from the history of star formation, reionization and galaxy evolution to measuring baryon acoustic oscillations at high redshift and constraining theories of dark matter, modified gravity and dark energy. After reviewing the first detections reported to date, we survey the experimental landscape, presenting the parameters and capabilities of relevant instruments such as COMAP, mmIMe, AIM-CO, CCAT-p, TIME, CONCERTO, CHIME, HIRAX, HERA, STARFIRE, MeerKAT/SKA and SPHEREx. Finally, we describe recent theoretical advances: different approaches to modeling line luminosity functions, several techniques to separate the desired signal from foregrounds, statistical methods to analyze the data, and frameworks to generate realistic intensity map simulations.

We demonstrate a new method to constrain gravity on the largest cosmological scales by combining measurements of cosmic microwave background (CMB) lensing and the galaxy velocity field. $E_G$ is a statistic, constructed from a gravitational lensing tracer and a measure of velocities such as redshift-space distortions (RSD), that can discriminate between gravity models while being independent of clustering bias and $\sigma_8$. While traditionally, the lensing field for $E_G$ has been probed through galaxy lensing, CMB lensing has been proposed as a more robust tracer of the lensing field for $E_G$ at higher redshifts while avoiding intrinsic alignments. We perform the largest-scale measurement of $E_G$ ever, up to 150 Mpc/$h$, by cross-correlating the Planck CMB lensing map with the Sloan Digital Sky Survey III (SDSS-III) CMASS galaxy sample and combining this with our measurement of the CMASS auto-power spectrum and the RSD parameter $\beta$. We report $E_G(z=0.57)=0.243\pm0.060$ (stat) $\pm0.013$ (sys), a measurement in tension with the general relativity prediction at a level of 2.6$\sigma$. Note that our $E_G$ measurement deviates from GR only at scales greater than 80 Mpc/$h$, scales which have not been probed by previous $E_G$ tests. Upcoming surveys, which will provide an order-of-magnitude reduction in statistical errors, can significantly constrain alternative gravity models when combined with better control of systematics.

We present Intercut, a Python-based program that applies secondary line identification and photometric cuts to mock galaxy surveys, in order to simulate interloper identification. This program can be used to optimize the removal of interloper contamination in upcoming surveys. Intercut reads a mock galaxy survey and an emission line sensitivity and simulates interloper removal through secondary line identification and broad-band photometry. This program is designed to use the COSMOS mock catalog, although the program can be modified for an alternative mock catalog. The output of the program returns an interloper fraction for each emission line, as well as the total fraction over all lines, as a function of redshift. We test Intercut by predicting interloper rates for the WFIRST emission line sensitivity, finding agreement with previous work. This program is publically available on Github

Next-generation spectroscopic surveys will map the large-scale structure of the observable universe, using emission line galaxies as tracers. While each survey will map the sky with a specific emission line, interloping emission lines can masquerade as the survey's intended emission line at different redshifts. Interloping lines from galaxies that are not removed can contaminate the power spectrum measurement, mixing correlations from various redshifts and diluting the true signal. We assess the potential for power spectrum contamination, finding that an interloper fraction worse than 0.2% could bias power spectrum measurements for future surveys by more than 10% of statistical errors, while also biasing power spectrum inferences. We also construct a formalism for predicting cosmological parameter bias, demonstrating that a 0.15%-0.3% interloper fraction could bias the growth rate by more than 10% of the error, which can affect constraints on gravity upcoming surveys. We use the COSMOS Mock Catalog (CMC), with the emission lines re-scaled to better reproduce recent data, to predict potential interloper fractions for the Prime Focus Spectrograph (PFS) and the Wide-Field InfraRed Survey Telescope (WFIRST). We find that secondary line identification, or confirming galaxy redshifts by finding correlated emission lines, can remove interlopers for PFS. For WFIRST, we use the CMC to predict that the 0.2% target can be reached for the WFIRST H$\alpha$ survey, but sensitive optical and near-infrared photometry will be required. For the WFIRST [OIII] survey, the predicted interloper fractions reach several percent and their effects will have to be estimated and removed statistically (e.g. with deep training samples). (Abridged)

We describe a methodology to probe gravity with the cosmic microwave background (CMB) lensing convergence $\kappa$, specifically by measuring $E_G$, the ratio of the Laplacian of the gravitational scalar potential difference with the velocity divergence. Using CMB lensing instead of galaxy-galaxy lensing avoids intrinsic alignments while also lacking a hard limit on the lens redshift and significant uncertainties in the source plane. We model $E_G$ for general relativity and modified gravity, finding that $E_G$ for $f(R)$ gravity should be scale-dependent due to the scale-dependence of the growth rate $f$. Next, we construct an estimator for $E_G$ in terms of the galaxy-CMB lensing and galaxy clustering angular power spectra, along with the RSD parameter $\beta$. We also forecast statistical errors for $E_G$ from the current Planck CMB lensing map and the spectroscopic galaxy and quasar samples from the Sloan Digital Sky Survey Data Release 11, being 9% with galaxies and 8% when quasars are included. We also find that upcoming spectroscopic and photometric surveys, combined with the final Planck lensing map, can measure precisely the redshift- and scale-dependence of $E_G$ out to redshifts $z=2$ and higher, with photometric surveys having an advantage due to their high number densities. Advanced ACTPol's lensing map will increase the $E_G$ sensitivity even further. Finally, we find that Advanced ACTPol cross-correlated with spectroscopic (photometric) surveys can differentiate between general relativity and $f(R)$ gravity at the level of $3\sigma$ ($13\sigma$). Performing a $<1$% measurement of $E_G$ requires a 10% precision in $\beta$ from Euclid or LSST, currently achievable with a spectroscopic survey but difficult with only a photometric survey.

SPHEREx (Spectro-Photometer for the History of the Universe, Epoch of Reionization, and Ices Explorer) ( http://spherex.caltech.edu ) is a proposed all-sky spectroscopic survey satellite designed to address all three science goals in NASA's Astrophysics Division: probe the origin and destiny of our Universe; explore whether planets around other stars could harbor life; and explore the origin and evolution of galaxies. SPHEREx will scan a series of Linear Variable Filters systematically across the entire sky. The SPHEREx data set will contain R=40 spectra fir 0.75$<\lambda<$4.1$\mu$m and R=150 spectra for 4.1$<\lambda<$4.8$\mu$m for every 6.2 arc second pixel over the entire-sky. In this paper, we detail the extra-galactic and cosmological studies SPHEREx will enable and present detailed systematic effect evaluations. We also outline the Ice and Galaxy Evolution Investigations.

We measure the cross-power spectra between luminous red galaxies (LRGs) from the Sloan Digital Sky Survey (SDSS)-III Data Release Eight (DR8) and cosmic infrared background (CIB) anisotropies from Planck and data from the Improved Reprocessing (IRIS) of the Infrared Astronomical Satellite (IRAS) at 353, 545, 857, and 3000 GHz, corresponding to 850, 550, 350 and 100 micron, respectively, in the multipole range 100<l<1000. Using approximately 6.5 10^5 photometrically determined LRGs in 7760 deg^2 of the northern hemisphere in the redshift range 0.45 < z < 0.65, we model the far-infrared background (FIRB) anisotropies with an extended version of the halo model. With these methods, we confirm the basic picture obtained from recent analyses of FIRB anisotropies with Herschel and Planck, that the most efficient halo mass at hosting star forming galaxies is log(M_ eff/M_\odot)=12.84+/-0.15. We estimate the percentage of FIRB anisotropies correlated with LRGs as approximately 11.8 %, 3.9 %, 1.8 %, and 1.0 % of the total at 3000, 857, 545, and 353 GHz, respectively. At redshift z~0.55, the bias of FIRB galaxies with respect to the dark matter density field has the value b_{FIRB}~1.45, and the mean dust temperature of FIRB galaxies is T_d=26 K. Finally, we discuss the impact of present and upcoming cross-correlations with far-infrared background anisotropies on the determination of the global star formation history and the link between galaxies and dark matter.

We describe the methodology to include nonlinear evolution, including tidal effects, in the computation of subhalo distribution properties in both cold (CDM) and warm (WDM) dark matter universes. Using semi-analytic modeling, we include effects from dynamical friction, tidal stripping, and tidal heating, allowing us to dynamically evolve the subhalo distribution. We calibrate our nonlinear evolution scheme to the CDM subhalo mass function in the Aquarius N-body simulation, producing a subhalo mass function within the range of simulations. We find tidal effects to be the dominant mechanism of nonlinear evolution in the subhalo population. Finally, we compute the subhalo mass function for $m_\chi=1.5$ keV WDM including the effects of nonlinear evolution, and compare radial number densities and mass density profiles of subhalos in CDM and WDM models. We show that all three signatures differ between the two dark matter models, suggesting that probes of substructure may be able to differentiate between them.

On the arcminute angular scales probed by Planck, the CMB anisotropies are gently perturbed by gravitational lensing. Here we present a detailed study of this effect, detecting lensing independently in the 100, 143, and 217GHz frequency bands with an overall significance of greater than 25sigma. We use the temperature-gradient correlations induced by lensing to reconstruct a (noisy) map of the CMB lensing potential, which provides an integrated measure of the mass distribution back to the CMB last-scattering surface. Our lensing potential map is significantly correlated with other tracers of mass, a fact which we demonstrate using several representative tracers of large-scale structure. We estimate the power spectrum of the lensing potential, finding generally good agreement with expectations from the best-fitting LCDM model for the Planck temperature power spectrum, showing that this measurement at z=1100 correctly predicts the properties of the lower-redshift, later-time structures which source the lensing potential. When combined with the temperature power spectrum, our measurement provides degeneracy-breaking power for parameter constraints; it improves CMB-alone constraints on curvature by a factor of two and also partly breaks the degeneracy between the amplitude of the primordial perturbation power spectrum and the optical depth to reionization, allowing a measurement of the optical depth to reionization which is independent of large-scale polarization data. Discarding scale information, our measurement corresponds to a 4% constraint on the amplitude of the lensing potential power spectrum, or a 2% constraint on the RMS amplitude of matter fluctuations at z~2.

The ESA's Planck satellite, dedicated to studying the early Universe and its subsequent evolution, was launched 14 May 2009 and has been scanning the microwave and submillimetre sky continuously since 12 August 2009. This paper gives an overview of the mission and its performance, the processing, analysis, and characteristics of the data, the scientific results, and the science data products and papers in the release. The science products include maps of the CMB and diffuse extragalactic foregrounds, a catalogue of compact Galactic and extragalactic sources, and a list of sources detected through the SZ effect. The likelihood code used to assess cosmological models against the Planck data and a lensing likelihood are described. Scientific results include robust support for the standard six-parameter LCDM model of cosmology and improved measurements of its parameters, including a highly significant deviation from scale invariance of the primordial power spectrum. The Planck values for these parameters and others derived from them are significantly different from those previously determined. Several large-scale anomalies in the temperature distribution of the CMB, first detected by WMAP, are confirmed with higher confidence. Planck sets new limits on the number and mass of neutrinos, and has measured gravitational lensing of CMB anisotropies at greater than 25 sigma. Planck finds no evidence for non-Gaussianity in the CMB. Planck's results agree well with results from the measurements of baryon acoustic oscillations. Planck finds a lower Hubble constant than found in some more local measures. Some tension is also present between the amplitude of matter fluctuations derived from CMB data and that derived from SZ data. The Planck and WMAP power spectra are offset from each other by an average level of about 2% around the first acoustic peak.

We present a quantitative model of Lyman-alpha (Ly-alpha) emission throughout cosmic history and determine the prospects for intensity mapping of spatial fluctuations in the Ly-alpha signal. Since i) our model assumes at z>6 the minimum star formation required to sustain reionization and ii) is based at z<6 on a luminosity function extrapolated from the few observed bright Ly-alpha emitters, this should be considered a lower limit. Mapping the line emission allows probes of reionization, star formation, and large-scale structure (LSS) as a function of redshift. While Ly-alpha emission during reionization has been studied, we also predict the post-reionization signal to test predictions of the intensity and motivate future intensity mapping probes of reionization. We include emission from massive dark matter halos and the intergalactic medium (IGM) in our model. We find agreement with current, measured luminosity functions of Ly-alpha emitters at z<8. However, diffuse IGM emission, not associated with Ly-alpha emitters, dominates the intensity up to z~10. While our model is applicable for deep-optical or near-infrared observers like the SuMIre survey or the James Webb Space Telescope, only intensity mapping will detect the diffuse IGM emission. We also construct a 3D power spectrum model of the Ly\alpha emission, and characterize possible foreground contamination. Finally, we study the prospects of an intensity mapper for measuring Ly-alpha fluctuations while identifying interloper contamination for removal. Our results suggest that while the reionization signal is challenging, Ly-alpha fluctuations can be an interesting new probe of LSS at late times when used in conjunction with other lines like, e.g., H-alpha, to monitor low-redshift foreground confusion.

Primordial non-Gaussianity of local type is predicted to lead to enhanced halo clustering on very large scales. Photometric quasars, which can be seen from cosmological redshifts z>2 even in wide-shallow optical surveys, are promising tracers for constraining non-Gaussianity using this effect. However, large-scale systematics can also mimic this signature of non-Gaussianity. In order to assess the contribution of systematic effects, we cross-correlate overdensity maps of photometric quasars from the Sloan Digital Sky Survey (SDSS) Data Release 6 (DR6) in different redshift ranges. We find that the maps are significantly correlated on large scales, even though we expect the angular distributions of quasars at different redshifts to be uncorrelated. This implies that the quasar maps are contaminated with systematic errors. We investigate the use of external templates that provide information on the spatial dependence of potential systematic errors to reduce the level of spurious clustering in the quasar data. We find that templates associated with stellar density, the stellar color locus, airmass, and seeing are major contaminants of the quasar maps, with seeing having the largest effect. Using template projection, we are able to decrease the significance of the cross-correlation measurement on the largest scales from 9.2-sigma to 5.4-sigma. Although this is an improvement, the remaining cross-correlation suggests the contamination in this quasar sample is too great to allow a competitive constraint on fNL by correlations internal to this sample. The SDSS quasar catalog exhibits spurious number density fluctuations of ~2% rms, and we need a contamination level less than 1% (0.6%) in order to measure values of fNL less than 100 (10). Properly dealing with these systematics will be paramount for future large scale structure surveys that seek to constrain non-Gaussianity.

We present a new procedure to measure the large-scale carbon monoxide (CO) emissions across cosmic history. As a tracer of large-scale structure (LSS), the CO gas content as a function of redshift can be quantified by its three-dimensional fluctuation power spectra. Furthermore, cross-correlating CO emission with other LSS tracers offers a way to measure the emission as a function of scale and redshift. Here we introduce the model relevant for such a cross-correlation measurement between CO and other LSS tracers, and between different CO rotational lines. We propose a novel use of Cosmic Microwave Background (CMB) data and attempt to extract redshifted CO emissions embedded in the Wilkinson Microwave Anisotropy Probe (WMAP) dataset. We cross-correlate the all-sky WMAP7 data with LSS data sets, namely, the photometric quasar sample and the luminous red galaxy sample from the Sloan Digital Sky Survey Data Release 6 and 7 respectively. We are unable to detect a cross-correlation signal with either CO(1-0) nor CO(2-1) lines, mainly due to the instrumental noise in the WMAP data. However, we are able to rule out models more than three times greater than our more optimistic model. We discuss the cross-correlation signal from the thermal Sunyaev-Zeldovich effect and dust as potential contaminants, and quantify their impact for our CO measurements. We discuss forecasts for current CMB experiments and a hypothetical future CO focused experiment, and propose to cross-correlate CO temperature data with the Hobby-Eberly Telescope Dark Energy Experiment Ly-alpha-emitter sample, for which a signal-to-noise ratio of 58 is possible.

We search a sample of photometric luminous red galaxies (LRGs) measured by the Sloan Digital Sky Survey (SDSS) for a quadrupolar anisotropy in the primordial power spectrum, in which P(\vec{k}) is an isotropic power spectrum P(k) multiplied by a quadrupolar modulation pattern. We first place limits on the 5 coefficients of a general quadrupole anisotropy. We also consider axisymmetric quadrupoles of the form P(\vec{k}) = P(k){1 + g_*[(\hat{k}\cdot\hat{n})^2-1/3]} where \hat{n} is the axis of the anisotropy. When we force the symmetry axis \hat{n} to be in the direction (l,b)=(94 degrees,26 degrees) identified in the recent Groeneboom et al. analysis of the cosmic microwave background, we find g_*=0.006+/-0.036 (1 sigma). With uniform priors on \hat{n} and g_* we find that -0.41<g_*<+0.38 with 95% probability, with the wide range due mainly to the large uncertainty of asymmetries aligned with the Galactic Plane. In none of these three analyses do we detect evidence for quadrupolar power anisotropy in large scale structure.

Statistical isotropy of primordial perturbations is a common assumption in cosmology, but it is an assumption that should be tested. To this end, we develop cosmic microwave background statistics for a primordial power spectrum that depends on the direction, as well as the magnitude, of the Fourier wavevector. We first consider a simple estimator that searches in a model-independent way for anisotropy in the square of the temperature (and/or polarization) fluctuation. We then construct the minimum-variance estimators for the coefficients of a spherical-harmonic expansion of the direction-dependence of the primordial power spectrum. To illustrate, we apply these statistics to an inflation model with a quadrupole dependence of the primordial power spectrum on direction and find that a power quadrupole as small as 2.0% can be detected with the Planck satellite.

We search data from the Energetic Gamma Ray Experiment Telescope (EGRET) for a gamma-ray line in the energy range 0.1-10 GeV from the 10 X 10 degree region around the Galactic center. Our null results lead to upper limits to the line flux from the Galactic center. Such lines may have appeared if the dark matter in the Galactic halo is composed of weakly-interacting massive particles (WIMPs) in the mass range 0.1-10 GeV. For a given dark-matter-halo model, our null search translates to upper limits to the WIMP two-photon annihilation cross section as a function of WIMP mass. We show that for a toy model in which Majorana WIMPs in this mass range annihilate only to electron-positron pairs, these upper limits supersede those derived from measurements of the 511-keV line and continuum photons from internal bremsstrahlung at the Galactic center.