Hierarchical structure formation implies that the number of subhalos within a dark matter halo depends not only on halo mass, but also on the formation history of the halo. This dependence on the formation history, which is highly correlated with hal
o concentration, can account for the super-Poissonian scatter in subhalo occupation at a fixed halo mass that has been previously measured in simulations. Here we propose a model to predict the subhalo abundance function for individual host halos, that incorporates both halo mass and concentration. We combine results of cosmological simulations with a new suite of zoom-in simulations of Milky Way-mass halos to calibrate our model. We show the model can successfully reproduce the mean and the scatter of subhalo occupation in these simulations. The implications of this correlation between subhalo abundance and halo concentration are further investigated. We also discuss cases in which inferences about halo properties can be affected if this correlation between subhalo abundance and halo concentration is ignored; in these cases our model would give a more accurate inference. We propose that with future deep surveys, satellite occupation in the low-mass regime can be used to verify the existence of halo assembly bias.
A dramatic increase in the accuracy and statistics of space-borne cosmic ray (CR) measurements has yielded several breakthroughs over the last several years. The most puzzling is the rise in the positron fraction above ~10 GeV over the predictions of
the propagation models assuming pure secondary production. The accuracy of the antiproton production cross section is critical for astrophysical applications and searches for new physics since antiprotons in CRs seem to hold the keys to many puzzles including the origin of those excess positrons. However, model calculations of antiproton production in CR interactions with interstellar gas are often employing parameterizations that are out of date or are using outdated physical concepts. That may lead to an incorrect interpretation of antiproton data which could have broad consequences for other areas of astrophysics. In this work, we calculate antiproton production in pp-, pA-, and AA-interactions using EPOS-LHC and QGSJET-II-04, two of the most advanced Monte Carlo (MC) generators tuned to numerous accelerator data including those from the Large Hadron Collider (LHC). We show that the antiproton yields obtained with these MC generators differ by up to an order of magnitude from yields of parameterizations commonly used in astrophysics.
We develop a general framework for data analysis and phenomenology of the CMB four-point function or trispectrum. To lowest order in the derivative expansion, the inflationary action admits three quartic operators consistent with symmetry: $dotsigma^
4$, $dotsigma^2 (partialsigma^2)$, and $(partialsigma)^4$. In single field inflation, only the first of these operators can be the leading non-Gaussian signal. A Fisher matrix analysis shows that there is one near-degeneracy among the three CMB trispectra, so we parameterize the trispectrum with two coefficients $g_{NL}^{dotsigma^4}$ and $g_{NL}^{(partialsigma)^4}$, in addition to the coefficient $g_{NL}^{rm loc}$ of $zeta^3$-type local non-Gaussianity. This three-parameter space is analogous to the parameter space $(f_{NL}^{rm loc}, f_{NL}^{rm equil}, f_{NL}^{rm orth})$ commonly used to parameterize the CMB three-point function. We next turn to data analysis and show how to represent these trispectra in a factorizable form which leads to computationally fast operations such as evaluating a CMB estimator or simulating a non-Gaussian CMB. We discuss practical issues in CMB analysis pipelines, and perform an optimal analysis of WMAP data. Our minimum-variance estimates are $g_{NL}^{rm loc} = (-3.80 pm 2.19) times 10^5$, $g_{NL}^{dotsigma^4} = (-3.20 pm 3.09) times 10^6$, and $g_{NL}^{(partialsigma)^4} = (-10.8 pm 6.33) times 10^5$ after correcting for the effects of CMB lensing. No evidence of a nonzero inflationary four-point function is seen.
X-ray surface brightness fluctuations in the core of the Perseus Cluster are analyzed, using deep observations with the Chandra observatory. The amplitude of gas density fluctuations on different scales is measured in a set of radial annuli. It varie
s from 8 to 12 per cent on scales of ~10-30 kpc within radii of 30-160 kpc from the cluster center and from 9 to 7 per cent on scales of ~20-30 kpc in an outer, 60-220 kpc annulus. Using a statistical linear relation between the observed amplitude of density fluctuations and predicted velocity, the characteristic velocity of gas motions on each scale is calculated. The typical amplitudes of the velocity outside the central 30 kpc region are 90-140 km/s on ~20-30 kpc scales and 70-100 km/s on smaller scales ~7-10 kpc. The velocity power spectrum is consistent with cascade of turbulence and its slope is in a broad agreement with the slope for canonical Kolmogorov turbulence. The gas clumping factor estimated from the power spectrum of the density fluctuations is lower than 7-8 per cent for radii ~30-220 kpc from the center, leading to a density bias of less than 3-4 per cent in the cluster core. Uncertainties of the analysis are examined and discussed. Future measurements of the gas velocities with the Astro-H, Athena and Smart-X observatories will directly measure the gas density-velocity perturbation relation and further reduce systematic uncertainties in these quantities.
The hot, X-ray-emitting intracluster medium (ICM) is the dominant baryonic constituent of clusters of galaxies. In the cores of many clusters, radiative energy losses from the ICM occur on timescales significantly shorter than the age of the system.
Unchecked, this cooling would lead to massive accumulations of cold gas and vigorous star formation, in contradiction to observations. Various sources of energy capable of compensating these cooling losses have been proposed, the most promising being heating by the supermassive black holes in the central galaxies through inflation of bubbles of relativistic plasma. Regardless of the original source of energy, the question of how this energy is transferred to the ICM has remained open. Here we present a plausible solution to this question based on deep Chandra X-ray observatory data and a new data-analysis method that enables us to evaluate directly the ICM heating rate due to the dissipation of turbulence. We find that turbulent heating is sufficient to offset radiative cooling and indeed appears to balance it locally at each radius - it might therefore be the key element in resolving the gas cooling problem in cluster cores and, more universally, in atmospheres of X-ray gas-rich systems.
In order to study the galaxy population of galaxy clusters with photometric data one must be able to accurately discriminate between cluster members and non-members. The redMaPPer cluster finding algorithm treats this problem probabilistically. Here,
we utilize SDSS and GAMA spectroscopic membership rates to validate the redMaPPer membership probability estimates for clusters with $zin[0.1,0.3]$. We find small - but correctable - biases, sourced by three different systematics. The first two were expected a priori, namely blue cluster galaxies and correlated structure along the line of sight. The third systematic is new: the redMaPPer template fitting exhibits a non-trivial dependence on photometric noise, which biases the original redMaPPer probabilities when utilizing noisy data. After correcting for these effects, we find exquisite agreement ($approx 1%$) between the photometric probability estimates and the spectroscopic membership rates, demonstrating that we can robustly recover cluster membership estimates from photometric data alone. As a byproduct of our analysis we find that on average unavoidable projection effects from correlated structure contribute $approx 6%$ of the richness of a redMaPPer galaxy cluster. This work also marks the second public release of the SDSS redMaPPer cluster catalog.
We have obtained a deep, simultaneous observation of the bright, nearby Seyfert galaxy IC 4329A with Suzaku and NuSTAR. Through a detailed spectral analysis, we are able to robustly separate the continuum, absorption and distant reflection components
in the spectrum. The absorbing column is found to be modest at $N_H = 6 times 10^{21}$ cm$^2$, and does not introduce any significant curvature in the Fe K band. We are able to place a strong constraint on the presence of a broadened Fe K{alpha} line: $E = 6.46^{+0.08}_{-0.07}$ keV rest frame with ${sigma} = 0.33^{+0.08}_{-0.07}$ keV and $EW = 34^{+8}_{-7}$ eV, though we are not able to constrain any of the parameters of a relativistic reflection model. These results highlight the range in broad Fe K{alpha} line strengths observed in nearby, bright AGN (roughly an order of magnitude), and imply a corresponding range in the physical properties of the inner accretion disk in these sources. We have also updated our previously reported measurement of the high-energy cutoff of the hard X-ray emission using both observatories rather than just NuSTAR alone: $E_{cut} = 186 pm 14$ keV. This high-energy cutoff acts as a proxy for the temperature of the coronal electron plasma, enabling us to further separate this parameter from the optical depth of the plasma and to update our results for these parameters as well. We derive $kT = 50^{+6}_{-3}$ keV with ${tau} = 2.34^{+0.16}_{-0.11}$ using a spherical geometry, $kT = 61 pm 1$ keV with ${tau} = 0.68 pm 0.02$ for a slab geometry, with both having an equivalent goodness-of-fit.
We address the problem of evaluating the power spectrum of the velocity field of the ICM using only information on the plasma density fluctuations, which can be measured today by Chandra and XMM-Newton observatories. We argue that for relaxed cluster
s there is a linear relation between the rms density and velocity fluctuations across a range of scales, from the largest ones, where motions are dominated by buoyancy, down to small, turbulent scales: $(deltarho_k/rho)^2 = eta_1^2 (V_{1,k}/c_s)^2$, where $deltarho_k/rho$ is the spectral amplitude of the density perturbations at wave number $k$, $V_{1,k}^2=V_k^2/3$ is the mean square component of the velocity field, $c_s$ is the sound speed, and $eta_1$ is a dimensionless constant of order unity. Using cosmological simulations of relaxed galaxy clusters, we calibrate this relation and find $eta_1approx 1 pm 0.3$. We argue that this value is set at large scales by buoyancy physics, while at small scales the density and velocity power spectra are proportional because the former are a passive scalar advected by the latter. This opens an interesting possibility to use gas density power spectra as a proxy for the velocity power spectra in relaxed clusters, across a wide range of scales.
We present an analysis of a ~160 ks NuSTAR observation of the nearby bright Seyfert galaxy IC4329A. The high-quality broadband spectrum enables us to separate the effects of distant reflection from the direct coronal continuum, and to therefore accur
ately measure the high-energy cutoff to be $E_{cut}=178^{+74}_{-40}$ keV. The coronal emission arises from accretion disk photons Compton up-scattered by a thermal plasma, with the spectral index and cutoff being due to a combination of the finite plasma temperature and optical depth. Applying standard Comptonization models, we measure both physical properties independently using the best signal-to-noise obtained to date in an AGN over the 3-79 keV band. We derive $kT_e=37^{+7}_{-6}$ keV with $tau=1.25^{+0.20}_{-0.10}$ assuming a slab geometry for the plasma, and $kT_e=33^{+6}_{-6}$ keV with $tau=3.41^{+0.58}_{-0.38}$ for a spherical geometry, with both having an equivalent goodness-of-fit.
We analyze the star formation properties of 16 infrared-selected, spectroscopically confirmed galaxy clusters at $1 < z < 1.5$ from the Spitzer/IRAC Shallow Cluster Survey (ISCS). We present new spectroscopic confirmation for six of these high-redshi
ft clusters, five of which are at $z>1.35$. Using infrared luminosities measured with deep Spitzer/MIPS observations at 24 $mu$m, along with robust optical+IRAC photometric redshifts and SED-fitted stellar masses, we present the dust-obscured star-forming fractions, star formation rates and specific star formation rates in these clusters as functions of redshift and projected clustercentric radius. We find that $zsim 1.4$ represents a transition redshift for the ISCS sample, with clear evidence of an unquenched era of cluster star formation at earlier times. Beyond this redshift the fraction of star-forming cluster members increases monotonically toward the cluster centers. Indeed, the specific star formation rate in the cores of these distant clusters is consistent with field values at similar redshifts, indicating that at $z>1.4$ environment-dependent quenching had not yet been established in ISCS clusters. Combining these observations with complementary studies showing a rapid increase in the AGN fraction, a stochastic star formation history, and a major merging episode at the same epoch in this cluster sample, we suggest that the starburst activity is likely merger-driven and that the subsequent quenching is due to feedback from merger-fueled AGN. The totality of the evidence suggests we are witnessing the final quenching period that brings an end to the era of star formation in galaxy clusters and initiates the era of passive evolution.
