Toward an Improved Analytical Description of Lagrangian Bias

We carry out a detailed numerical investigation of the spatial correlation function of the initial positions of cosmological dark matter halos. In this Lagrangian coordinate system, which is especially useful for analytic studies of cosmological feedback, we are able to construct cross-correlation functions of objects with varying masses and formation redshifts and compare them with a variety of analytical approaches. For the case in which both formation redshifts are equal, we find good agreement between our numerical results and the bivariate model of Scannapieco & Barkana (2002; SB02) at all masses, redshifts, and separations, while the model of Porciani et al. (1998) does well for all parameters except for objects with different masses at small separations. We find that the standard mapping between Lagrangian and Eulerian bias performs well for rare objects at all separations, but fails if the objects are highly-nonlinear (low-sigma) peaks. In the Lagrangian case in which the formation redshifts differ, the SB02 model does well for all separations and combinations of masses, apart from a discrepancy at small separations in situations in which the smaller object is formed earlier and the difference between redshifts or masses is large. As this same limitation arises in the standard approach to the single-point progenitor distribution developed by Lacey & Cole (1993), we conclude that a more complete understanding of the progenitor distribution is the most important outstanding issue in the analytic modeling of Lagrangian bias.Comment: 22 pages, 8 figures, ApJ, in pres

A Test of the Collisional Dark Matter Hypothesis from Cluster Lensing

Spergel & Steinhardt proposed the possibility that the dark matter particles are self-interacting, as a solution to two discrepancies between the predictions of cold dark matter models and the observations: first, the observed dark matter distribution in some dwarf galaxies has large, constant-density cores, as opposed to the predicted central cusps; and second, small satellites of normal galaxies are much less abundant than predicted. The dark matter self-interaction would produce isothermal cores in halos, and would also expel the dark matter particles from dwarfs orbiting within large halos. However, another inevitable consequence of the model is that halos should become spherical once most particles have interacted. Here, I rule out this model by the fact that the innermost regions of dark matter halos in massive clusters of galaxies are elliptical, as shown by gravitational lensing and other observations. The absence of collisions in the lensing cores of massive clusters implies that any dark matter self-interaction is too weak to have affected the observed density profiles in the dark-matter dominated dwarf galaxies, or to have eased the destruction of dwarf satellites in galactic halos. If $s_x$ is the cross section and $m_x$ the mass of the dark matter particle, then s_x/m_x < 10^{-25.5} \cm^2/\gev.Comment: to appear in ApJ, January 1 200

Mobile radio propagation prediction using ray tracing methods

The basic problem is to solve the two-dimensional scalar Helmholtz equation for a point source (the antenna) situated in the vicinity of an array of scatterers (such as the houses and any other relevant objects in 1 square km of urban environment). The wavelength is a few centimeters and the houses a few metres across, so there are three disparate length scales in the problem. The question posed by BT concerned ray counting on the assumptions that: (i) rays were subject to a reflection coefficient of about 0.5 when bouncing off a house wall and (ii) that diffraction at corners reduced their energy by 90%. The quantity of particular interest was the number of rays that need to be accounted for at any particular point in order for those neglected to only contribute 10% of the field at that point; a secondary question concerned the use of rays to predict regions where the field was less than 1% of that in the region directly illuminated by the antenna. The progress made in answering these two questions is described in the next two sections and possibly useful representations of the solution of the Helmholtz equations in terms other than rays are given in the final section

Triggering the Formation of Halo Globular Clusters with Galaxy Outflows

We investigate the interactions of high-redshift galaxy outflows with low-mass virialized (Tvir < 10,000K) clouds of primordial composition. While atomic cooling allows star formation in larger primordial objects, such "minihalos" are generally unable to form stars by themselves. However, the large population of high-redshift starburst galaxies may have induced widespread star formation in these objects, via shocks that caused intense cooling both through nonequilibrium H2 formation and metal-line emission. Using a simple analytic model, we show that the resulting star clusters naturally reproduce three key features of the observed population of halo globular clusters (GCs). First, the 10,000 K maximum virial temperature corresponds to the ~ 10^6 solar mass upper limit on the stellar mass of GCs. Secondly, the momentum imparted in such interactions is sufficient to strip the gas from its associated dark matter halo, explaining why GCs do not reside in dark matter potential wells. Finally, the mixing of ejected metals into the primordial gas is able to explain the ~ 0.1 dex homogeneity of stellar metallicities within a given GC, while at the same time allowing for a large spread in metallicity between different clusters. To study this possibility in detail, we use a simple 1D numerical model of turbulence transport to simulate mixing in cloud-outflow interactions. We find that as the shock shears across the side of the cloud, Kelvin-Helmholtz instabilities arise, which cause mixing of enriched material into > 20% of the cloud. Such estimates ignore the likely presence of large-scale vortices, however, which would further enhance turbulence generation. Thus quantitative mixing predictions must await more detailed numerical studies.Comment: 21 pages, 11 figures, Apj in pres

Transfer of autocollimator calibration for use with scanning gantry profilometers for accurate determination of surface slope and curvature of state of the art x ray mirrors

X ray optics, desired for beamlines at free electron laser and diffraction limited storage ring x ray light sources, must have almost perfect surfaces, capable of delivering light to experiments without significant degradation of brightness and coherence. To accurately characterize such optics at an optical metrology lab, two basic types of surface slope profilometers are used the long trace profilers LTPs and nanometer optical measuring NOM like angular deflectometers, based on electronic autocollimator AC ELCOMAT 3000. The inherent systematic errors of the instrument s optical sensors set the principle limit to their measuring performance. Where autocollimator of a NOM like profiler may be calibrated at a unique dedicated facility, this is for a particular configuration of distance, aperture size, and angular range that does not always match the exact use in a scanning measurement with the profiler. Here we discuss the developed methodology, experimental set up, and numerical methods of transferring the calibration of one reference AC to the scanning AC of the Optical Surface Measuring System OSMS , recently brought to operation at the ALS Xray Optics Laboratory. We show that precision calibration of the OSMS performed in three steps, allows us to provide high confidence and accuracy low spatial frequency metrology and not print into measurements the inherent systematic error of tool in use. With the examples of the OSMS measurements with a state of the art x ray aspherical mirror, available from one of the most advanced vendors of X ray optics, we demonstrate the high efficacy of the developed calibration procedure. The results of our work are important for obtaining high reliability data, needed for sophisticated numerical simulations of beamline performance and optimization of beamline usage of the optics. This work was supported by the U. S. Department of Energy under contract number DE AC02 05CH1123

Evolution of perturbations in distinct classes of canonical scalar field models of dark energy

Dark energy must cluster in order to be consistent with the equivalence principle. The background evolution can be effectively modelled by either a scalar field or by a barotropic fluid.The fluid model can be used to emulate perturbations in a scalar field model of dark energy, though this model breaks down at large scales. In this paper we study evolution of dark energy perturbations in canonical scalar field models: the classes of thawing and freezing models.The dark energy equation of state evolves differently in these classes.In freezing models, the equation of state deviates from that of a cosmological constant at early times.For thawing models, the dark energy equation of state remains near that of the cosmological constant at early times and begins to deviate from it only at late times.Since the dark energy equation of state evolves differently in these classes,the dark energy perturbations too evolve differently. In freezing models, since the equation of state deviates from that of a cosmological constant at early times, there is a significant difference in evolution of matter perturbations from those in the cosmological constant model.In comparison, matter perturbations in thawing models differ from the cosmological constant only at late times. This difference provides an additional handle to distinguish between these classes of models and this difference should manifest itself in the ISW effect.Comment: 11 pages, 6 figures, accepted for publication in Phys. Rev.

The Axis-Ratio Distribution of Galaxy Clusters in the SDSS-C4 Catalog as a New Cosmological Probe

We analyze the C4 catalog of galaxy clusters from the Sloan Digital Sky Survey (SDSS) to investigate the axis-ratio distribution of the projected two dimensional cluster profiles. We consider only those objects in the catalog whose virial mass is close to 10^{14}h^{-1}M_{sun}, with member galaxies within the scale radius 1000 kpc. The total number of such objects turns out to be 336. We also derive a theoretical distribution by incorporating the effect of projection onto the sky into the analytic formalism proposed recently by Lee, Jing, & Suto. The theoretical distribution of the cluster axis-ratios is shown to depend on the amplitude of the linear power spectrum (sigma_8) as well as the density parameter (Omega_{m}). Finally, fitting the observational data to the analytic distribution with Omega_{m} and sigma_{8} as two adjustable free parameters, we find the best-fitting value of sigma_{8}=(1.01 +/- 0.09)(Omega_{m}/0.6)^{(0.07 +/- 0.02) +0.1 Omega_{m}}$. It is a new sigma_{8}-Omega_{m} relation, different from the previous one derived from the local abundance of X-ray clusters. We expect that the axis-ratio distribution of galaxy clusters, if combined with the local abundance of clusters, may put simultaneous constraints on sigma_{8} and Omega_{m}.Comment: accepted for publication in ApJ, 17 pages, 3 figures, improved analysis, more discussion on the validity and the caveats of the mode

Double Distribution of Dark Matter Halos with respect to Mass and Local Overdensity

We present a double distribution function of dark matter halos, with respect to both object mass and local over- (or under-) density. This analytical tool provides a statistical treatment of the properties of matter surrounding collapsed objects, and can be used to study environmental effects on hierarchical structure formation. The size of the "local environment" of a collapsed object is defined to depend on the mass of the object. The Press-Schechter mass function is recovered by integration of our double distribution over the density contrast. We also present a detailed treatment of the evolution of overdensities and underdensities in Einstein-deSitter and flat LCDM universes, according to the spherical evolution model. We explicitly distinguish between true and linearly extrapolated overdensities and provide conversion relations between the two quantities.Comment: 25 pages, 10 figures, comments welcom

Probing the structure of the cold dark matter halo with ancient mica

Mica can store (for >1 Gy) etchable tracks caused by atoms recoiling from WIMPs. Ancient mica is a directional detector despite the complex motions it makes with respect to the WIMP "wind". We can exploit the properties of directionality and long integration time to probe for structure in the dark matter halo of our galaxy. We compute a sample of possible signals in mica for a plausible model of halo structure.Comment: 7 pages, 2 figure

The Robustness of Dark Matter Density Profiles in Dissipationless Mergers

We present a comprehensive series of dissipationless N-body simulations to investigate the evolution of density distribution in equal-mass mergers between dark matter (DM) halos and multicomponent galaxies. The DM halo models are constructed with various asymptotic power-law indices ranging from steep cusps to core-like profiles and the structural properties of the galaxy models are motivated by the LCDM paradigm of structure formation. The adopted force resolution allows robust density profile estimates in the inner ~1% of the virial radii of the simulated systems. We demonstrate that the central slopes and overall shapes of the remnant density profiles are virtually identical to those of the initial systems suggesting that the remnants retain a remarkable memory of the density structure of their progenitors, despite the relaxation that accompanies merger activity. We also find that halo concentrations remain approximately constant through hierarchical merging involving identical systems and show that remnants contain significant fractions of their bound mass well beyond their formal virial radii. These conclusions hold for a wide variety of initial asymptotic density slopes, orbital energies, and encounter configurations, including sequences of consecutive merger events, simultaneous mergers of severals ystems, and mergers of halos with embedded cold baryonic components in the form of disks, spheroids, or both. As an immediate consequence, the net effect of gas cooling, which contracts and steepens the inner density profiles of DM halos, should be preserved through a period of dissipationless major merging. Our results imply that the characteristic universal shape of DM density profiles may be set early in the evolution of halos.Comment: Accepted for publication in ApJ, 20 pages, 10 figures, LaTeX (uses emulateapj.cls