PACS and SPIRE photometer maps of M 33: First results of the HERschel M 33 Extended Survey (HERM33ES)

Context. Within the framework of the HERM33ES key program, we are studying the star forming interstellar medium in the nearby, metal-poor spiral galaxy M 33, exploiting the high resolution and sensitivity of Herschel. Aims. We use PACS and SPIRE maps at 100, 160, 250, 350, and 500 μm wavelength, to study the variation of the spectral energy distributions (SEDs) with galacto-centric distance. Methods. Detailed SED modeling is performed using azimuthally averaged fluxes in elliptical rings of 2 kpc width, out to 8 kpc galacto-centric distance. Simple isothermal and two-component grey body models, with fixed dust emissivity index, are fitted to the SEDs between 24 μm and 500 μm using also MIPS/Spitzer  data, to derive first estimates of the dust physical conditions. Results. The far-infrared and submillimeter maps reveal the branched, knotted spiral structure of M 33. An underlying diffuse disk is seen in all SPIRE maps (250–500 μm). Two component fits to the SEDs agree better than isothermal models with the observed, total and radially averaged flux densities. The two component model, with β fixed at 1.5, best fits the global and the radial SEDs. The cold dust component clearly dominates; the relative mass of the warm component is less than 0.3% for all the fits. The temperature of the warm component is not well constrained and is found to be about 60 K ± 10 K. The temperature of the cold component drops significantly from ~24 K in the inner 2 kpc radius to 13 K beyond 6 kpc radial distance, for the best fitting model. The gas-to-dust ratio for β = 1.5, averaged over the galaxy, is higher than the solar value by a factor of 1.5 and is roughly in agreement with the subsolar metallicity of M 33

Robust, data-driven inference in non-linear cosmostatistics

We discuss two projects in non-linear cosmostatistics applicable to very large surveys of galaxies. The first is a Bayesian reconstruction of galaxy redshifts and their number density distribution from approximate, photometric redshift data. The second focuses on cosmic voids and uses them to construct cosmic spheres that allow reconstructing the expansion history of the Universe using the Alcock-Paczynski test. In both cases we find that non-linearities enable the methods or enhance the results: non-linear gravitational evolution creates voids and our photo-z reconstruction works best in the highest density (and hence most non-linear) portions of our simulations.Comment: 14 pages, 10 figures. Talk given at "Statistical Challenges in Modern Astronomy V," held at Penn Stat

Constraints on the small-scale power spectrum of density fluctuations from high-redshift gamma-ray bursts

Cosmological models that include suppression of the power spectrum of density fluctuations on small scales exhibit an exponential reduction of high-redshift, non-linear structures, including a reduction in the rate of gamma ray bursts (GRBs). Here we quantify the constraints that the detection of distant GRBs would place on structure formation models with reduced small-scale power. We compute the number of GRBs that could be detectable by the Swift satellite at high redshifts (z > 6), assuming that the GRBs trace the cosmic star formation history, which itself traces the formation of non-linear structures. We calibrate simple models of the intrinsic luminosity function of the bursts to the number and flux distribution of GRBs observed by the Burst And Transient Source Experiment (BATSE). We find that a discovery of high-z GRBs would imply strong constraints on models with reduced small-scale power. For example, a single GRB at z > 10, or 10 GRBs at z > 5, discovered by Swift during its scheduled two-year mission, would rule out an exponential suppression of the power spectrum on scales below R_c=0.09 Mpc (exemplified by warm dark matter models with a particle mass of m_x=2 keV). Models with a less sharp suppression of small-scale power, such as those with a red tilt or a running scalar index, n_s, are more difficult to constrain, because they are more degenerate with an increase in the power spectrum normalization, sigma_8, and with models in which star-formation is allowed in low-mass minihalos. We find that a tilt of \delta n_s ~ 0.1 is difficult to detect; however, an observed rate of 1 GRB/yr at z > 12 would yield an upper limit on the running of the spectral index, alpha = d(n_s)/d(ln k) > -0.05.Comment: 10 pages, 6 figures; Minor changes to match version published in Ap

Conditional Mass Functions and Merger Rates of Dark Matter Halos in the Ellipsoidal Collapse Model

Analytic models based on spherical and ellipsoidal gravitational collapse have been used to derive the mass functions of dark matter halos and their progenitors (the conditional mass function). The ellipsoidal model generally provides a better match to simulation results, but there has been no simple analytic expression in this model for the conditional mass function that is accurate for small time steps, a limit that is important for generating halo merger trees and computing halo merger rates. We remedy the situation by deriving accurate analytic formulae for the first-crossing distribution, the conditional mass function, and the halo merger rate in the ellipsoidal collapse model in the limit of small look-back times. We show that our formulae provide a closer match to the Millennium simulation results than those in the spherical collapse model and the ellipsoidal model of Sheth & Tormen (2002).Comment: 5 pages, 3 figures, accepted by MNRAS letter

Precision Determination of the Mass Function of Dark Matter Halos

The predicted mass function of dark matter halos is essential in connecting observed galaxy cluster counts and models of galaxy clustering to the properties of the primordial density field. We determine the mass function in the concordance $\Lambda$CDM cosmology, as well as its uncertainty, using sixteen $1024^3$-particle nested-volume dark-matter simulations, spanning a mass range of over five orders of magnitude. Using the nested volumes and single-halo tests, we find and correct for a systematic error in the friends-of-friends halo-finding algorithm. We find a fitting form and full error covariance for the mass function that successfully describes the simulations' mass function and is well-behaved outside the simulations' resolutions. Estimated forecasts of uncertainty in cosmological parameters from future cluster count surveys have negligible contribution from remaining statistical uncertainties in the central cosmology multiplicity function. There exists a potentially non-negligible cosmological dependence (non-universality) of the halo multiplicity function.Comment: 4 pages, 3 figures, submitted to ApJ

Extraskeletal myxoid chondrosarcoma of thigh: a rare case report

Extraskeletal myxoid chondrosarcoma (EMC) are a rare entity of soft tissue tumors that occur predominantly in soft tissue of lower extremities. Here we present a case of 45-year-old female presented with left thigh swelling. MRI finding suggested primary neoplastic lesion. Fine needle aspiration cytology (FNAC) suggested myxoid soft tissue neoplasm/sarcoma. Morphological examination revealed typical extraskeletal myxoid chondrosarcoma with strong immunoreactivity for vimentin and focal immunoreactivity for epithelial membrane antigen (EMA)

Redshift space correlations and scale-dependent stochastic biasing of density peaks

We calculate the redshift space correlation function and the power spectrum of density peaks of a Gaussian random field. In the linear regime k < 0.1 h/Mpc, the redshift space power spectrum is P^s_{pk}(k,u) = exp(-f^2 s_{vel}^2 k^2 u^2) * [b_{pk}(k) + b_{vel}(k) f u^2]^2 * P_m(k), where u is the angle with respect to the line of sight, s_{vel} is the one-dimensional velocity dispersion, f is the growth rate, and b_{pk}(k) and b_{vel}(k) are k-dependent linear spatial and velocity bias factors. For peaks, the value of s_{vel} depends upon the functional form of b_{vel}. The peaks model is remarkable because it has unbiased velocities -- peak motions are driven by dark matter flows -- but, in order to achieve this, b_{vel} is k-dependent. We speculate that this is true in general: k-dependence of the spatial bias will lead to k-dependence of b_{vel} even if the biased tracers flow with the dark matter. Because of the k-dependence of the linear bias parameters, standard manipulations applied to the peak model will lead to k-dependent estimates of the growth factor that could erroneously be interpreted as a signature of modified dark energy or gravity. We use the Fisher formalism to show that the constraint on the growth rate f is degraded by a factor of two if one allows for a k-dependent velocity bias of the peak type. We discuss a simple estimate of nonlinear evolution and illustrate the effect of the peak bias on the redshift space multipoles. For k < 0.1 h/Mpc, the peak bias is deterministic but k-dependent, so the configuration space bias is stochastic and scale dependent, both in real and redshift space. We provide expressions for this stochasticity and its evolution (abridged).Comment: 24 pages, 7 figures (v3): references added (v4): added figure+appendix. In press in PR

Close Pairs as Proxies for Galaxy Cluster Mergers

Galaxy cluster merger statistics are an important component in understanding the formation of large-scale structure. Unfortunately, it is difficult to study merger properties and evolution directly because the identification of cluster mergers in observations is problematic. We use large N-body simulations to study the statistical properties of massive halo mergers, specifically investigating the utility of close halo pairs as proxies for mergers. We examine the relationship between pairs and mergers for a wide range of merger timescales, halo masses, and redshifts (0<z<1). We also quantify the utility of pairs in measuring merger bias. While pairs at very small separations will reliably merge, these constitute a small fraction of the total merger population. Thus, pairs do not provide a reliable direct proxy to the total merger population. We do find an intriguing universality in the relation between close pairs and mergers, which in principle could allow for an estimate of the statistical merger rate from the pair fraction within a scaled separation, but including the effects of redshift space distortions strongly degrades this relation. We find similar behavior for galaxy-mass halos, making our results applicable to field galaxy mergers at high redshift. We investigate how the halo merger rate can be statistically described by the halo mass function via the merger kernel (coagulation), finding an interesting environmental dependence of merging: halos within the mass resolution of our simulations merge less efficiently in overdense environments. Specifically, halo pairs with separations less than a few Mpc/h are more likely to merge in underdense environments; at larger separations, pairs are more likely to merge in overdense environments.Comment: 12 pages, 9 figures; Accepted for publication in ApJ. Significant additions to text and two figures changed. Added new findings on the universality of pair mergers and added analysis of the effect of FoF linking length on halo merger

Searching for the earliest galaxies in the 21 cm forest

We use a model developed by Xu et al. (2010) to compute the 21 cm line absorption signatures imprinted by star-forming dwarf galaxies (DGs) and starless minihalos (MHs). The method, based on a statistical comparison of the equivalent width (W_\nu) distribution and flux correlation function, allows us to derive a simple selection criteria for candidate DGs at very high (z >= 8) redshift. We find that ~ 18% of the total number of DGs along a line of sight to a target radio source (GRB or quasar) can be identified by the condition W_\nu < 0; these objects correspond to the high-mass tail of the DG distribution at high redshift, and are embedded in large HII regions. The criterion W_\nu > 0.37 kHz instead selects ~ 11% of MHs. Selected candidate DGs could later be re-observed in the near-IR by the JWST with high efficiency, thus providing a direct probe of the most likely reionization sources.Comment: 8 pages, 3 figures. Accepted for publication in Science in China Series

The Clustering of Massive Halos

The clustering properties of dark matter halos are a firm prediction of modern theories of structure formation. We use two large volume, high-resolution N-body simulations to study how the correlation function of massive dark matter halos depends upon their mass and formation history. We find that halos with the lowest concentrations are presently more clustered than those of higher concentration, the size of the effect increasing with halo mass; this agrees with trends found in studies of lower mass halos. The clustering dependence on other characterizations of the full mass accretion history appears weaker than the effect with concentration. Using the integrated correlation function, marked correlation functions, and a power-law fit to the correlation function, we find evidence that halos which have recently undergone a major merger or a large mass gain have slightly enhanced clustering relative to a randomly chosen population with the same mass distribution.Comment: 10 pages, 8 figures; text improved, references and one figure added; accepted for publication in Ap