Source author record

Joseph Smidt

Joseph Smidt appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

44works
7topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

44 published item(s)

preprint2019arXiv

Outflows from inflows: the nature of Bondi-like accretion

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ρ}/\partialθ|$ 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$.

preprint2016arXiv

Ab Initio Cosmological Simulations of CR7 as an Active Black Hole

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-$α$ 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-$α$ 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 $μ$Jy observations and would discriminate this scenario from one where the luminosity is driven by Population III stars.

preprint2016arXiv

Axion decay and anisotropy of near-IR extragalactic background light

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.

preprint2015arXiv

Emission from Pair-Instability Supernovae with Rotation

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.

preprint2015arXiv

Planck Trispectrum Constraints on Primordial Non-Gaussianity at Cubic Order

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 $τ_{\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 $τ_{\rm{NL}}=0.3 \pm 0.9 \times 10^4$ and $g_{\rm{NL}}=-1.2 \pm 2.8 \times 10^5$. If $τ_{\rm NL}=0$, we find $g_{\rm NL}= -1.3\pm 1.8 \times 10^5$.

preprint2015arXiv

Population III Hypernovae

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.

preprint2015arXiv

The First Gamma-Ray Bursts in the Universe

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},γ} =$ 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.

preprint2015arXiv

Ultraviolet Luminosity Density of the Universe During the Epoch of Reionization

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 $μ$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 ρ_{\rm UV} = 27.4^{+0.2}_{-1.2}$ erg s$^{-1}$ Hz$^{-1}$ Mpc$^{-3}$ $(1σ)$. 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.

preprint2014arXiv

Cross-correlating Sunyaev-Zel'dovich and Weak Lensing Maps

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 (κ_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_{κy}(κ_s,y_s) in terms of individual one-point PDFs (p_κ(κ_s), p_y(y_s)) and the relevant bias functions (b_κ(κ_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.

preprint2014arXiv

Cross-Correlation of Near and Far-Infrared Background Anisotropies as Traced by Spitzer and Herschel

We present the cross-correlation between the far-infrared background fluctuations as measured with the Herschel Space Observatory at 250, 350, and 500 μm and the near-infrared background fluctuations with Spitzer Space Telescope at 3.6 μ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 μm to 500 μ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 μm, and (c) intra-halo light, or diffuse stars in dark matter halos, that likely dominates fluctuations at 3.6 μm. The model is able to reasonably reproduce the auto correlations at each of the far-IR wavelengths and at 3.6 μ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 μ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 μm to 500 μm comes from the fact that the galaxies that are primarily contributing to 500 μm fluctuations peak at a higher redshift than those at 250 μm.

preprint2014arXiv

Finding the First Cosmic Explosions. IV. 90 - 140 M$_{\odot}$ Pair-Instability Supernovae

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.

preprint2014arXiv

On the Origin of Near-Infrared Extragalactic Background Light Anisotropy

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.

preprint2014arXiv

Pair-Instability Supernovae in the Local Universe

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

preprint2013arXiv

Constraints on Spatial Variations in the Fine-Structure constant from Planck

We use the Cosmic Microwave Background (CMB) anisotropy data from Planck to constrain the spatial fluctuations of the fine-structure constant α. Through Thompson scattering of CMB photons, spatial anisotropies of α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 δα/α_0= (1.34 +/- 5.82) x 10^-2 at the 95% confidence level with respect to the standard value α_0. We find no evidence for a spatially varying αat a redshift of 10^3.

preprint2013arXiv

Finding the First Cosmic Explosions I: Pair-Instability Supernovae

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.

preprint2013arXiv

Finding the First Cosmic Explosions. III. Pulsational Pair-Instability Supernovae

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.

preprint2013arXiv

New Approaches to Probing Minkowski Functionals

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.

preprint2013arXiv

Supermassive Population III Supernovae and the Birth of the First Quasars

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.

preprint2013arXiv

The Biggest Explosions in the Universe

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.

preprint2013arXiv

The Biggest Explosions in the Universe. II

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.

preprint2013arXiv

The Supernova that Destroyed a Protogalaxy: Prompt Chemical Enrichment and Supermassive Black Hole Growth

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.

preprint2012arXiv

A measurement of the intrahalo light fraction with near-infrared background anisotropies

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.

preprint2012arXiv

Statistical Properties of Thermal Sunyaev-Zel'dovich Maps

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.

preprint2012arXiv

The Impact of Assuming Flatness in the Determination of Neutrino Properties from Cosmological Data

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_ν, and the curvature of the universe Ω_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 Ω_k, with N_{eff} = 4.03 +/- 0.45 and 10^3 Ω_k = -4.46 +/- 5.24. Thus, even when Ω_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_ν< 0.45 eV to \sum m_ν< 0.95 eV. Thus, the impact of assuming flatness in neutrino cosmology is significant and an essential consideration with future experiments.

preprint2012arXiv

The Near-IR Background Intensity and Anisotropies During The Epoch of Reionization

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

preprint2011arXiv

CANDELS: The Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey

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σ 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σ 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.

preprint2011arXiv

CANDELS: The Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey - The Hubble Space Telescope Observations, Imaging Data Products and Mosaics

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.

preprint2011arXiv

Constraints on the SZ Power Spectrum on Degree Angular Scales in WMAP Data

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.

preprint2011arXiv

From Weak Lensing to non-Gaussianity via Minkowski Functionals

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.

preprint2011arXiv

Lensed galaxies in CANDELS

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.

preprint2011arXiv

Measuring Light from the Epoch of Reionization with CIBER, the Cosmic Infrared Background Experiment

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.

preprint2011arXiv

The Morphology of the Thermal Sunyaev-Zel'dovich Sky

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.

preprint2010arXiv

A Constraint On the Integrated Mass Power Spectrum out to z = 1100 from Lensing of the Cosmic Microwave Background

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.

preprint2010arXiv

A Measurement of Cubic-Order Primordial Non-Gaussianity (g_{NL} and τ_{NL}) With WMAP 5-Year Data

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 τ_{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 < τ_{\rm NL}/10^4 < 3.3 improving the previous COBE-based limit on τ_{\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 τ_{\rm Nl}/(6f_{\rm NL}/5)^2 ranges from -3 to 21 at the 95% confidence level.

preprint2010arXiv

CMB Constraints on Primordial non-Gaussianity from the Bispectrum (f_{NL}) and Trispectrum (g_{NL} and τ_{NL}) and a New Consistency Test of Single-Field Inflation

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 τ_{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 τ_{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} = τ_{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 τ_{NL}, and g_{NL} and find -7.4 < g_{NL}/10^5 < 8.2 and -0.6 < τ_{NL}/10^4 < 3.3 improving the previous COBE-based limit on τ_{NL} < 10^8 nearly four orders of magnitude with WMAP.

preprint2010arXiv

Higher-order Statistics of Weak Lensing Shear and Flexion

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.

preprint2010arXiv

Primordial Non-Gaussianity from a Joint Analysis of Cosmic Microwave Background Temperature and Polarization

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.

preprint2010arXiv

The Herschel-SPIRE Legacy Survey (HSLS): the scientific goals of a shallow and wide submillimeter imaging survey with SPIRE

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.

preprint2010arXiv

The Impact of Secondary non-Gaussianities in the CMB on Cosmological Parameter Estimation

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.

preprint2009arXiv

A Measurement of Primordial Non-Gaussianity Using WMAP 5-Year Temperature Skewness Power Spectrum

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σlevel.

preprint2009arXiv

Cosmological Radiative Transfer Comparison Project II: The Radiation-Hydrodynamic Tests

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.

preprint2009arXiv

New Optimised Estimators for the Primordial Trispectrum

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.

preprint2009arXiv

Non-Gaussianity in WMAP Data Due to the Correlation of CMB Lensing Potential with Secondary Anisotropies

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.

preprint2009arXiv

Parameter Space of General Gauge Mediation

We study a subspace of General Gauge Mediation (GGM) models which generalize models of gauge mediation. We find superpartner spectra that are markedly different from those of typical gauge and gaugino mediation scenarios. While typical gauge mediation predictions of either a neutralino or stau next-to-lightest supersymmetric particle (NLSP) are easily reproducible with the GGM parameters, chargino and sneutrino NLSPs are generic for many reasonable choices of GGM parameters.