The classic Bondi solution remains a common starting point both for studying black hole growth across cosmic time in cosmological simulations and for smaller scale simulations of AGN feedback. In nature, however, there will be inhomogenous distributions of rotational velocity and density along the outer radius ($R_o$) marking the sphere of influence of a black hole. While there have been many studies of how the Bondi solution changes with a prescribed angular momentum boundary condition, they have all assumed a constant density at $R_o$. In this Letter, we show that a non-uniform density at $R_o$ causes a meridional flow and due to conservation of angular momentum, the Bondi solution qualitatively changes into an inflow-outflow solution. Using physical arguments, we analytically identify the critical logarithmic density gradient $|\partial{\ln{\rho}}/\partial{\theta}|$ above which this change of the solution occurs. For realistic $R_o$, this critical gradient is less than 0.01 and tends to 0 as $R_o \rightarrow \infty$. We show using numerical simulations that, unlike for solutions with an imposed rotational velocity, the accretion rate for solutions under an inhomogenous density boundary condition remains constant at nearly the Bondi rate $\dot{M}_B$, while the outflow rate can greatly exceed $\dot{M}_B$.

Supermassive black holes (SMBHs) are the central engines of luminous quasars and are found in most massive galaxies today. But the recent discoveries of ULAS J1120+0641, a $2 \times 10^9$ M$_{\odot}$ BH at $z =$ 7.1, and ULAS J1342+0928, a $8.0 \times 10^{8}$ M$_{\odot}$ BH at $z =$ 7.5, now push the era of quasar formation up to just 690 Myr after the Big Bang. Here we report new cosmological simulations of SMBHs with X-rays fully coupled to primordial chemistry and hydrodynamics that show that J1120 and J1342 can form from direct collapse black holes (DCBHs) if their growth is fed by cold, dense accretion streams, like those thought to fuel rapid star formation in some galaxies at later epochs. Our models reproduce all of the observed properties of J1120: its mass, luminosity, and H II region as well as star formation rates and metallicities in its host galaxy. They also reproduce the dynamical mass of the innermost 1.5 kpc of its emission region recently measured by ALMA and J-band magnitudes that are in good agreement with those found by the VISTA Hemisphere Survey.

The simultaneous detection of gravitational and electromagnetic waves from a binary neutron star merger has both solidified the link between neutron star mergers and short-duration gamma-ray bursts (GRBs) and demonstrated the ability of astronomers to follow-up the gravitational wave detection to place constraints on the ejecta from these mergers as well as the nature of the GRB engine and its surroundings. As the sensitivity of aLIGO and VIRGO increases, it is likely that a growing number of such detections will occur in the next few years, leading to a sufficiently-large number of events to constrain the populations of these GRB events. While long-duration GRBs originate from massive stars and thus are located near their stellar nurseries, binary neutron stars may merge on much longer timescales, and thus may have had time to migrate appreciably. The strength and character of the electromagnetic afterglow emission of binary neutron star mergers is a sensitive function of the circum-merger environment. Though the explosion sites of short GRBs have been explored in the literature, the question has yet to be fully addressed in its cosmological context. We present cosmological simulations following the evolution of a galaxy cluster including star formation combined with binary population synthesis models to self-consistently track the locations and environmental gas densities of compact binary merger sites throughout the cosmic web. We present probability distributions for densities as a function of redshift and discuss model sensitivity to population synthesis model assumptions.

We present a photometric catalog for Spitzer Space Telescope warm mission observations of the North Ecliptic Pole (NEP; centered at $\rm R.A.=18^h00^m00^s$, $\rm Decl.=66^d33^m38^s.552$). The observations are conducted with IRAC in 3.6 $\mu$m and 4.5 $\mu$m bands over an area of 7.04 deg$^2$ reaching 1$\sigma$ depths of 1.29 $\mu$Jy and 0.79 $\mu$Jy in the 3.6 $\mu$m and 4.5 $\mu$m bands respectively. The photometric catalog contains 380,858 sources with 3.6 $\mu$m and 4.5 $\mu$m band photometry over the full-depth NEP mosaic. Point source completeness simulations show that the catalog is 80% complete down to 19.7 AB. The accompanying catalog can be utilized in constraining the physical properties of extra-galactic objects, studying the AGN population, measuring the infrared colors of stellar objects, and studying the extra-galactic infrared background light.

Supernovae are important probes of the properties of stars at high redshifts because they can be detected at early epochs and their masses can be inferred from their light curves. Finding the first cosmic explosions in the universe will only be possible with the James Webb Space Telescope, the Wide-Field Infrared Survey Telescope and the next generation of extremely large telescopes. But strong gravitational lensing by massive clusters, like those in the Cluster Lensing and Supernova Survey with Hubble (CLASH), could reveal such events now by magnifying their flux by factors of 10 or more. We find that CLASH will likely discover at least 2 - 3 core-collapse supernovae at 5 < z < 12 and perhaps as many as ten. Future surveys of cluster lenses similar in scope to CLASH by the James Webb Space Telescope might find hundreds of these events out to z ~ 15 - 17. Besides revealing the masses of early stars, these ancient supernovae will also constrain cosmic star formation rates in the era of first galaxy formation.

We present near-infrared (0.8-1.8 $\mu$m) spectra of 105 bright (${m_{J}}$ $<$ 10) stars observed with the low resolution spectrometer on the rocket-borne Cosmic Infrared Background Experiment (CIBER). As our observations are performed above the earth's atmosphere, our spectra are free from telluric contamination, which makes them a unique resource for near-infrared spectral calibration. Two-Micron All Sky Survey (2MASS) photometry information is used to identify cross-matched stars after reduction and extraction of the spectra. We identify the spectral types of the observed stars by comparing them with spectral templates from the Infrared Telescope Facility (IRTF) library. All the observed spectra are consistent with late F to M stellar spectral types, and we identify various infrared absorption lines.

We present the first ab initio cosmological simulations of a CR7-like object which approximately reproduce the observed line widths and strengths. In our model, CR7 is powered by a massive ($3.23 \times 10^7$ $M_\odot$) black hole (BH) the accretion rate of which varies between $\simeq$ 0.25 and $\simeq$ 0.9 times the Eddington rate on timescales as short as 10$^3$ yr. Our model takes into account multi-dimensional effects, X-ray feedback, secondary ionizations and primordial chemistry. We estimate Ly-$\alpha$ line widths by post-processing simulation output with Monte Carlo radiative transfer and calculate emissivity contributions from radiative recombination and collisional excitation. We find the luminosities in the Lyman-$\alpha$ and He II 1640 angstrom lines to be $5.0\times10^{44}$ and $2.4\times10^{43}$ erg s$^{-1}$, respectively, in agreement with the observed values of $>$ $8.3\times10^{43}$ and $2.0\times10^{43}$ erg s$^{-1}$. We also find that the black hole heats the halo and renders it unable to produce stars as required to keep the halo metal free. These results demonstrate the viability of the BH hypothesis for CR7 in a cosmological context. Assuming the BH mass and accretion rate that we find, we estimate the synchrotron luminosity of CR7 to be $P \simeq 10^{40} - 10^{41}$ erg s$^{-1}$, which is sufficiently luminous to be observed in $\mu$Jy observations and would discriminate this scenario from one where the luminosity is driven by Population III stars.

The extragalactic background light (EBL) is comprised of the cumulative radiation from all galaxies and active galactic nuclei over the cosmic history. In addition to point sources, EBL also contains information from diffuse sources of radiation. The angular power spectra of the near-infrared intensities could contain additional signals and a complete understanding of the nature of the IR background is still lacking in the literature. Here we explore the constraints that can be placed on particle decays, especially candidate dark matter models involving axions that trace dark matter halos of galaxies. Axions with a mass around a few eV will decay via two photons with wavelengths in the near-IR band, and will leave a signature in the IR background intensity power spectrum. Using recent power spectra measurements from the Hubble Space Telescope (HST) and Cosmic Infrared Background Experiment (CIBER), we find that the 0.6 to 1.6 micron power spectra can be explained by axions with masses around 4 eV. The total axion abundance Omega_a~0.05, and it is comparable to the baryon density of the Universe. The suggested mean axion mass and abundance are not ruled out by existing cosmological observations. Interestingly, the axion model with a mass distribution is preferred by the data, which cannot be explained by the standard quantum chromodynamics (QCD) theory and needs further discussion.

Questions surround the connection of luminous extragalactic masers to galactic processes. The observation that water and hydroxyl megamasers rarely coexist in the same galaxy has given rise to a hypothesis that the two species appear in different phases of nuclear activity. The detection of simultaneous hydroxyl and water megamaser emission toward IC694 has called this hypothesis into question but, because many megamasers have not been surveyed for emission in the other molecule, it remains unclear whether IC694 occupies a narrow phase of galaxy evolution or whether the relationship between megamaser species and galactic processes is more complicated than previously believed. In this paper, we present results of a systematic search for 22 GHz water maser emission among OH megamaser hosts to identify additional objects hosting both megamaser. Our work roughly doubles the number of galaxies searched for emission in both molecules which host at least one confirmed maser. We confirm with high degree of confidence ($> 8 \sigma$) the detection of water emission toward IIZw96, firmly establishing it as the second object to co-host both water and hydroxyl megamasers after IC694. We find high luminosity, narrow features in the water feature in IIZw96. All dual megamaser candidates appear in merging galaxy systems suggestive that megamaser coexistance may signal a brief phase along the merger sequence. A statistical analysis of the results of our observations provide possible evidence for an exclusion of H$_2$O kilomasers among OH megamaser hosts.

Pair Instability Supernovae have been suggested as candidates for some Super Luminous Supernovae, such as SN 2007bi, and as one of the dominant types of explosion occurring in the early Universe from massive, zero-metallicity Population III stars. The progenitors of such events can be rapidly rotating, therefore exhibiting different evolutionary properties due to the effects of rotationally-induced mixing and mass-loss. Proper identification of such events requires rigorous radiation hydrodynamics and radiative transfer calculations that capture not only the behavior of the light curve but also the spectral evolution of these events. We present radiation hydrodynamics and radiation transport calculations for 90-300 Msun rotating pair-instability supernovae covering both the shock break-out and late light curve phases. We also investigate cases of different initial metallicity and rotation rate to determine the impact of these parameters on the detailed spectral characteristics of these events. In agreement with recent results on non-rotating pair instability supernovae, we find that for a range of progenitor masses and rotation rates these events have intrinsically red colors in contradiction with observations of super-luminous supernovae. The spectroscopic properties of rotating pair instability supernovae are similar to those of non-rotating events with stripped hydrogen and helium envelopes. We find that the progenitor metallicity and rotation rate properties are erased after the explosion and cannot be identified in the resulting model spectra. It is the combined effects of pre-supernova mass-loss and the basic properties of the supernova ejecta such as mass, temperature and velocity that have the most direct impact in the model spectra of pair instability supernovae.

Population III supernovae have been of growing interest of late for their potential to directly probe the properties of the first stars, particularly the most energetic events that are visible near the edge of the observable universe. But until now, hypernovae, the unusually energetic Type Ib/c supernovae that are sometimes associated with gamma-ray bursts, have been overlooked as cosmic beacons at the highest redshifts. In this, the latest of a series of studies on Population III supernovae, we present numerical simulations of 25 - 50 M$_{\odot}$ hypernovae and their light curves done with the Los Alamos RAGE and SPECTRUM codes. We find that they will be visible at z = 10 - 15 to the James Webb Space Telescope (JWST) and z = 4 - 5 to the Wide-Field Infrared Survey Telescope (WFIRST), tracing star formation rates in the first galaxies and at the end of cosmological reionization. If, however, the hypernova crashes into a dense shell ejected by its progenitor, it is expected that a superluminous event will occur that may be seen at z ~ 20, in the first generation of stars.

Gamma-ray bursts (GRBs) are the ultimate cosmic lighthouses, capable of illuminating the universe at its earliest epochs. Could such events probe the properties of the first stars at z $\sim$ 20, the end of the cosmic Dark Ages? Previous studies of Population III GRBs only considered explosions in the diffuse relic H II regions of their progenitors, or bursts that are far more more energetic than those observed to date. But the processes that produce GRBs at the highest redshifts likely reset their local environments, creating much more complicated structures than those in which relativistic jets have been modeled so far. These structures can greatly affect the luminosity of the afterglow, and hence the redshift at which it can be detected. We have now simulated Population III GRB afterglows in H II regions, winds, and dense shells ejected by the star during the processes that produce the burst. Our model, which has been used in previous work, has been extended to include contributions from reverse shocks, inverse Compton cooling and the effects of sphericity and beaming in the blast wave, and is valid in a variety of circumjet density profiles. We find that GRBs with E$_{\mathrm{iso},\gamma} =$ 10$^{ 51}$ - 10$^{53}$ erg will be visible at z $\gtrsim$ 20 to the next generation of near infrared and radio observatories. In many cases, the environment of the burst, and hence progenitor type, can be inferred from the afterglow light curve. Although some Population III GRBs are visible to Swift and the Very Large Array now, the optimal strategy for their detection will be future missions like EXIST and JANUS, which have large survey areas and onboard X-ray and infrared telescopes that can track their near infrared flux from the moment of the burst, thereby identifying its redshift.

The spatial fluctuations of the extragalactic background light trace the total emission from all stars and galaxies in the Universe. A multi-wavelength study can be used to measure the integrated emission from first galaxies during reionization when the Universe was about 500 million years old. Here we report arcminute-scale spatial fluctuations in one of the deepest sky surveys with the Hubble Space Telescope in five wavebands between 0.6 and 1.6 $\mu$m. We model-fit the angular power spectra of intensity fluctuation measurements to find the ultraviolet luminosity density of galaxies at $z$ > 8 to be $\log \rho_{\rm UV} = 27.4^{+0.2}_{-1.2}$ erg s$^{-1}$ Hz$^{-1}$ Mpc$^{-3}$ $(1\sigma)$. This level of integrated light emission allows for a significant surface density of fainter primeval galaxies that are below the point source detection level in current surveys.

The era of the universe's first (Population III) stars is essentially unconstrained by observation. Ultra-luminous and massive stars from this time altered the chemistry of the cosmos, provided the radiative scaffolding to support the formation of the first protogalaxies, and facilitated the creation and growth of now-supermassive black holes. Unfortunately, because these stars lie literally at the edge of the observable universe, they will remain beyond the reach of even the next generation of telescopes such as the James Webb Space Telescope and the Thirty-Meter Telescope. In this paper, we provide a primer to supernovae modeling and the first stars to make our discussion accessible to those new to or outside our field. We review recent work of the Los Alamos Supernova Light Curve Project and Brigham Young University to explore the possibility of probing this era through observations of the spectacular deaths of the first stars. We find that many such brilliant supernova explosions will be observable as far back as $\sim 99$% of the universe's current age, tracing primordial star formation rates and the locations of their protogalaxies on the sky. The observation of Population III supernovae will be among the most spectacular discoveries in observational astronomy in the coming decade.

We report measurements of the Diffuse Galactic Light (DGL) spectrum in the near-infrared, spanning the wavelength range 0.95-1.65 {\mu}m by the Cosmic Infrared Background ExpeRiment (CIBER). Using the low-resolution spectrometer (LRS) calibrated for absolute spectro-photometry, we acquired long-slit spectral images of the total diffuse sky brightness towards four high-latitude fields spread over four sounding rocket flights. To separate the DGL spectrum from the total sky brightness, we correlated the spectral images with a 100 {\mu}m intensity map, which traces the dust column density in optically thin regions. The measured DGL spectrum shows no resolved features and is consistent with other DGL measurements in the optical and at near-infrared wavelengths longer than 1.8 {\mu}m. Our result implies that the continuum is consistently reproduced by models of scattered starlight in the Rayleigh scattering regime with a few large grains.

Non-Gaussianity of the primordial density perturbations provides an important measure to constrain models of inflation. At cubic order the non-Gaussianity is captured by two parameters $\tau_{\rm NL}$ and $g_{\rm NL}$ that determine the amplitude of the density perturbation trispectrum. Here we report measurements of the kurtosis power spectra of the cosmic microwave background (CMB) temperature as mapped by Planck by making use of correlations between square temperature-square temperature and cubic temperature-temperature anisotropies. In combination with noise simulations, we find the best joint estimates to be $\tau_{\rm{NL}}=0.3 \pm 0.9 \times 10^4$ and $g_{\rm{NL}}=-1.2 \pm 2.8 \times 10^5$. If $\tau_{\rm NL}=0$, we find $g_{\rm NL}= -1.3\pm 1.8 \times 10^5$.

We present the cross-correlation between the far-infrared background fluctuations as measured with the Herschel Space Observatory at 250, 350, and 500 {\mu}m and the near-infrared background fluctuations with Spitzer Space Telescope at 3.6 {\mu}m. The cross-correlation between far and near-IR background anisotropies are detected such that the correlation coefficient at a few to ten arcminute angular scales decreases from 0.3 to 0.1 when the far-IR wavelength increases from 250 {\mu}m to 500 {\mu}m. We model the cross-correlation using a halo model with three components: (a) far-IR bright or dusty star-forming galaxies below the masking depth in Herschel maps, (b) near-IR faint galaxies below the masking depth at 3.6 {\mu}m, and (c) intra-halo light, or diffuse stars in dark matter halos, that likely dominates fluctuations at 3.6 {\mu}m. The model is able to reasonably reproduce the auto correlations at each of the far-IR wavelengths and at 3.6 {\mu}m and their corresponding cross-correlations. While the far and near-IR auto-correlations are dominated by faint dusty, star-forming galaxies and intra-halo light, respectively, we find that roughly half of the cross-correlation between near and far-IR backgrounds is due to the same galaxies that remain unmasked at 3.6 {\mu}m. The remaining signal in the cross-correlation is due to intra-halo light present in the same dark matter halos as those hosting the same faint and unmasked galaxies. In this model, the decrease in the cross-correlation signal from 250 {\mu}m to 500 {\mu}m comes from the fact that the galaxies that are primarily contributing to 500 {\mu}m fluctuations peak at a higher redshift than those at 250 {\mu}m.

The discovery of 150 - 300 M$_{\odot}$ stars in the Local Group and pair-instability supernova candidates at low redshifts has excited interest in this exotic explosion mechanism. Realistic light curves for pair-instability supernovae at near-solar metallicities are key to identifying and properly interpreting these events as more are found. We have modeled pair-instability supernovae of 150 - 500 M$_{\odot}$ Z $\sim$ 0.1 - 0.4 Z$_{\odot}$ stars. These stars lose up to 80% of their mass to strong line-driven winds and explode as bare He cores. We find that their light curves and spectra are quite different from those of Population III pair-instability explosions, which therefore cannot be used as templates for low-redshift events. Although non-zero metallicity pair-instability supernovae are generally dimmer than their Population III counterparts, in some cases they will be bright enough to be detected at the earliest epochs at which they can occur, the formation of the first galaxies at $z \sim$ 10 - 15. Others can masquerade as dim, short duration supernovae that are only visible in the local universe and that under the right conditions could be hidden in a wide variety of supernova classes. We also report for the first time that some pair-instability explosions can create black holes with masses of $\sim$ 100 M$_{\odot}$.

Population III stars that die as pair-instability supernovae are usually thought to fall in the mass range of 140 - 260 M$_{\odot}$. But several lines of work have now shown that rotation can build up the He cores needed to encounter the pair instability at stellar masses as low as 90 $_{\odot}$. Depending on the slope of the initial mass function of Population III stars, there could be 4 - 5 times as many stars from 90 - 140 $_{\odot}$ in the primordial universe than in the usually accepted range. We present numerical simulations of the pair-instability explosions of such stars performed with the MESA, FLASH and RAGE codes. We find that they will be visible to supernova factories such as Pan-STARRS and LSST in the optical out to z $\sim$ 1 - 2 and to JWST and the 30 m-class telescopes in the NIR out to $z \sim$ 7 - 10. Such explosions will thus probe the stellar populations of the first galaxies and cosmic star formation rates in the era of cosmological reionization. These supernovae are also easily distinguished from more massive pair-instability explosions, underscoring the fact that there is far greater variety to the light curves of these events than previously understood.

Extragalactic background light (EBL) anisotropy traces variations in the total production of photons over cosmic history, and may contain faint, extended components missed in galaxy point source surveys. Infrared EBL fluctuations have been attributed to primordial galaxies and black holes at the epoch of reionization (EOR), or alternately, intra-halo light (IHL) from stars tidally stripped from their parent galaxies at low redshift. We report new EBL anisotropy measurements from a specialized sounding rocket experiment at 1.1 and 1.6 micrometers. The observed fluctuations exceed the amplitude from known galaxy populations, are inconsistent with EOR galaxies and black holes, and are largely explained by IHL emission. The measured fluctuations are associated with an EBL intensity that is comparable to the background from known galaxies measured through number counts, and therefore a substantial contribution to the energy contained in photons in the cosmos.

We present novel statistical tools to cross-correlate frequency cleaned thermal Sunyaev-Zel'dovich (tSZ) maps and tomographic weak lensing (wl) convergence maps. Moving beyond the lowest order cross-correlation, we introduce a hierarchy of mixed higher-order statistics, the cumulants and cumulant correlators, to analyze non-Gaussianity in real space, as well as corresponding polyspectra in the harmonic domain. Using these moments, we derive analytical expressions for the joint two-point probability distribution function (2PDF) for smoothed tSZ (y_s) and convergence (\kappa_s) maps. The presence of tomographic information allows us to study the evolution of higher order {\em mixed} tSZ-weak lensing statistics with redshift. We express the joint PDFs p_{\kappa y}(\kappa_s,y_s) in terms of individual one-point PDFs (p_{\kappa}(\kappa_s), p_y(y_s)) and the relevant bias functions (b_{\kappa}(\kappa_s), b_y(y_s)). Analytical results for two different regimes are presented that correspond to the small and large angular smoothing scales. Results are also obtained for corresponding {\em hot spots} in the tSZ and convergence maps. In addition to results based on hierarchical techniques and perturbative methods, we present results of calculations based on the lognormal approximation. The analytical expressions derived here are generic and applicable to cross-correlation studies of arbitrary tracers of large scale structure including e.g. that of tSZ and soft X-ray background.

Population III supernovae have been the focus of growing attention because of their potential to directly probe the properties of the first stars, particularly the most energetic events that can be seen at the edge of the observable universe. But until now pair-pulsation supernovae, in which explosive thermonuclear burning in massive stars fails to unbind them but can eject their outer layers into space, have been overlooked as cosmic beacons at the earliest redshifts. These shells can later collide and, like Type IIn supernovae, produce superluminous events in the UV at high redshifts that could be detected in the near infrared today. We present numerical simulations of a 110 M$_{\odot}$ pair-pulsation explosion done with the Los Alamos radiation hydrodynamics code RAGE. We find that collisions between consecutive pair pulsations are visible in the near infrared out to z $\sim$ 15 - 20 and can probe the earliest stellar populations at cosmic dawn.

The existence of supermassive black holes as early as z ~ 7 is one of the great unsolved problems in cosmological structure formation. One leading theory argues that they are born during catastrophic baryon collapse in z ~ 15 protogalaxies in strong Lyman-Werner UV backgrounds. Atomic line cooling in such galaxies fragments baryons into massive clumps that are thought to directly collapse to 10^4 - 10^5 solar-mass black holes. We have now discovered that some of these fragments can instead become supermassive stars that eventually explode as pair-instability supernovae with energies of ~ 10^55 erg, the most energetic explosions in the universe. We have calculated light curves and spectra for supermassive Pop III PI SNe with the Los Alamos RAGE and SPECTRUM codes. We find that they will be visible in NIR all-sky surveys by Euclid out to z ~ 10 - 15 and by WFIRST and WISH out to z ~ 15 - 20, perhaps revealing the birthplaces of the first quasars.

One of the leading contenders for the origin of supermassive black holes at $z \gtrsim$ 7 is catastrophic baryon collapse in atomically-cooled halos at $z \sim$ 15. In this scenario, a few protogalaxies form in the presence of strong Lyman-Werner UV backgrounds that quench H$_2$ formation in their constituent halos, preventing them from forming stars or blowing heavy elements into the intergalactic medium prior to formation. At masses of 10$^ 8$ \Ms\ and virial temperatures of 10$^4$ K, gas in these halos rapidly cools by H lines, in some cases forming 10$^4$ - 10$^6$ \Ms\ Pop III stars and, a short time later, the seeds of supermassive black holes. Instead of collapsing directly to black holes some of these stars died in the most energetic thermonuclear explosions in the universe. We have modeled the explosions of such stars in the dense cores of line-cooled protogalaxies in the presence of cosmological flows. In stark contrast to the explosions in diffuse regions in previous simulations, these SNe briefly engulf the protogalaxy but then collapse back into its dark matter potential. Fallback drives turbulence that efficiently distributes metals throughout the interior of the halo and fuels the rapid growth of nascent black holes at its center. The accompanying starburst and x-ray emission from these line-cooled galaxies easily distinguish them from more slowly evolving neighbors and might reveal the birthplaces of supermassive black holes on the sky.

The first stars are the key to the formation of primitive galaxies, early cosmological reionization and chemical enrichment, and the origin of supermassive black holes. Unfortunately, in spite of their extreme luminosities, individual Population III stars will likely remain beyond the reach of direct observation for decades to come. However, their properties could be revealed by their supernova explosions, which may soon be detected by a new generation of NIR observatories such as JWST and WFIRST. We present light curves and spectra for Pop III pair-instability supernovae calculated with the Los Alamos radiation hydrodynamics code RAGE. Our numerical simulations account for the interaction of the blast with realistic circumstellar envelopes, the opacity of the envelope, and Lyman absorption by the neutral IGM at high redshift, all of which are crucial to computing the NIR signatures of the first cosmic explosions. We find that JWST will detect pair-instability supernovae out to z > 30, WFIRST will detect them in all-sky surveys out to z ~ 15 - 20 and LSST and Pan-STARRS will find them at z ~ 7 - 8. The discovery of these ancient explosions will probe the first stellar populations and reveal the existence of primitive galaxies that might not otherwise have been detected.

Supermassive primordial stars are expected to form in a small fraction of massive protogalaxies in the early universe, and are generally conceived of as the progenitors of the seeds of supermassive black holes (BHs). Supermassive stars with masses of ~55,000 M_Sun, however, have been found to explode and completely disrupt in a supernova (SN) with an energy of up to ~10^55 erg instead of collapsing to a BH. Such events, ~10,000 times more energetic than typical SNe today, would be among the biggest explosions in the history of the universe. Here we present a simulation of such a SN in two stages. Using the RAGE radiation hydrodynamics code we first evolve the explosion from an early stage through the breakout of the shock from the surface of the star until the blast wave has propagated out to several parsecs from the explosion site, which lies deep within an atomic cooling dark matter (DM) halo at z ~ 15. Then, using the GADGET cosmological hydrodynamics code we evolve the explosion out to several kiloparsecs from the explosion site, far into the low-density intergalactic medium. The host DM halo, with a total mass of 4 x 10^7 M_Sun, much more massive than typical primordial star-forming halos, is completely evacuated of high density gas after < 10 Myr, although dense metal-enriched gas recollapses into the halo, where it will likely form second-generation stars with metallicities of ~ 0.05 Z_Sun after > 70 Myr. The chemical signature of supermassive star explosions may be found in such long-lived second-generation stars today.

The first primitive galaxies formed from accretion and mergers by z ~ 15, and were primarily responsible for cosmological reionization and the chemical enrichment of the early cosmos. But a few of these galaxies may have formed in the presence of strong Lyman-Werner UV fluxes that sterilized them of H_2, preventing them from forming stars or expelling heavy elements into the IGM prior to assembly. At masses of 10^8 Ms and virial temperatures of 10^4 K, these halos began to rapidly cool by atomic lines, perhaps forming 10^4 - 10^6 Ms Pop III stars and, later, the seeds of supermassive black holes. We have modeled the explosion of a supermassive Pop III star in the dense core of a line-cooled protogalaxy with the ZEUS-MP code. We find that the supernova (SN) expands to a radius of ~ 1 kpc, briefly engulfing the entire galaxy, but then collapses back into the potential well of the dark matter. Fallback fully mixes the interior of the protogalaxy with metals, igniting a violent starburst and fueling the rapid growth of a massive black hole at its center. The starburst would populate the protogalaxy with stars in greater numbers and at higher metallicities than in more slowly-evolving, nearby halos. The SN remnant becomes a strong synchrotron source that can be observed with eVLA and eMERLIN and has a unique signature that easily distinguishes it from less energetic SN remnants. Such explosions, and their attendant starbursts, may well have marked the birthplaces of supermassive black holes on the sky.

We generalize the concept of the ordinary skew-spectrum to probe the effect of non-Gaussianity on the morphology of Cosmic Microwave Background (CMB) maps in several domains: in real-space (where they are commonly known as cumulant-correlators), and in harmonic and needlet bases. The essential aim is to retain more information than normally contained in these statistics, in order to assist in determining the source of any measured non-Gaussianity, in the same spirit as Munshi & Heavens (2010) skew-spectra were used to identify foreground contaminants to the CMB bispectrum in Planck data. Using a perturbative series to construct the Minkowski Functionals (MFs), we provide a pseudo-Cl based approach in both harmonic and needlet representations to estimate these spectra in the presence of a mask and inhomogeneous noise. Assuming homogeneous noise we present approx- imate expressions for error covariance for the purpose of joint estimation of these spectra. We present specific results for four different models of primordial non-Gaussianity local, equilateral, orthogonal and enfolded models, as well as non-Gaussianity caused by unsubtracted point sources. Closed form results of next-order corrections to MFs too are obtained in terms of a quadruplet of kurt-spectra. We also use the method of modal decomposition of the bispectrum and trispectrum to reconstruct the MFs as an alternative method of reconstruction of morphological properties of CMB maps. Finally, we introduce the odd-parity skew-spectra to probe the odd-parity bispectrum and its impact on the morphology of the CMB sky. Although developed for the CMB, the generic results obtained here can be useful in other areas of cosmology.

We use the Cosmic Microwave Background (CMB) anisotropy data from Planck to constrain the spatial fluctuations of the fine-structure constant \alpha. Through Thompson scattering of CMB photons, spatial anisotropies of \alpha lead to higher-order correlations in the CMB anisotropies. We use a quadratic estimator based on the four-point correlation function of the CMB temperature anisotropy to extract the angular power spectrum of the spatial variation of the fine-structure constant projected along the line of sight at the last scattering surface. At tens of degree angular scales and above, we constrain the rms fluctuations of the fine structure constant to be \delta \alpha/\alpha_0= (1.34 +/- 5.82) x 10^-2 at the 95% confidence level with respect to the standard value \alpha_0. We find no evidence for a spatially varying \alpha at a redshift of 10^3.

We present a measurement of the angular power spectrum of the cosmic far-infrared background (CFIRB) anisotropies in one of the extragalactic fields of the Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS) at 250, 350 and 500 \mu m bands. Consistent with recent measurements of the CFIRB power spectrum in Herschel-SPIRE maps, we confirm the existence of a clear one-halo term of galaxy clustering on arcminute angular scales with large-scale two-halo term of clustering at 30 arcminutes to angular scales of a few degrees. The power spectrum at the largest angular scales, especially at 250 \mu m, is contaminated by the Galactic cirrus. The angular power spectrum is modeled using a conditional luminosity function approach to describe the spatial distribution of unresolved galaxies that make up the bulk of the CFIRB. Integrating over the dusty galaxy population responsible for the background anisotropies, we find that the cosmic abundance of dust, relative to the critical density, to be between \Omega_dust=10^{-6} and 8 x 10^{-6} in the redshift range z ~ 0-3. This dust abundance is consistent with estimates of the dust content in the Universe using quasar reddening and magnification measurements in the SDSS.

At high angular frequencies, beyond the damping tail of the primary power spectrum, the dominant contribution to the power spectrum of cosmic microwave background (CMB) temperature fluctuations is the thermal Sunyaev-Zel'dovich (tSZ) effect. We investigate various important statistical properties of the Sunyaev-Zel'dovich maps, using well-motivated models for dark matter clustering to construct statistical descriptions of the tSZ effect to all orders enabling us to determine the entire probability distribution function (PDF). Any generic deterministic biasing scheme can be incorporated in our analysis and the effects of projection, biasing and the underlying density distribution can be analysed separately and transparently in this approach. We introduce the cumulant correlators as tools to analyse tSZ catalogs and relate them to corresponding statistical descriptors of the underlying density distribution. The statistics of hot spots in frequency-cleaned tSZ maps are also developed in a self-consistent way to an arbitrary order, to obtain results complementary to those found using the halo model. We also consider different beam sizes, to check the extent to which the PDF can be extracted from various observational configurations. The formalism is presented with two specific models for underlying matter clustering: (1) the hierarchical ansatz; and (2) the lognormal distribution. We find both models to be in very good agreement with the simulation results, though the lognormal model has an edge over the hierarchical model.

Unresolved near-infrared background anisotropies are expected to have contributions from the earliest galaxies during reionization and faint, dwarf galaxies at intermediate redshifts. Previous measurements were unable to conclusively pinpoint the dominant origin because they did not sample spatial scales that were sufficiently large to distinguish between these two possibilities. Here we report a measurement of the anisotropy power spectrum from sub-arcminute to one degree angular scales and find the clustering amplitude to be larger than the model predictions involving the two existing explanations. As the shot-noise level of the power spectrum is consistent with that expected from faint galaxies, a new source population on the sky is not necessary to explain the observations. A physical mechanism that increases the clustering amplitude, however, is needed. Motivated by recent results related to the extended stellar light profile in dark matter halos, we consider the possibility that the fluctuations originate from diffuse intrahalo stars of all galaxies. We find that the measured power spectrum can be explained by an intrahalo light fraction of 0.07 to 0.2 % relative to the total luminosity in dark matter halos of masses log(M/M_Sun) ~ 9 to 12 at redshifts of ~ 1 to 4.

A fraction of the extragalactic near-infrared (near-IR) background light involves redshifted photons from the ultraviolet (UV) emission from galaxies present during reionization at redshifts above 6. The absolute intensity and the anisotropies of the near-IR background provide an observational probe of the first-light galaxies and their spatial distribution. We estimate the extragalactic background light intensity during reionization by accounting for the stellar and nebular emission from first-light galaxies. We require the UV photon density from these galaxies to generate a reionization history that is consistent with the optical depth to electron scattering from cosmic microwave background measurements. We also require the bright-end luminosity function of galaxies in our models to reproduce the measured Lyman drop-out luminosity functions at redshifts of 6 to 8. The absolute intensity is about 0.1 to 0.3 nW m$^{-2}$ sr$^{-1}$ at the peak of its spectrum at $\sim$ 1.1 $\mu$m. We also discuss the anisotropy power spectrum of the near-IR background using a halo model to describe the galaxy distribution. We compare our predictions for the anisotropy power spectrum to existing measurements from deep near-IR imaging data from {\it Spitzer}/IRAC, {\it Hubble}/NICMOS, and {\it AKARI}. The predicted rms fluctuations at tens of arcminute angular scales are roughly an order of magnitude smaller than the existing measurements. While strong arguments have been made that the measured fluctuations do not have an origin involving faint low-redshift galaxies, we find that the existing measurements are also incompatible with an origin during the era of reionization. The measured near-IR background anisotropies remain unexplained and could be associated with an unidentified non-astrophysical origin.

We have developed and characterized an imaging instrument to measure the spatial properties of the diffuse near-infrared extragalactic background light in a search for fluctuations from z > 6 galaxies during the epoch of reionization. The instrument is part of the Cosmic Infrared Background Experiment (CIBER), designed to observe the extragalactic background light above the Earth's atmosphere during a suborbital sounding rocket flight. The imaging instrument incorporates a 2x2 degree field of view, to measure fluctuations over the predicted peak of the spatial power spectrum at 10 arcminutes, and 7"x7" pixels, to remove lower redshift galaxies to a depth sufficient to reduce the low-redshift galaxy clustering foreground below instrumental sensitivity. The imaging instrument employs two cameras with \Delta \lambda / \lambda ~0.5 bandpasses centered at 1.1 and 1.6 microns to spectrally discriminate reionization extragalactic background fluctuations from local foreground fluctuations. CIBER operates at wavelengths where the electromagnetic spectrum of the reionization extragalactic background is thought to peak, and complements fluctuations measurements by AKARI and Spitzer at longer wavelengths. We have characterized the instrument in the laboratory, including measurements of the sensitivity, flat-field response, stray light performance, and noise properties. Several modifications were made to the instrument following a first flight in 2009 February. The instrument performed to specifications in subsequent flights in 2010 July and 2012 March, and the scientific data are now being analyzed.

Cosmological data have provided new constraints on the number of neutrino species and the neutrino mass. However these constraints depend on assumptions related to the underlying cosmology. Since a correlation is expected between the number of effective neutrinos N_{eff}, the neutrino mass \sum m_\nu, and the curvature of the universe \Omega_k, it is useful to investigate the current constraints in the framework of a non-flat universe. In this paper we update the constraints on neutrino parameters by making use of the latest cosmic microwave background (CMB) data from the ACT and SPT experiments and consider the possibility of a universe with non-zero curvature. We first place new constraints on N_{eff} and \Omega_k, with N_{eff} = 4.03 +/- 0.45 and 10^3 \Omega_k = -4.46 +/- 5.24. Thus, even when \Omega_k is allowed to vary, N_{eff} = 3 is still disfavored with 95% confidence. We then investigate the correlation between neutrino mass and curvature that shifts the 95% upper limit of \sum m_\nu < 0.45 eV to \sum m_\nu < 0.95 eV. Thus, the impact of assuming flatness in neutrino cosmology is significant and an essential consideration with future experiments.

The wide-area imaging surveys with the {\it Herschel} Space Observatory at sub-mm wavelengths have now resulted in catalogs of order one hundred thousand dusty, star-burst galaxies. We make a statistical estimate of $N(z)$ using a clustering analysis of sub-mm galaxies detected at each of 250, 350 and 500 $\mu$m from the Herschel Multi-tiered Extragalactic Survey (HerMES) centered on the Bo\"{o}tes field. We cross-correlate {\it Herschel} galaxies against galaxy samples at optical and near-IR wavelengths from the Sloan Digital Sky Survey (SDSS), the NOAO Deep Wide Field Survey (NDWFS) and the Spitzer Deep Wide Field Survey (SDWFS). We create optical and near-IR galaxy samples based on their photometric or spectroscopic redshift distributions and test the accuracy of those redshift distributions with similar galaxy samples defined with catalogs of the Cosmological Evolution Survey (COSMOS), as the COSMOS field has superior spectroscopy coverage. We model-fit the clustering auto and cross-correlations of {\it Herschel} and optical/IR galaxy samples to estimate $N(z)$ and clustering bias factors. The $S_{350} > 20$ mJy galaxies have a bias factor varying with redshift as $b(z)=1.0^{+1.0}_{-0.5}(1+z)^{1.2^{+0.3}_{-0.7}}$. This bias and the redshift dependence is broadly in agreement with galaxies that occupy dark matter halos of mass in the range of 10$^{12}$ to 10$^{13}$ M$_{\sun}$. We find that the redshift distribution peaks around $z \sim 0.5$ to 1 for galaxies selected at 250 $\mu$m with an average redshift of $< z > = 1.8 \pm 0.2$. For 350 and 500 $\mu$m-selected SPIRE samples the peak shifts to higher redshift, but the average redshift remains the same with a value of $1.9 \pm 0.2$.

This paper describes the Hubble Space Telescope imaging data products and data reduction procedures for the Cosmic Assembly Near-IR Deep Extragalactic Legacy Survey (CANDELS). This survey is designed to document the evolution of galaxies and black holes at $z\sim1.5-8$, and to study Type Ia SNe beyond $z>1.5$. Five premier multi-wavelength sky regions are selected, each with extensive multiwavelength observations. The primary CANDELS data consist of imaging obtained in the Wide Field Camera 3 / infrared channel (WFC3/IR) and UVIS channel, along with the Advanced Camera for Surveys (ACS). The CANDELS/Deep survey covers \sim125 square arcminutes within GOODS-N and GOODS-S, while the remainder consists of the CANDELS/Wide survey, achieving a total of \sim800 square arcminutes across GOODS and three additional fields (EGS, COSMOS, and UDS). We summarize the observational aspects of the survey as motivated by the scientific goals and present a detailed description of the data reduction procedures and products from the survey. Our data reduction methods utilize the most up to date calibration files and image combination procedures. We have paid special attention to correcting a range of instrumental effects, including CTE degradation for ACS, removal of electronic bias-striping present in ACS data after SM4, and persistence effects and other artifacts in WFC3/IR. For each field, we release mosaics for individual epochs and eventual mosaics containing data from all epochs combined, to facilitate photometric variability studies and the deepest possible photometry. A more detailed overview of the science goals and observational design of the survey are presented in a companion paper.

The Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS) is designed to document the first third of galactic evolution, over the approximate redshift (z) range 8--1.5. It will image >250,000 distant galaxies using three separate cameras on the Hubble Space Telescope, from the mid-ultraviolet to the near-infrared, and will find and measure Type Ia supernovae at z>1.5 to test their accuracy as standardizable candles for cosmology. Five premier multi-wavelength sky regions are selected, each with extensive ancillary data. The use of five widely separated fields mitigates cosmic variance and yields statistically robust and complete samples of galaxies down to a stellar mass of 10^9 M_\odot to z \approx 2, reaching the knee of the ultraviolet luminosity function (UVLF) of galaxies to z \approx 8. The survey covers approximately 800 arcmin^2 and is divided into two parts. The CANDELS/Deep survey (5\sigma\ point-source limit H=27.7 mag) covers \sim 125 arcmin^2 within GOODS-N and GOODS-S. The CANDELS/Wide survey includes GOODS and three additional fields (EGS, COSMOS, and UDS) and covers the full area to a 5\sigma\ point-source limit of H \gtrsim 27.0 mag. Together with the Hubble Ultra Deep Fields, the strategy creates a three-tiered "wedding cake" approach that has proven efficient for extragalactic surveys. Data from the survey are nonproprietary and are useful for a wide variety of science investigations. In this paper, we describe the basic motivations for the survey, the CANDELS team science goals and the resulting observational requirements, the field selection and geometry, and the observing design. The Hubble data processing and products are described in a companion paper.

We present results from a search for gravitationally lensed galaxies present in the Hubble Space Telescope (HST)/Wide Field Camera-3 (WFC3) images of the Cosmic Assembly Near-IR Deep Extragalactic Legacy Survey (CANDELS). We present one bona fide lens system in UDS and two compact lens candidates in the GOODS-S field. The lensing system in UDS involves two background galaxies, one at z=1.847 lensed to an arc and a counterimage, and the second at a photometric redshift of z=2.32^{+0.10}_{-0.06} lensed to a double image. We reconstruct the lensed sources in the source plane and find in each of the two cases the sources can be separated to a pair of galaxies. The sources responsible for the arc are compact with effective radii of 0.3 to 0.4 kpc in WFC3 J_{125}-band and a total stellar mass and a star-formation rate of 2.1_{-0.4}^{+2.4} times 10^7 M_sun and 2.3_{-1.7}^{+ 0.6} M_sun yr^{-1}, respectively.The abnormally high H_{160}-band flux of this source is likely due to OIII emission lines with a rest-frame equivalent width about 700 Angstroms for OIII 5007 Angstroms. The sources responsible for the double image have corresponding values of about 0.4 to 0.5 kpc, 1.4_{-0.8}^{+1.9} times 10^9 M_sun, and 8.7_{-7.0}^{+11.1} M_sun yr^{-1}. Once completed CANDELS is expected to contain about 15 lensing systems and will allow statistical studies on both lensing mass profiles and z ~ 2 lensed galaxies.

The Sunyaev-Zel'dovich (SZ) effect has a distinct spectral signature that allows its separation from fluctuations in the cosmic microwave background (CMB) and foregrounds. Using CMB anisotropies measured in Wilkinson Microwave Anisotropy Probe's five-year maps, we constrain the SZ fluctuations at large, degree angular scales corresponding to multipoles in the range from 10 to 400. We provide upper bounds on SZ fluctuations at multipoles greater than 50, and find evidence for a hemispherically asymmetric signal at ten degrees angular scales. The amplitude of the detected signal cannot be easily explained with the allowed number density and temperature of electrons in the Galactic halo. We have failed to explain the excess signal as a residual from known Galactic foregrounds or instrumental uncertainties such as 1/f-noise.

At high angular frequencies the thermal Sunyaev-Zel'dovich (tSZ) effect constitutes the dominant signal in the CMB sky. The tSZ effect is caused by large scale pressure fluctuations in the baryonic distribution in the Universe so its statistical properties provide estimates of corresponding properties of the projected 3D pressure fluctuations. It's power spectrum is a sensitive probe of the density fluctuations, and the bispectrum can be used to separate the bias associated with pressure. The bispectrum is often probed with a one-point real-space analogue, the skewness. In addition to the skewness the morphological properties, as probed by the well known Minkowski Functionals (MFs), also require the generalized one-point statistics, which at the lowest order are identical to the skewness parameters. The concept of generalized skewness parameters can be extended to define a set of three associated generalized skew-spectra. We use these skew-spectra to probe the morphology of the tSZ sky or the y-sky. We show how these power spectra can be recovered from the data in the presence of arbitrary mask and noise templates using the well known Pseudo-Cl (PCL) approach for arbitrary beam shape. We also employ an approach based on the halo model to compute the tSZ bispectrum. The bispectrum from each of these models is then used to construct the generalized skew-spectra. We consider the performance of an all-sky survey with Planck-type noise and compare the results against a noise-free ideal experiment using a range of smoothing angles. We find that the skew-spectra can be estimated with very high signal-to-noise ratio from future frequency cleaned tSZ maps that will be available from experiments such as Planck. This will allow their mode by mode estimation for a wide range of angular frequencies and will help us to differentiate them from various other sources of non-Gaussianity.

We present a new harmonic-domain approach for extracting morphological information, in the form of Minkowski Functionals (MFs), from weak lensing (WL) convergence maps. Using a perturbative expansion of the MFs, which is expected to be valid for the range of angular scales probed by most current weak-lensing surveys, we show that the study of three generalized skewness parameters is equivalent to the study of the three MFs defined in two dimensions. We then extend these skewness parameters to three associated skew-spectra which carry more information about the convergence bispectrum than their one-point counterparts. We discuss various issues such as noise and incomplete sky coverage in the context of estimation of these skew-spectra from realistic data. Our technique provides an alternative to the pixel-space approaches typically used in the estimation of MFs, and it can be particularly useful in the presence of masks with non-trivial topology. Analytical modeling of weak lensing statistics relies on an accurate modeling of the statistics of underlying density distribution. We apply three different formalisms to model the underlying dark-matter bispectrum: the hierarchical ansatz, halo model and a fitting function based on numerical simulations; MFs resulting from each of these formalisms are computed and compared. We investigate the extent to witch late-time gravity-induced non-Gaussianity (to which weak lensing is primarily sensitive) can be separated from primordial non-Gaussianity and how this separation depends on source redshift and angular scale.

Ultraviolet emission from the first generation of stars in the Universe ionized the intergalactic medium in a process which was completed by z~6; the wavelength of these photons has been redshifted by (1+z) into the near infrared today and can be measured using instruments situated above the Earth's atmosphere. First flying in February 2009, the Cosmic Infrared Background Experiment (CIBER) comprises four instruments housed in a single reusable sounding rocket borne payload. CIBER will measure spatial anisotropies in the extragalactic IR background caused by cosmological structure from the epoch of reionization using two broadband imaging instruments, make a detailed characterization of the spectral shape of the IR background using a low resolution spectrometer, and measure the absolute brightness of the Zodical light foreground with a high resolution spectrometer in each of our six science fields. This paper presents the scientific motivation for CIBER and details of its first two flights, including a review of the published scientific results from the first flight and an outlook for future reionization science with CIBER data.

The temperature fluctuations and polarization of the Cosmic Microwave Background (CMB) are now a well-known probe of the Universe at an infant age of 400,000 years. During the transit to us from the surface of last scattering, the CMB photons are expected to undergo modifications induced by the intervening large-scale structure. Among the expected secondary effects is the weak gravitational lensing of the CMB by the foreground dark matter distribution. We derive a quadratic estimator that uses the non-Gaussianities generated by the lensing effect at the four-point function level to extract the power spectrum of lensing potential fluctuations integrated out to z ~ 1100 with peak contributions from potential fluctuations at z of 2 to 3. Using WMAP 7-year temperature maps, we report the first direct constraints of this lensing potential power spectrum and find that it has an amplitude of A_L = 0.96 +/- 0.60, 1.06 +/- 0.69 and 0.97 +/- 0.47 using the W, V and W+V bands, respectively.

We outline the expected constraints on non-Gaussianity from the cosmic microwave background (CMB) with current and future experiments, focusing on both the third (f_{NL}) and fourth-order (g_{NL} and \tau_{NL}) amplitudes of the local configuration or non-Gaussianity. The experimental focus is the skewness (two-to-one) and kurtosis (two-to-two and three-to-one) power spectra from weighted maps. In adition to a measurement of \tau_{NL} and g_{NL} with WMAP 5-year data, our study provides the first forecasts for future constraints on g_{NL}. We describe how these statistics can be corrected for the mask and cut-sky through a window function, bypassing the need to compute linear terms that were introduced for the previous-generation non-Gaussianity statistics, such as the skewness estimator. We discus the ratio A_{NL} = \tau_{NL}/(6f_{NL}/5)^2 as an additional test of single-field inflationary models and discuss the physical significance of each statistic. Using these estimators with WMAP 5-Year V+W-band data out to l_{max}=600 we constrain the cubic order non-Gaussianity parameters \tau_{NL}, and g_{NL} and find -7.4 < g_{NL}/10^5 < 8.2 and -0.6 < \tau_{NL}/10^4 < 3.3 improving the previous COBE-based limit on \tau_{NL} < 10^8 nearly four orders of magnitude with WMAP.

A large sub-mm survey with Herschel will enable many exciting science opportunities, especially in an era of wide-field optical and radio surveys and high resolution cosmic microwave background experiments. The Herschel-SPIRE Legacy Survey (HSLS), will lead to imaging data over 4000 sq. degrees at 250, 350, and 500 micron. Major Goals of HSLS are: (a) produce a catalog of 2.5 to 3 million galaxies down to 26, 27 and 33 mJy (50% completeness; 5 sigma confusion noise) at 250, 350 and 500 micron, respectively, in the southern hemisphere (3000 sq. degrees) and in an equatorial strip (1000 sq. degrees), areas which have extensive multi-wavelength coverage and are easily accessible from ALMA. Two thirds of the of the sources are expected to be at z > 1, one third at z > 2 and about a 1000 at z > 5. (b) Remove point source confusion in secondary anisotropy studies with Planck and ground-based CMB data. (c) Find at least 1200 strongly lensed bright sub-mm sources leading to a 2% test of general relativity. (d) Identify 200 proto-cluster regions at z of 2 and perform an unbiased study of the environmental dependence of star formation. (e) Perform an unbiased survey for star formation and dust at high Galactic latitude and make a census of debris disks and dust around AGB stars and white dwarfs.

We consider corrections to the underlying cosmology due to secondary contributions from weak gravitational lensing, the integrated Sachs-Wolfe effect, and the Sunyaev-Zel'dovich effect contained in the trispectrum. We incorporate these additional contributions to the covariance of a binned angular power spectrum of temperature anisotropies in the analysis of current and prospective data sets. Although recent experiments such as ACBAR and CBI are not particularly sensitive to these additional non-Gaussian effects, the interpretation of Planck and CMBPol anisotropy spectra will require an accounting of non-Gaussian covariance leading to a degradation in cosmological parameter estimates by up to 20% and 30%, respectively.

We measure two higher-order power spectra involving weighted cubic and squared temperature anisotropy maps from WMAP 5-year data to study the trispectrum generated by primordial non-Gaussianity. Using these measurements combined with Gaussian and noise simulations, we constrain the cubic order non-Gaussianity parameters \tau_{NL}, and g_{NL}. With V+W-band data out to l_{max}=600, we find -7.4 < g_{\rm NL}/10^5 < 8.2 and -0.6 < \tau_{\rm NL}/10^4 < 3.3 improving the previous COBE-based limit on \tau_{\rm NL} < 10^8 nearly four orders of magnitude with WMAP. We find that the ratio of trispectrum to bispectrum amplitude as captured by the ratio of \tau_{\rm Nl}/(6f_{\rm NL}/5)^2 ranges from -3 to 21 at the 95% confidence level.

Owing to their more extensive sky coverage and tighter control on systematic errors, future deep weak lensing surveys should provide a better statistical picture of the dark matter clustering beyond the level of the power spectrum. In this context, the study of non-Gaussianity induced by gravity can help tighten constraints on the background cosmology by breaking parameter degeneracies, as well as throwing light on the nature of dark matter, dark energy or alternative gravity theories. Analysis of the shear or flexion properties of such maps is more complicated than the simpler case of the convergence due to the spinorial nature of the fields involved. Here we develop analytical tools for the study of higher-order statistics such as the bispectrum (or trispectrum) directly using such maps at different source redshift. The statistics we introduce can be constructed from cumulants of the shear or flexions, involving the cross-correlation of squared and cubic maps at different redshifts. Typically, the low signal-to-noise ratio prevents recovery of the bispectrum or trispectrum mode by mode. We define power spectra associated with each multi- spectra which compresses some of the available information of higher order multispectra. We show how these can be recovered from a noisy observational data even in the presence of arbitrary mask, which introduces mixing between Electric (E-type) and Magnetic (B-type) polarization, in an unbiased way. We also introduce higher order cross-correlators which can cross-correlate lensing shear with different tracers of large scale structures.

We explore a systematic approach to the analysis of primordial non-Gaussianity using fluctuations in temperature and polarization of the Cosmic Microwave Background (CMB). Following Munshi & Heavens (2009), we define a set of power-spectra as compressed forms of the bispectrum and trispectrum derived from CMB temperature and polarization maps; these spectra compress the information content of the corresponding full multispectra and can be useful in constraining early Universe theories. We generalize the standard pseudo-C_l estimators in such a way that they apply to these spectra involving both spin-0 and spin-2 fields, developing explicit expressions which can be used in the practical implementation of these estimators. While these estimators are suboptimal, they are nevertheless unbiased and robust hence can provide useful diagnostic tests at a relatively small computational cost. We next consider approximate inverse-covariance weighting of the data and construct a set of near-optimal estimators based on that approach. Instead of combining all available information from the entire set of mixed bi- or trispectra, i.e multispectra describing both temperature and polarization information, we provide analytical constructions for individual estimators, associated with particular multispectra. The bias and scatter of these estimators can be computed using Monte-Carlo techniques. Finally, we provide estimators which are completely optimal for arbitrary scan strategies and involve inverse covariance weighting; we present the results of an error analysis performed using a Fisher-matrix formalism at both the one-point and two-point level.

Cosmic microwave background studies of non-Gaussianity involving higher-order multispectra can distinguish between early universe theories that predict nearly identical power spectra. However, the recovery of higher-order multispectra is difficult from realistic data due to their complex response to inhomogeneous noise and partial sky coverage, which are often difficult to model analytically. A traditional alternative is to use one-point cumulants of various orders, which collapse the information present in a multispectrum to one number. The disadvantage of such a radical compression of the data is a loss of information as to the source of the statistical behaviour. A recent study by Munshi & Heavens (2009) has shown how to define the skew spectrum (the power spectra of a certain cubic field, related to the bispectrum) in an optimal way and how to estimate it from realistic data. The skew spectrum retains some of the information from the full configuration-dependence of the bispectrum, and can contain all the information on non-Gaussianity. In the present study, we extend the results of the skew spectrum to the case of two degenerate power-spectra related to the trispectrum. We also explore the relationship of these power-spectra and cumulant correlators previously used to study non-Gaussianity in projected galaxy surveys or weak lensing surveys. We construct nearly optimal estimators for quick tests and generalise them to estimators which can handle realistic data with all their complexity in a completely optimal manner. We show how these higher-order statistics and the related power spectra are related to the Taylor expansion coefficients of the potential in inflation models, and demonstrate how the trispectrum can constrain both the quadratic and cubic terms.

We measure the skewness power spectrum of the Cosmic Microwave Background (CMB) anisotropies optimized for a detection of the secondary bispectrum generated by the correlation of the CMB lensing potential with integrated Sachs-Wolfe effect and the Sunyaev-Zel'dovich effect. The covariance of our measurements is generated by Monte-Carlo simulations of Gaussian CMB fields with noise properties consistent with Wilkinson Microwave Anisotropy Probe (WMAP) 5-year data. When interpreting multi-frequency measurements we also take into account the confusion resulting from unresolved radio point sources. We analyze Q, V and W-band WMAP 5-year raw and foreground-cleaned maps using the KQ75 mask out to l_max=600. We find no significant evidence for a non-zero non-Gaussian signal from the lensing-secondary correlation in all three bands and we constrain the overall amplitude of the cross power spectrum between CMB lensing potential and the sum of SZ and ISW fluctuations to be 0.42 \pm 0.86 and 1.19 \pm 0.86 in combined V and W-band raw and foreground-cleaned maps provided by the WMAP team, respectively. The point source amplitude at the bispectrum level measured with this skewness power spectrum is higher than previous measurements of point source non-Gaussianity. We also consider an analysis where we also account for the primordial non-Gaussianity in addition to lensing-secondary bispectrum and point sources. The focus of this paper is on secondary anisotropies. Consequently the estimator is not optimised for primordial non-Gaussianity and the limit we find on local non-Gaussianity from the foreground-cleaned V+W maps is f_NL = -13 \pm 62, when marginalized over point sources and lensing-ISW/SZ constributions to the total bispectrum.

We constrain the primordial non-Gaussianity parameter of the local model f_{NL} using the skewness power spectrum associated with the two-to-one cumulant correlator of cosmic microwave background temperature anisotropies. This bispectrum-related power spectrum was constructed after weighting the temperature map with the appropriate window functions to form an estimator that probes the multipolar dependence of the underlying bispectrum associated with the primordial non-Gaussianity. We also estimate a separate skewness power spectrum sensitive more strongly to unresolved point sources. When compared to previous attempts at measuring the primordial non-Gaussianity with WMAP data, our estimators have the main advantage that we do not collapse information to a single number. When model fitting the two-to-one skewness power spectrum we make use of bispectra generated by the primordial non-Gaussianity, radio point sources, and lensing-secondary correlation. We analyze Q, V and W-band WMAP 5-year data using the KQ75 mask out to l_{max}=600. Using V and W-band data and marginalizing over model parameters related to point sources and lensing-secondary bispectrum, our overall and preferred constraint on f_{NL} is 11.0\pm23.7 at the 68% confidence level (-36.4 < f_{NL} < 58.4 at 95% confidence). We find no evidence for a non-zero value of f_{NL} even marginally at the 1\sigma level.

The development of radiation hydrodynamical methods that are able to follow gas dynamics and radiative transfer self-consistently is key to the solution of many problems in numerical astrophysics. Such fluid flows are highly complex, rarely allowing even for approximate analytical solutions against which numerical codes can be tested. An alternative validation procedure is to compare different methods against each other on common problems, in order to assess the robustness of the results and establish a range of validity for the methods. Previously, we presented such a comparison for a set of pure radiative transfer tests (i.e. for fixed, non-evolving density fields). This is the second paper of the Cosmological Radiative Transfer (RT) Comparison Project, in which we compare 9 independent RT codes directly coupled to gasdynamics on 3 relatively simple astrophysical hydrodynamics problems: (5) the expansion of an H II region in a uniform medium; (6) an ionization front (I-front) in a 1/r^2 density profile with a flat core, and (7), the photoevaporation of a uniform dense clump. Results show a broad agreement between the different methods and no big failures, indicating that the participating codes have reached a certain level of maturity and reliability. However, many details still do differ, and virtually every code has showed some shortcomings and has disagreed, in one respect or another, with the majority of the results. This underscores the fact that no method is universal and all require careful testing of the particular features which are most relevant to the specific problem at hand.

We present numerical simulations of how a 120 M$_\odot$ primordial star regulates star formation in nearby cosmological halos at $z \sim$ 20 by photoevaporation. Our models include nine-species primordial chemistry and self-consistent multifrequency conservative transfer of UV photons with all relevant radiative processes. Whether or not new stars form in halos clustered around a Population III star ultimately depends on their core densities and proximity to the star. Diffuse halos with central densities below 2 - 3 cm$^{-3}$ are completely ionized and evaporated anywhere in the cluster. Evolved halos with core densities above 2000 cm$^{-3}$ are impervious to both ionizing and Lyman-Werner flux at most distances from the star and collapse as quickly as they would in its absence. Star formation in halos of intermediate density can be either promoted or suppressed depending on how the I-front remnant shock compresses, deforms and enriches the core with H$_2$. We find that the 120 M$_\odot$ star photodissociates H$_2$ in most halos in the cluster but that catalysis by H- restores it a few hundred kyr after the death of the star, with little effect on star formation. Our models exhibit significant departures from previous one-dimensional spherically-symmetric simulations, which are prone to serious errors due to unphysical geometric focusing effects.

We present numerical simulations of the photoevaporation of cosmological halos clustered around a 120 M$_\odot$ primordial star, confining our study to structures capable of hosting Population III star formation. The calculations include self-consistent multifrequency conservative transfer of UV photons together with nine-species primordial chemistry and all relevant radiative processes. The ultimate fates of these halos varies with central density and proximity to the central source but generally fall into one of four categories. Diffuse halos with central densities below 2 - 3 cm$^{-3}$ are completely ionized and evaporated by the central star anywhere in the cluster. More evolved halo cores at densities above 2000 cm$^{-3}$ are impervious to both ionizing and Lyman-Werner flux at most distances from the star and collapse of their cores proceeds without delay. Radiative feedback in halos of intermediate density can be either positive or negative, depending on how the I-front remnant shock both compresses and deforms the core and enriches it with H$_2$. We find that the 120 M$_\odot$ star photodissociates H$_2$ in most halos within the cluster but that catalysis by H- rapidly restores molecular hydrogen within a few hundred Kyr after the death of the star, with little delay in star formation. Our models exhibit significant departures from previous one-dimensional spherically-symmetric simulations, which are prone to serious errors due to unphysical geometric focusing effects.