Search

We present the first large-scale radiative transfer simulations of cosmic reionization, in a simulation volume of (100/h Mpc)^3, while at the same time capturing the dwarf galaxies which are primarily responsible for reionization. We achieve this by combining the results from extremely large, cosmological, N-body simulations with a new, fast and efficient code for 3D radiative transfer, C^2-Ray. The resulting electron-scattering optical depth is in good agreement with the first-year WMAP polarization data. We show that reionization clearly proceeded in an inside-out fashion, with the high-density regions being ionized earlier, on average, than the voids. Ionization histories of smaller-size (5 to 10 comoving Mpc) subregions exibit a large scatter about the mean and do not describe the global reionization history well. The minimum reliable volume size for such predictions is ~30 Mpc. We derive the power-spectra of the neutral, ionized and total gas density fields and show that there is a significant boost of the density fluctuations in both the neutral and the ionized components relative to the total at arcminute and larger scales. We find two populations of HII regions according to their size, numerous, mid-sized (~10 Mpc) regions and a few, rare, very large regions tens of Mpc in size. We derive the statistical distributions of the ionized fraction and ionized gas density at various scales and for the first time show that both distributions are clearly non-Gaussian. (abridged)

We show how precision lensing measurements can be obtained through the lensing magnification effect in high redshift 21cm emission from galaxies. Normally, cosmic magnification measurements have been seriously complicated by galaxy clustering. With precise redshifts obtained from 21cm emission line wavelength, one can correlate galaxies at different source planes, or exclude close pairs to eliminate such contaminations. We provide forecasts for future surveys, specifically the SKA and CLAR. SKA can achieve percent precision on the dark matter power spectrum and the galaxy dark matter cross correlation power spectrum, while CLAR can measure an accurate cross correlation power spectrum. The neutral hydrogen fraction was most likely significantly higher at high redshifts, which improves the number of observed galaxies significantly, such that also CLAR can measure the dark matter lensing power spectrum. SKA can also allow precise measurement of lensing bispectrum., Comment: 11 pages, 8 figures. Accepted to MNRAS. We deleted two figures and shortened the paper to meet MNRAS's requirement. All main results remain unchanged

The relic neutrinos are expected to acquire a bulk relative velocity with respect to the dark matter at low redshifts, and neutrino wakes are expected to develop downstream of the dark matter halos. We propose a method of measuring the neutrino mass based on this mechanism. This neutrino wake will cause a dipole distortion of the galaxy-galaxy lensing pattern. This effect could be detected by combining upcoming lensing surveys with a low redshift galaxy survey or a 21 cm intensity mapping survey, which can map the neutrino flow field. The data obtained with LSST and Euclid should enable us to make a positive detection if the three neutrino masses are quasidegenerate with each neutrino mass of $\sim$0.1 eV, and a future high precision 21 cm lensing survey would allow the normal hierarchy and inverted hierarchy cases to be distinguished, and even the right-handed Dirac neutrinos may be detectable., 5 pages, 4 figures, v2 matches the published version

All recent numerical simulations agree that stars in the main sequence mass range of 9-40 solar masses do not produce a prompt hydrodynamic ejection of the outer layers after core collapse and bounce. Rather they suggest that stellar core collapse and supernova explosion are dynamically distinct astrophysical events, separated by an unspectacular accretion phase of at least ~40 ms duration. As long as the neutrinospheres remain convectively stable, the explosion dynamics is determined by the neutrons, protons, electrons and neutrinos in the layer of impact-heated matter piling up on the protoneutron star. The crucial role of neutrino transport in this regime has been emphasized in many previous investigations. Here, we search for efficient means to address the role of magnetic fields and fluid instabilities in stellar core collapse and the postbounce phase., Comment: 4 pages, contribution to Nuclei in the Cosmos VIII, Jul. 19-23, submitted to Nucl. Phys. A

We discuss the nature of the possible high-l excess in the Cosmic Microwave Background (CMB) anisotropy power spectrum observed by the Cosmic Background Imager (CBI). We probe the angular structure of the excess in the CBI deep fields and investigate whether it could be due to the scattering of CMB photons by hot electrons within clusters, the Sunyaev-Zeldovich (SZ) effect. We estimate the density fluctuation parameters for amplitude, sigma_8, and shape, Gamma, from CMB primary anisotropy data and other cosmological data. We use the results of two separate hydrodynamical codes for Lambda-CDM cosmologies, consistent with the allowed sigma_8 and Gamma values, to quantify the expected contribution from the SZ effect to the bandpowers of the CBI experiment and pass simulated SZ effect maps through our CBI analysis pipeline. The result is very sensitive to the value of sigma_8, and is roughly consistent with the observed power if sigma_8 ~ 1. We conclude that the CBI anomaly could be a result of the SZ effect for the class of Lambda-CDM concordance models if sigma_8 is in the upper range of values allowed by current CMB and Large Scale Structure (LSS) data., Comment: Accepted by The Astrophysical Journal; 17 pages including 12 color figures. v2 matches accepted version. Additional information at http://www.astro.caltech.edu/~tjp/CBI/

We present a new parallel PM N-body code named PMFAST that is freely available to the public. PMFAST is based on a two-level mesh gravity solver where the gravitational forces are separated into long and short range components. The decomposition scheme minimizes communication costs and allows tolerance for slow networks. The code approaches optimality in several dimensions. The force computations are local and exploit highly optimized vendor FFT libraries. It features minimal memory overhead, with the particle positions and velocities being the main cost. The code features support for distributed and shared memory parallelization through the use of MPI and OpenMP, respectively. The current release version uses two grid levels on a slab decomposition, with periodic boundary conditions for cosmological applications. Open boundary conditions could be added with little computational overhead. We present timing information and results from a recent cosmological production run of the code using a 3712^3 mesh with 6.4 x 10^9 particles. PMFAST is cost-effective, memory-efficient, and is publicly available., Comment: 18 pages, 11 figures

The propagation of cosmological ionization fronts (I-fronts) during reionization is strongly influenced by small-scale structure. Here we summarize our recent attemps to understand the effect of this small-scale structure. We present high resolution cosmological N-body simulations at high-z (z>6) which resolve a wide range of halo mass, from minihalos to clusters of large, rare halos. We also study how minihalos affect I-fronts, through simulations of minihalo photoevaporation. Furthermore, we modify the I-front propagation equations to account for evolving small-scale structure, and incorporate these results into a semi-analytical reionization model. When intergalactic medium clumping and minihalo clustering around sources are included, small-scale structure affects reionization by slowing it down and extending it in time. This helps to explain observations by the Wilkinson Microwave Anisotropy Probe, which imply an early and extended reionization epoch. We also study how source clustering affects the evolution and size of H II regions, finding, in agreement with simulations, that H II regions usually expand, rarely shrinking. Hence, "relic H II regions'' are an exception, rather than the rule. When the suppression of small-mass sources in already-ionized regions by Jeans-mass filtering is accounted for, H II regions are smaller, delaying overlap. We also present a new numerical method for radiative transfer which is fast, efficient, and easily coupled to hydrodynamics and N-body codes, along with sample tests and applications., 6 pages, 6 figures, revised to match the published version. To appear in Proceedings of IAU Colloquium No. 199, 2005, P.R. Williams, C. Shu, and B. Menard, eds

(Abridged) The precision modeling of the thermal Sunyaev Zeldovich (SZ) effect and future high accuracy measurements will provide a powerful way to constrain the thermal history of the universe. In this paper, we combine high resolution adiabatic hydro simulations and analytical models to provide a precision modeling of the gas density weighted temperature $\bar{T}_g$ and the mean SZ Compton $y$ parameter. Simulation artifacts are quantified and corrected and our analytical model passes tests of various simulations and is able to predict $\bar{T}_g$ and $\bar{y}$ to several percent accuracy. For low matter density $\Lambda$CDM cosmology, the present $\bar{T}_g$ is $0.32 (\sigma_8/0.84)^{3.05-0.15\Omega_m}(\Omega_m/0.268)^{1.28-0.2\sigma_8} {\rm keV}$. The mean $y$ parameter is $2.6\times 10^{-6} (\sigma_8/0.84)^{4.1-2\Omega_m}(\Omega_m/0.268)^{1.28-0.2\sigma_8}$. The current upper limit of $y

The PrimevAl Structure Telescope (PAST), will be used to locate and study the era the of the first luminous objects, the epoch of reionization. The first stars ionized the gas around them producing a pattern of ionization that reflects the large scale density structure present at the time. The PAST array will be used to sense and study this ionization, by mapping the brightness of 21 cm neutral Hydrogen Cosmic Background (HCB) at redshift from 6 to 25. The HCB disappears on ionization, allowing the study of large scale structure and of star formation at this very early epoch.

We describe a numerical method to solve the magnetohydrodynamic (MHD) equations. The fluid variables are updated along each direction using the flux conservative, 2nd order, total variation diminishing (TVD), upwind scheme of Jin and Xin. The magnetic field is updated separately in two-dimensional advection-constraint steps. The electromotive force (EMF) is computed in the advection step using the TVD scheme, and this same EMF is used immediately in the constraint step in order to preserve \grad.B=0 without the need to store intermediate fluxes. Operator splitting is used to extend the code to three dimensions, and Runge-Kutta is used to get second order accuracy in time. The advantages of this code are high resolution per grid cell, second order accuracy in space and time, enforcement of the \grad.B=0 constraint to machine precision, no memory overhead, speed, and simplicity. A 3-D Fortran implementation less than 400 lines long is made freely available. We also implemented a fully scalable message-passing parallel MPI version. We present tests of the code on MHD waves and shocks.

We present a new algorithm to rapidly and optimally compute power spectra. This new algorithm is based on a generalization of iterative multigrid, and has computational cost O(N log N), compared to the standard brute force approach which costs O(N^3). The procedure retains this speed on the full sky and for ill-conditioned matrices. It is applicable to galaxy power spectra, CMB, polarization and weak lensing data. We present a mathematical convergence analysis, and performance results.

Weak gravitational lensing provides a direct statistical measure of the dark matter distribution. The variance is easiest to measure, which constrains the degenerate product \sigma_8\Omega^0.6. The degeneracy is broken by measuring the skewness arising from the fact that densities must remain positive, which is not possible when the initially symmetric perturbations become non-linear. Skewness measures the non-linear mass scale, which in combination with the variance measures \Omega directly. We present the first detection of dark matter skewness from the Virmos-Decart survey. We have measured the full three point function, and its projections onto windowed skewness. We separate the lensing mode and the B mode. The lensing skewness is detected for a compensated Gaussian on scales of 5.37 arc minutes to be \kappa^3=1.06+/-0.06x10^-6. The B-modes are consistent with zero at this scale. The variance for the same window function is \kappa^2= 5.32+/-0.62+/-0.98x10^-5, resulting in S_3=375^{+342}_{-124}. Comparing to N-body simulations, we find \Omega_0, Comment: 26 pages 8 figures, accepted by ApJ, minor revisions, corrected eqn (21)

We present a survey of the cosmological applications of the next generation of weak lensing surveys, paying special attention to the computational challenges presented by the number of galaxies, $N_{gal} ~$ 10$^{5}$. We focus on optimal methods with no pixelization and derive a multigrid $P^3M$ algorithm that performs the relevant computations in $O(N_{gal} \log N_{gal})$ time. We test the algorithm by studying three applications of weak lensing surveys - convergence map reconstruction, cluster detection and $E$ and $B$ power spectrum estimation using realistic 1 deg^{2} simulations derived from N-body simulations. The map reconstruction is able to reconstruct large scale features without artifacts. Detecting clusters using only weak lensing is difficult because of line of sight contamination and noise, with low completeness if one desires low contamination of the sample. A power spectrum analysis of the convergence field is more promising and we are able to reconstruct the convergence spectrum with no loss of information down to the smallest scales. The numerical methods used here can be applied to other data sets with same $O(N\log N)$ scaling and can be generalised to a sphere., Submitted to New Astronomy. High resolution version available at http://feynman.princeton.edu/~paddy

We report measurements of anisotropy in the cosmic microwave background radiation over the multipole range l ~ 200 - 3500 with the Cosmic Background Imager based on deep observations of three fields. These results confirm the drop in power with increasing l first reported in earlier measurements with this instrument, and extend the observations of this decline in power out to l \~ 2000. The decline in power is consistent with the predicted damping of primary anisotropies. At larger multipoles, l = 2000 - 3500, the power is 3.1 sigma greater than standard models for intrinsic microwave background anisotropy in this multipole range, and 3.5 sigma greater than zero. This excess power is not consistent with expected levels of residual radio source contamination but, for sigma_8 >~ 1, is consistent with predicted levels due to a secondary Sunyaev-Zeldovich anisotropy. Further observations are necessary to confirm the level of this excess and, if confirmed, determine its origin., Accepted by The Astrophysical Journal; 15 pages including 7 color figures. Additional information at http://www.astro.caltech.edu/~tjp/CBI/

We describe an algorithm for the extraction of the angular power spectrum of an intensity field, such as the cosmic microwave background (CMB), from interferometer data. This new method, based on the gridding of interferometer visibilities in the aperture plane followed by a maximum likelihood solution for bandpowers, is much faster than direct likelihood analysis of the visibilities, and deals with foreground radio sources, multiple pointings, and differencing. The gridded aperture-plane estimators are also used to construct Wiener-filtered images using the signal and noise covariance matrices used in the likelihood analysis. Results are shown for simulated data. The method has been used to determine the power spectrum of the cosmic microwave background from observations with the Cosmic Background Imager, and the results are given in companion papers.

We present a pedagogical review of some of the methods employed in Eulerian computational fluid dynamics (CFD). Fluid mechanics is governed by the Euler equations, which are conservation laws for mass, momentum, and energy. The standard approach to Eulerian CFD is to divide space into finite volumes or cells and store the cell-averaged values of conserved hydro quantities. The integral Euler equations are then solved by computing the flux of the mass, momentum, and energy across cell boundaries. We review both first-order and second-order flux assignment schemes. All linear schemes are either dispersive or diffusive. The nonlinear, second-order accurate total variation diminishing (TVD) approach provides high resolution capturing of shocks and prevents unphysical oscillations. We review the relaxing TVD scheme, a simple and robust method to solve systems of conservation laws like the Euler equations. A 3-D relaxing TVD code is applied to the Sedov-Taylor blast wave test. The propagation of the blast wave is accurately captured and the shock front is sharply resolved. We apply a 3-D self-gravitating hydro code to simulating the formation of blue straggler stars through stellar mergers and present some numerical results. A sample 3-D relaxing TVD code is provided in the appendix., Comment: 23 pages, 12 figures; includes sample 3-D hydro code and new section on stellar mergers; accepted by PASP

We analyze the weak lensing data of the VIRMOS imaging survey using projections (called E and B-modes) of the two independents observed correlation functions. The E-mode contains all the lensing signal, while noise and systematics contribute equally to the E and B modes provided that intrinsic alignment is negligible. The mode separation allows a measurement of the signal with a \sqrt{2} smaller error bars, and a separate channel to test for systematic errors. We apply various transformations, including a spherical harmonic space power spectrum C^E_l and C^B_l, which provides a direct measurement of the projected dark matter distribution for 500, Comment: accepted version, minor changes, 18 pages including 6 figures

We report the first direct observational evidence for the alignment of galaxy spin axes with the local tidal shear field. We measure quantitatively the strength of this directional correlation of disk galaxies from the Tully catalog, with the local shear field reconstructed from Point Source Catalog Redshift survey galaxies. We demonstrate that the null hypothesis of random galaxy alignments relative to the shear frame can be ruled out with more than 99.98% confidence. The observed intrinsic correlation averaged over the censored samples that have detected nonzero signals is measured in terms of the alignment parameter = 0.17 ± 0.04, which includes only statistical errors of the censored data but not the cosmic variance error. The reconstruction procedure is likely to underestimate systematically. Our result is consistent with the linear tidal-torque model and supports the idea that the present galaxy spins may be used as a probe of primordial tidal shear and mass density fields. The intrinsic alignments of galaxy spins dominate over weak gravitational lensing for shallow surveys such as the Sloan Digital Sky Survey, while it should be negligible for deeper surveys at z ~ 1.

The black hole information paradox presumes that quantum field theory in curved space–time can provide unitary propagation from a near-horizon mode to an asymptotic Hawking quantum. Instead of invoking conjectural quantum-gravity effects to modify such an assumption, we propose a self-consistency check. We establish an analogy to Feynman’s analysis of a double-slit experiment. Feynman showed that unitary propagation of the interfering particles, namely ignoring the entanglement with the double-slit, becomes an arbitrarily reliable assumption when the screen upon which the interference pattern is projected is infinitely far away. We argue for an analogous self-consistency check for quantum field theory in curved space–time. We apply it to the propagation of Hawking quanta and test whether ignoring the entanglement with the geometry also becomes arbitrarily reliable in the limit of a large black hole. We present curious results to suggest a negative answer, and we discuss how this loss of naive unitarity in QFT might be related to a solution of the paradox based on the soft-hair-memory effect.

We apply the nonlinear reconstruction method to simulated halo fields. For halo number density $2.77\times 10^{-2}$ $(h^{-1} {\rm Mpc})^{-3}$ at $z=0$, corresponding to the SDSS main sample density, we find the scale where the noise saturates the linear signal is improved to $k\gtrsim0.36\ h {\rm Mpc}^{-1}$, a factor of $2.29$ improvement in scale, or $12$ in number of linear modes. The improvement is less for higher redshift or lower halo density. We expect this to substantially improve the BAO accuracy of dense, low redshift surveys, including the SDSS main sample, 6dFGS and 21cm intensity mapping initiatives., Comment: 10 pages, 7 figures, published on ApJ. Match to published version. Two subsections and two figures are added

Neutral hydrogen clouds are known to exist in the Universe, however their spatial distributions and physical properties are poorly understood. Such missing information can be studied by the new generation Chinese radio telescopes through a blind searching of 21-cm absorption systems. We forecast the capabilities of surveys of 21-cm absorption systems by two representative radio telescopes in China -- Five-hundred-meter Aperture Spherical radio Telescope (FAST) and Tianlai 21-cm cosmology experiment (Tianlai). Facilitated by either the high sensitivity (FAST) or the wide field of view (Tianlai) of these telescopes, more than a thousand 21-cm absorption systems can be discovered in a few years, representing orders of magnitude improvement over the cumulative discoveries in the past half a century., 4 pages. Accepted for publication in Research in Astronomy and Astrophysics

We propose a possibility to improve the current precision measurements on compact binaries. When the orbital axis is almost perpendicular to our line of sight, a pulsar behind its companion can form two strong-lensing images. These images cannot be resolved, but we can use multi-wavelength interferometry to accurately determine the passage through superior conjunction. This method does not depend strongly on the stability of the pulse profile, and applies equally well to both slow and fast pulsars. We discuss the possible improvement this can bring to the bound on stochastic gravitational wave background and to determine black hole spin. We also discuss the possibility of discovering a suitable binary system by the Square Kilometer Array that our method can apply to., 5 pages, 5 figures

So far there is only indirect evidence that the Universe is undergoing an accelerated expansion. The evidence for cosmic acceleration is based on the observation of different objects at different distances, and requires invoking the Copernican cosmological principle, and Einstein's equations of motion. We examine the direct observability using recession velocity drifts (Sandage-Loeb effect) of 21cm hydrogen absorption systems in upcoming radio surveys. This measures the change in velocity of the {\it same} objects separate by a time interval and is a model-independent measure of acceleration. We forecast that for a CHIME-like survey with a decade time span, we can detect the acceleration of a $\Lambda$CDM Universe with $5\sigma$ confidence. This acceleration test requires modest data analysis and storage changes from the normal processing, and cannot be recovered retroactively., Comment: 4 pages, 1 figure. Phys. Rev. Lett. 113, 041303, 2014

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

Among the next generation of radio interferometers are a class of transit telescopes such as the Canadian Hydrogen Intensity Mapping Experiment (CHIME) which are designed for surveying cosmological 21 cm emission. Analysing the data from these instruments is challenging: they are instantaneously very wide-field, map most of the observable sky every day, and require the removal of foregrounds many orders of magnitude brighter than the 21 cm signal.

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

Polarization leakage of foreground synchrotron emission is a critical issue in HI intensity mapping experiments. While the sought-after HI emission is unpolarized, polarized foregrounds such as Galactic and extragalactic synchrotron radiation, if coupled with instrumental impurity, can mimic or overwhelm the HI signals. In this paper we present the methodology for polarization calibration at 700-900 MHz, applied on data obtained from the Green Bank Telescope (GBT). We use astrophysical sources, both polarized and unpolarized sources including quasars and pulsars, as calibrators to characterize the polarization leakage and control systematic effects in our GBT HI intensity mapping project. The resulting fractional errors on polarization measurements on boresight are well controlled to within 0.6%-0.8% of their total intensity. The polarized beam patterns are measured by performing spider scans across both polarized quasars and pulsars. A dominant Stokes I to V leakage feature and secondary features of Stokes I to Q and I to U leakages in the 700-900 MHz frequency range are identified. These characterizations are important for separating foreground polarization leakage from the HI 21 cm signal., Comment: 12 pages, 12 figures, accepted for publication in ApJ

Observations of 21-cm radio emission by neutral hydrogen at redshifts z approximately 0.5 to approximately 2.5 are expected to provide a sensitive probe of cosmic dark energy. This is particularly true around the onset of acceleration at z approximately 1, where traditional optical cosmology becomes very difficult because of the infrared opacity of the atmosphere. Hitherto, 21-cm emission has been detected only to z = 0.24. More distant galaxies generally are too faint for individual detections but it is possible to measure the aggregate emission from many unresolved galaxies in the 'cosmic web'. Here we report a three-dimensional 21-cm intensity field at z = 0.53 to 1.12. We then co-add neutral-hydrogen (H i) emission from the volumes surrounding about 10,000 galaxies (from the DEEP2 optical galaxy redshift survey). We detect the aggregate 21-cm glow at a significance of approximately 4sigma.

Recent implications of results from quantum information theory applied to black holes has led to the confusing conclusions that requires either abandoning the equivalence principle (e.g. the firewall picture), or the no-hair theorem (e.g. the fuzzball picture), or even more unpalatable options. The recent discovery of a pulsar orbiting a black hole opens up new possibilities for tests of theories of gravity. We examine possible observational effects of semiclassical quantum gravity in the vicinity of black holes, as probed by pulsars and event horizon telescope imaging of flares. Pulsar radiation is observable at wavelengths only two orders of magnitude shorter than the Hawking radiation, so precision interferometry of lensed pulsar images may shed light on the quantum gravitational processes and interaction of Hawking radiation with the spacetime near the black hole. This paper discusses the impact on the pulsar radiation interference pattern, which is observable through the modulation index in the foreseeable future, and discusses a possible classical limit of BHC., added figure to illustrate observable consequences

We report a tentative detection of spin correlations in the Tully catalogue of nearby galaxies. We define a simple but nontrivial spin correlation function, and find an analytic estimate of it in the frame of the linear perturbation theory. Then, we present the observed spin correlation signal from the Tully galaxies with error bars. The three dimensional spin correlation turns out to be significant at the 97% confidence level, detected out to a few $h^{-1}$ Mpc. This observed correlation is consistent with the theoretical prediction based on the gravitational instability picture of galaxy formation. An analysis of systematic errors is also presented. The observed strength of correlation may be sufficient to significantly affect blank field of weak lensing searches., Accepted by the ApJ Letters, revised Physical Analysis, LaTex file, 10 pages, 1 ps figure

We present a general linear algorithm for measuring the surface mass density 1-\kappa from the observable reduced shear g=\gamma/(1-\kappa) in the strong lensing regime. We show that in general, the observed polarization field can be decomposed into ``electric'' and ``magnetic'' components, which have independent and redundant solutions, but perfectly orthogonal noise properties. By combining these solutions, one can increase the signal-to-noise ratio by \sqrt{2}. The solutions allow dynamic optimization of signal and noise, both in real and Fourier space (using arbitrary smoothing windows). Boundary conditions have no effect on the reconstructions, apart from its effect on the signal-to-noise. Many existing reconstruction techniques are recovered as special cases of this framework. The magnetic solution has the added benefit of yielding the global and local parity of the reconstruction in a single step., Comment: final accepted version for ApJL

Using the cosmic virial theorem, Press-Schechter analysis and numerical simulations, we compute the expected X-ray background (XRB) from the diffuse IGM with the clumping factor expected from gravitational shock heating. The predicted fluxes and temperatures are excluded from the observed XRB. The predicted clumping can be reduced by entropy injection. The required energy is computed from the two-point correlation function, as well as from Press-Schechter formalisms. The minimal energy injection of 1 keV/nucleon excludes radiative or gravitational heating as a primary energy source. We argue that the intergalactic medium (IGM) must have been heated through violent processes such as massive supernova bursts. If the heating proceeded through supernova explosions, it likely proceeded in bursts which may be observable in high redshift supernova searches. Within our model we reproduce the observed cluster luminosity-temperature relation with energy injection of 1 keV/nucleon if this injection is assumed to be uncorrelated with the local density. These parameters predict that the diffuse IGM soft XRB has a temperature of ~1 keV with a flux near 10 keV/cm^2 s str keV, which may be detectable in the near future.

The standard unbroken electroweak theory is known to erase baryon number. The baryon number symmetry can be restored in the core of a neutron star as its density diverges via gravitational instability due to a binary merger event. We argue that for certain double Higgs models with discrete symmetries, this process may result in an expanding self-sustained burning front which would convert the entire neutron matter into radiation. This process would release ~10^{54} ergs in electromagnetic radiation over ~10^{-4} sec, with negligible baryonic contamination. The resulting fireball would have all the properties necessary to produce a gamma-ray burst as a result of its interaction with ambient interstellar gas. The subsequent Higgs decay would produce a millisecond burst of ~10^{52} ergs in ~100 GeV neutrinos which might be observable. The above mechanism may have also caused electroweak baryogenesis in the early universe, giving rise to the observed matter-antimatter asymmetry today., 15 pages including 1 figure

An algorithm for simulating self-gravitating cosmological astrophysical fluids is presented. The advantages include a large dynamic range, parallelizability, high resolution per grid element and fast execution speed. The code is based on a finite volume flux conservative Total-Variation-Diminishing (TVD) scheme for the shock capturing hydro, and an iterative multigrid solver for the gravity. The grid is a time dependent field, whose motion is described by a generalized potential flow. Approximately constant mass per cell can be obtained, providing all the advantages of a Lagrangian scheme. The grid deformation combined with appropriate limiting and smoothing schemes guarantees a regular and well behaved grid geometry, where nearest neighbor relationships remain constant. The full hydrodynamic fluid equations are implemented in the curvilinear moving grid, allowing for arbitrary fluid flow relative to the grid geometry. This combination retains all the advantages of the grid based schemes including high speed per fluid element and a rapid gravity solver. The current implementation is described, and empirical simulation results are presented. Accurate execution speed calculations are given in terms of floating point operations per time step per grid cell. This code is freely available to the community., Comment: 53 pages including 14 figures, submitted to ApJS

The string-dominated universe locally resembles an open universe, and fits dynamical measures of power spectra, cluster abundances, redshift distortions, lensing constraints, luminosity and angular diameter distance relations and microwave background observations. We show examples of networks which might give rise to recent string-domination without requiring any fine-tuned parameters. We discuss how future observations can distinguish this model from other cosmologies., Comment: 17 pages including 4 figures, of which one is in color

We show how the virial theorem can be applied to the hot gas in clusters of galaxies to obtain a yardstick, which could then be used to determine cosmological parameters. This yardstick relies on the assumptions of hydrostatic equilibrium and that the gas fraction is approximately constant. The constancy is checked empirically from a local population of clusters. By using the observed parameters consisting of temperature, surface brightness and radial profile $\beta$, one can calculate the expected core radius. Comparing it to the observed angular size, one can in principle calibrate the cosmological deceleration parameter $q_0$. We test this method on a small sample of 6 clusters, and show its promise and accuracy. The preliminary implications would be to suggest $q_0 \approx 0.85\pm 0.29$ with $1-\sigma$ statistical error bars, with several systematic uncertainties remaining. Taken at face value, this would argue against a cosmological constant. The method is robust to errors in the measurement of the core radius as long as the product of the central density and the core radius squared $\rho_0 r_c^2$ are well determined. New lensing and X-ray data can dramatically improve on the statistics.

We present a new technique to measure neutrino masses using their flow field relative to dark matter. Present day streaming motions of neutrinos relative to dark matter and baryons are several hundred km/s, comparable with their thermal velocity dispersion. This results in a unique dipole anisotropic distortion of the matter-neutrino cross power spectrum, which is observable through the dipole distortion in the cross correlation of different galaxy populations. Such a dipole vanishes if not for this relative velocity and so it is a clean signature for neutrino mass. We estimate the size of this effect and find that current and future galaxy surveys may be sensitive to these signature distortions., 5 pages, 4 figures, 1 table, version2 with supplemental materials, v3 matches the published version

We propose that interstellar extreme scattering events, usually observed as pulsar scintillations, may be caused by a coherent agent rather than the usually assumed turbulence of $\rm H_2$ clouds. We find that the penetration of a flux of ionizing, positively charged strangelets or quark nuggets into a dense interstellar hydrogen cloud may produce ionization trails. Depending on the specific nature and energy of the incoming droplets, diffusive propagation or even capture in the cloud are possible. As a result, enhanced electron densities may form and constitute a lens-like scattering screen for radio pulsars and possibly for quasars., 5 pages, 2 figures

The large-scale distribution of neutral hydrogen in the Universe will be luminous through its 21 cm emission. Here, for the first time, we use the auto-power spectrum of 21 cm intensity fluctuations to constrain neutral hydrogen fluctuations at z~0.8. Our data were acquired with the Green Bank Telescope and span the redshift range 0.6 < z < 1 over two fields totalling ~41 deg. sq. and 190 h of radio integration time. The dominant synchrotron foregrounds exceed the signal by ~10^3, but have fewer degrees of freedom and can be removed efficiently. Even in the presence of residual foregrounds, the auto-power can still be interpreted as an upper bound on the 21 cm signal. Our previous measurements of the cross-correlation of 21 cm intensity and the WiggleZ galaxy survey provide a lower bound. Through a Bayesian treatment of signal and foregrounds, we can combine both fields in auto- and cross-power into a measurement of Omega_HI b_HI = [0.62^{+0.23}_{-0.15}] * 10^{-3} at 68% confidence with 9% systematic calibration uncertainty, where Omega_HI is the neutral hydrogen (HI) fraction and b_HI is the HI bias parameter. We describe observational challenges with the present data set and plans to overcome them., 5 pages, 3 figures. v2 as published; MNRASL (2013)

We use VLBI imaging of the interstellar scattering speckle pattern associated with the pulsar PSR 0834+06 to measure the astrometric motion of its emission. The ~ 5AU interstellar baselines, provided by interference between speckles spanning the scattering disk, enable us to detect motions with sub nanoarcsecond accuracy. We measure a small pulse deflection of ~18+/-2 km (not including geometric uncertainties), which is 100 times smaller than the native resolution of this interstellar interferometer. This implies that the emission region is small, and at an altitude of a few hundred km, with the exact value depending on field geometry. This is substantially closer to the star than to the light cylinder. Future VLBI measurements can improve on this finding. This new regime of ultra-precise astrometry may enable precision parallax distance determination of pulsar binary displacements., Comment: 5 pages, corrected type in arxiv abstract

This paper presents CUBEP3M, a publicly-available high performance cosmological N-body code and describes many utilities and extensions that have been added to the standard package. These include a memory-light runtime SO halo finder, a non-Gaussian initial conditions generator, and a system of unique particle identification. CUBEP3M is fast, its accuracy is tuneable to optimize speed or memory, and has been run on more than 27,000 cores, achieving within a factor of two of ideal weak scaling even at this problem size. The code can be run in an extra-lean mode where the peak memory imprint for large runs is as low as 37 bytes per particles, which is almost two times leaner than other widely used N-body codes. However, load imbalances can increase this requirement by a factor of two, such that fast configurations with all the utilities enabled and load imbalances factored in require between 70 and 120 bytes per particles. CUBEP3M is well designed to study large scales cosmological systems, where imbalances are not too large and adaptive time-stepping not essential. It has already been used for a broad number of science applications that require either large samples of non-linear realizations or very large dark matter N-body simulations, including cosmological reionization, halo formation, baryonic acoustic oscillations, weak lensing or non-Gaussian statistics. We discuss the structure, the accuracy, known systematic effects and the scaling performance of the code and its utilities, when applicable., Comment: 20 pages, 17 figures, added halo profiles, updated to match MNRAS accepted version

We present the largest-volume (425 Mpc/h=607 Mpc on a side) full radiative transfer simulation of cosmic reionization to date. We show that there is significant additional power in density fluctuations at very large scales. We systematically investigate the effects this additional power has on the progress, duration and features of reionization, as well as on selected reionization observables. We find that comoving simulation volume of ~100 Mpc/h per side is sufficient for deriving a convergent mean reionization history, but that the reionization patchiness is significantly underestimated. We use jackknife splitting to quantify the convergence of reionization properties with simulation volume for both mean-density and variable-density sub-regions. We find that sub-volumes of ~100 Mpc/h per side or larger yield convergent reionization histories, except for the earliest times, but smaller volumes of ~50 Mpc/h or less are not well converged at any redshift. Reionization history milestones show significant scatter between the sub-volumes, of Delta z=0.6-1 for ~50 Mpc/h volumes, decreasing to Delta z=0.3-0.5 for ~100 Mpc/h volumes, and $\Delta z$~0.1 for ~200 Mpc/h volumes. If we only consider mean-density sub-regions the scatter decreases, but remains at Delta z~0.1-0.2 for the different size sub-volumes. Consequently, many potential reionization observables like 21-cm rms, 21-cm PDF skewness and kurtosis all show good convergence for volumes of ~200 Mpc/h, but retain considerable scatter for smaller volumes. In contrast, the three-dimensional 21-cm power spectra at large scales (k, Comment: 20 pages, 19 figures (most in color)

We present a momentum-space solution for the time-dependent Schr{umlt o}dinger equation of a realistic hydrogen atom in a strong laser pulse. The method can integrate the nonperturbative system to thousands of optical cycles, previously thought not yet feasible. A hydrogen atom under a 285-nm (4.357-eV), 97.5-TW/cm{sup 2} short pulse is found to produce two prominent photoelectron peaks above threshold, and significant intermediate bound-state resonance with the 4{ital p} state. With a 248-nm (5-eV) probe laser pulse that is ten times longer, the resonant structure is explored. {copyright} {ital 1996 The American Physical Society.}

We present a general formalism for computing Cosmic Background Radiation (CBR) and density fluctuations in open models with stiff sources. We find analytic Green's functions for the linearized Einstein equations in the presence of stiff sources and use this formalism to estimate the amplitude and harmonic spectrum of microwave background fluctuations produced by topological defects in an open universe. Unlike inflationary models that predict a flat universe and a spectrum of CBR fluctuations that is enhanced at large angular scales, defect models predict that CBR fluctuations are suppressed on angular scales larger than that subtended by the curvature scale. In an $\Omega = 0.2 - 0.4$ universe, these models, when normalized to the amplitude of CBR fluctuations observed by COBE, require a moderate bias factor, $2-3$, to be compatible with the observed fluctuations in galaxy counts. In these models, accurate predictions can be made which are testable through CBR experiments in the near future. A CBR measurement of $\Omega$ would then be possible, up to the limit imposed by cosmic variance. We discuss some of the philosophical implications of an open model and propose a solution to the flatness problem., Comment: 35 RevTeX 3.0 pages without 9 figures. Full postscript file including figures available by anonymous ftp from astro.princeton.edu:/upen/Papers/open.ps, submitted to Phys Rev D

We develop a finite dimensional matrix method to treat wavepacket propagation in potential scattering problems. The method is useful for general quantum mechanical problems. The scattering dynamics of a minimum uncertainty wavepacket in the Coulomb potential is presented.

We consider the generic scenario of dark energy which arises through the latent heat of a hidden sector first order cosmological phase transition. This field could account for the extra radiation degree of freedom suggested by the CMB. We present the bubble nucleation solution for the viscous limit. The decay rate of the field is constrained by published KSZ data, and may be an explanation of current excess ISW correlations. Cross correlation of current and future surveys can further constrain or test the parameter space. The decay model is plausibly in the observable range, and avoids anthropic problems. This class of models is not well constrained by the popular dark energy figure of merit., 4 pages, 4 figures