Precision measurement of cosmic magnification from 21 cm emitting galaxies

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 unchange

Sleeper end resistance of ballasted railway tracks

This paper describes model tests used to investigate how ballast shoulder width and height contribute to a railway sleeper’s resistance to lateral movement for a range of shoulder widths and heights. Deflection and resistance were measured and photographs taken during the tests.The photographs were analyzed using a digital image correlation technique to identify the zones of ballast surface disturbance, which demonstrated that a bulbed failure volume was mobilized at the ultimate limit state. An idealized three-dimensional failure mechanism is proposed, and resistances are calculated using the limit equilibrium approach. The calculation provides a reliable estimate of the measured resistance. The work identifies the optimum shoulder width and height. The calculations are extended to demonstrate that when a number of sleepers are moved simultaneously, the sleeper end resistance may be one-third less per sleeper than that indicated in tests on an isolated sleeper. Image analysis and limit equilibrium calculations show that this is caused by overlapping of mobilized failure volumes from adjacent sleepers

Efficient approximations of neutrino physics for three-dimensional simulations of stellar core collapse

Neutrino transport in spherically symmetric models of stellar core collapse and bounce has achieved a technically complete level, rewarded by the agreement among independent groups that a multi-dimensional treatment of the fluid-instabilities in the post-bounce phase is indispensable to model supernova explosions. While much effort is required to develop a reliable neutrino transport technique in axisymmetry, we explore neutrino physics approximations and parameterizations for an efficient three-dimensional simulation of the fluid-instabilities in the shock-heated matter that accumulates between the accretion shock and the protoneutron star. We demonstrate the reliability of a simple parameterization scheme in the collapse phase and extend our 3D magneto-hydrodynamical collapse simulations to a preliminary postbounce evolution. The growth of magnetic fields is investigated.Comment: 5 pages, 4 figures, in Proceedings of "Nuclei in the Cosmos IX, Geneva, Jun 25-30", associated movies are displayed at http://www.physik.unibas.ch/~liebend/displa

Precision era of the kinetic Sunyaev-Zeldovich effect: simulations, analytical models and observations and the power to constrain reionization

The kinetic SZ effect, which is the dominant CMB source at arc-minute scales and $\nu \sim 217$ Ghz, probes the ionized gas peculiar momentum up to the epoch of reionization and is a sensitive measure of the reionization history. We ran high resolution self-similar and $\Lambda$CDM hydro simulations and built an analytical model to study this effect. Our model reproduces the $\Lambda$CDM simulation results to several percent accuracy, passes various tests against self-similar simulations, and shows a wider range of applicability than previous analytical models. Our model in its continuous version is free of simulation limitations such as finite simulation box and finite resolution and allows an accurate prediction of the kinetic SZ power spectrum $C_l$. For the WMAP cosmology, we find $l^2C_l/(2\pi)\simeq 0.91 \times 10^{-12} [(1+z_{\rm reion})/10]^{0.34}(l/5000)^{0.23-0.015(z_{\rm reion}-9)}$ for the reionization redshift $6<z_{\rm reion}<20$ and $3000<l<9000$. The corresponding temperature fluctuation is several $\mu$K at these ranges. The dependence of $C_l$ on the reionization history allows an accurate measurement of the reionization epoch. For the Atacama cosmology telescope experiment, $C_l$ can be measured with $\sim 1$ accuracy. $C_l$ scales as $(\Omega_b h)^2 \sigma_8^{4\sim 6}$. Given cosmological parameters, ACT would be able to constrain $z_{\rm reion}$ with several percent accuracy. Some multi-reionization scenarios degenerate in the primary CMB temperature and TE measurement can be distinguished with $\sim 10 \sigma$ confidence.Comment: 14 pages, 7 figures. Accepted by MNRAS. We corrected the primary CMB power spectrum we used. We added discussions about the effects of lensing and relativistic SZ correctio. We withdraw a claim about the patchy reionizatio

The Power Spectrum of the Sunyaev-Zel'dovich Effect

(Abridged) The hot gas in the IGM produces anisotropies in the Cosmic Microwave Background (CMB) through the thermal Sunyaev-Zel'dovich (SZ) effect. The SZ effect is a powerful probe of large-scale structure in the universe and must be carefully subtracted from measurements of the primary CMB anisotropies. We use moving-mesh hydrodynamical simulations to study the 3-dimensional statistics of the gas, and compute the mean comptonization parameter and the angular power spectrum of the SZ fluctuations, for different cosmologies. We compare these results with predictions using the Press-Schechter formalism. We find that the two methods agree approximately, but differ in details. We discuss this discrepancy, and show that resolution limits the reliability of our results to the 200<l<2000 range. For cluster- normalized CDM models, the SZ power spectrum is comparable to the primordial power spectrum around l=2000. We show that groups and filaments (kT<5 keV) contribute about 50% of the SZ power spectrum at l=500. About half of the SZ power spectrum on these scales is produced at redshifts z<0.1, and can thus be detected and removed using existing catalogs of galaxies and X-ray clusters. We discuss the implications of these results for the future MAP and Planck Surveyor missions.Comment: 21 revtex pages, including 2 tables and 12 figures. To appear in PRD. Minor revisions to match accepted version. Also available at http://www.astro.princeton.edu/~refre

FISH: A 3D parallel MHD code for astrophysical applications

FISH is a fast and simple ideal magneto-hydrodynamics code that scales to ~10 000 processes for a Cartesian computational domain of ~1000^3 cells. The simplicity of FISH has been achieved by the rigorous application of the operator splitting technique, while second order accuracy is maintained by the symmetric ordering of the operators. Between directional sweeps, the three-dimensional data is rotated in memory so that the sweep is always performed in a cache-efficient way along the direction of contiguous memory. Hence, the code only requires a one-dimensional description of the conservation equations to be solved. This approach also enable an elegant novel parallelisation of the code that is based on persistent communications with MPI for cubic domain decomposition on machines with distributed memory. This scheme is then combined with an additional OpenMP parallelisation of different sweeps that can take advantage of clusters of shared memory. We document the detailed implementation of a second order TVD advection scheme based on flux reconstruction. The magnetic fields are evolved by a constrained transport scheme. We show that the subtraction of a simple estimate of the hydrostatic gradient from the total gradients can significantly reduce the dissipation of the advection scheme in simulations of gravitationally bound hydrostatic objects. Through its simplicity and efficiency, FISH is as well-suited for hydrodynamics classes as for large-scale astrophysical simulations on high-performance computer clusters. In preparation for the release of a public version, we demonstrate the performance of FISH in a suite of astrophysically orientated test cases.Comment: 27 pages, 11 figure

Halo stochasticity in global clustering analysis

In the present work we study the statistics of haloes, which in the halo model determines the distribution of galaxies. Haloes are known to be biased tracer of dark matter, and at large scales it is usually assumed there is no intrinsic stochasticity between the two fields. Following the work of Seljak & Warren (2004), we explore how correct this assumption is and, moving a step further, we try to qualify the nature of stochasticity. We use Principal Component Analysis applied to the outputs of a cosmological N-body simulation to: (1) explore the behaviour of stochasticity in the correlation between haloes of different masses; (2) explore the behaviour of stochasticity in the correlation between haloes and dark matter. We show results obtained using a catalogue with 2.1 million haloes, from a PMFAST simulation with box size of 1000h^{-1}Mpc. In the relation between different populations of haloes we find that stochasticity is not-negligible even at large scales. In agreement with the conclusions of Tegmark & Bromley (1999) who studied the correlations of different galaxy populations, we found that the shot-noise subtracted stochasticity is qualitatively different from `enhanced' shot noise and, specifically, it is dominated by a single stochastic eigenvalue. We call this the `minimally stochastic' scenario, as opposed to shot noise which is `maximally stochastic'. In the correlation between haloes and dark matter, we find that stochasticity is minimized, as expected, near the dark matter peak (k ~ 0.02 h Mpc^{-1} for a LambdaCDM cosmology) and, even at large scales, it is of the order of 15 per cent above the shot noise. Moreover, we find that the reconstruction of the dark matter distribution is improved when we use eigenvectors as tracers of the bias. [Abridged]Comment: 9 pages, 12 figures. Submitted to MNRA

Gravity and Non-gravity Modes in the VIRMOS-DESCART Weak Lensing Survey

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<l<10^4.Comment: accepted version, minor changes, 18 pages including 6 figure

Spin and orbital excitation spectrum in the Kugel-Khomskii model

We discuss spin and orbital ordering in the twofold orbital degenerate superexchange model in three dimensions relevant to perovskite transition metal oxides. We focus on the particular point on the classical phase diagram where orbital degeneracy is lifted by quantum effects exclusively. Dispersion and damping of the spin and orbital excitations are calculated at this point taking into account their mutual interaction. Interaction corrections to the mean-field order parameters are found to be small. We conclude that quasi-one-dimensional Neel spin order accompanied by the uniform d_{3z^2-r^2}-type orbital ordering is stable against quantum fluctuations.Comment: 4 pages with 3 PS figures, 1 table, RevTeX, accepted to Phys. Rev. B. Rapid Communicatio