Astrophysical Uncertainties in the Cosmic Ray Electron and Positron Spectrum From Annihilating Dark Matter

In recent years, a number of experiments have been conducted with the goal of studying cosmic rays at GeV to TeV energies. This is a particularly interesting regime from the perspective of indirect dark matter detection. To draw reliable conclusions regarding dark matter from cosmic ray measurements, however, it is important to first understand the propagation of cosmic rays through the magnetic and radiation fields of the Milky Way. In this paper, we constrain the characteristics of the cosmic ray propagation model through comparison with observational inputs, including recent data from the CREAM experiment, and use these constraints to estimate the corresponding uncertainties in the spectrum of cosmic ray electrons and positrons from dark matter particles annihilating in the halo of the Milky Way.Comment: 21 pages, 9 figure

A New Approach to Searching for Dark Matter Signals in Fermi-LAT Gamma Rays

Several cosmic ray experiments have measured excesses in electrons and positrons, relative to standard backgrounds, for energies from ~ 10 GeV - 1 TeV. These excesses could be due to new astrophysical sources, but an explanation in which the electrons and positrons are dark matter annihilation or decay products is also consistent. Fortunately, the Fermi-LAT diffuse gamma ray measurements can further test these models, since the electrons and positrons produce gamma rays in their interactions in the interstellar medium. Although the dark matter gamma ray signal consistent with the local electron and positron measurements should be quite large, as we review, there are substantial uncertainties in the modeling of diffuse backgrounds and, additionally, experimental uncertainties that make it difficult to claim a dark matter discovery. In this paper, we introduce an alternative method for understanding the diffuse gamma ray spectrum in which we take the intensity ratio in each energy bin of two different regions of the sky, thereby canceling common systematic uncertainties. For many spectra, this ratio fits well to a power law with a single break in energy. The two measured exponent indices are a robust discriminant between candidate models, and we demonstrate that dark matter annihilation scenarios can predict index values that require "extreme" parameters for background-only explanations.Comment: v1: 11 pages, 7 figures, 1 table, revtex4; v2: 13 pages, 8 figures, 1 table, revtex4, Figure 4 added, minor additions made to text, references added, conclusions unchanged, published versio

Event Reconstruction in the PHENIX Central Arm Spectrometers

The central arm spectrometers for the PHENIX experiment at the Relativistic Heavy Ion Collider have been designed for the optimization of particle identification in relativistic heavy ion collisions. The spectrometers present a challenging environment for event reconstruction due to a very high track multiplicity in a complicated, focusing, magnetic field. In order to meet this challenge, nine distinct detector types are integrated for charged particle tracking, momentum reconstruction, and particle identification. The techniques which have been developed for the task of event reconstruction are described.Comment: Accepted for publication in Nucl. Instrum. A. 34 pages, 23 figure

Measurements of double-helicity asymmetries in inclusive J/ψJ/\psi production in longitudinally polarized p+pp+p collisions at s=510\sqrt{s}=510 GeV

We report the double helicity asymmetry, $A_{LL}^{J/\psi}$, in inclusive $J/\psi$ production at forward rapidity as a function of transverse momentum $p_T$ and rapidity $|y|$. The data analyzed were taken during $\sqrt{s}=510$ GeV longitudinally polarized $p$$+$$p$ collisions at the Relativistic Heavy Ion Collider (RHIC) in the 2013 run using the PHENIX detector. At this collision energy, $J/\psi$ particles are predominantly produced through gluon-gluon scatterings, thus $A_{LL}^{J/\psi}$ is sensitive to the gluon polarization inside the proton. We measured $A_{LL}^{J/\psi}$ by detecting the decay daughter muon pairs $\mu^+ \mu^-$ within the PHENIX muon spectrometers in the rapidity range $1.2<|y|<2.2$. In this kinematic range, we measured the $A_{LL}^{J/\psi}$ to be $0.012 \pm 0.010$~(stat)~$\pm$~$0.003$(syst). The $A_{LL}^{J/\psi}$ can be expressed to be proportional to the product of the gluon polarization distributions at two distinct ranges of Bjorken $x$: one at moderate range $x \approx 0.05$ where recent RHIC data of jet and $\pi^0$ double helicity spin asymmetries have shown evidence for significant gluon polarization, and the other one covering the poorly known small-$x$ region $x \approx 2\times 10^{-3}$. Thus our new results could be used to further constrain the gluon polarization for $x< 0.05$.Comment: 335 authors, 10 pages, 4 figures, 3 tables, 2013 data. Version accepted for publication by Phys. Rev. D. Plain text data tables for the points plotted in figures for this and previous PHENIX publications are (or will be) publicly available at http://www.phenix.bnl.gov/papers.htm

Heterozygosity-fitness correlations in a wild mammal population: accounting for parental and environmental effects

HFCs (heterozygosity–fitness correlations) measure the direct relationship between an individual's genetic diversity and fitness. The effects of parental heterozygosity and the environment on HFCs are currently under-researched. We investigated these in a high-density U.K. population of European badgers (Meles meles), using a multimodel capture–mark–recapture framework and 35 microsatellite loci. We detected interannual variation in first-year, but not adult, survival probability. Adult females had higher annual survival probabilities than adult males. Cubs with more heterozygous fathers had higher first-year survival, but only in wetter summers; there was no relationship with individual or maternal heterozygosity. Moist soil conditions enhance badger food supply (earthworms), improving survival. In dryer years, higher indiscriminate mortality rates appear to mask differential heterozygosity-related survival effects. This paternal interaction was significant in the most supported model; however, the model-averaged estimate had a relative importance of 0.50 and overlapped zero slightly. First-year survival probabilities were not correlated with the inbreeding coefficient (f); however, small sample sizes limited the power to detect inbreeding depression. Correlations between individual heterozygosity and inbreeding were weak, in line with published meta-analyses showing that HFCs tend to be weak. We found support for general rather than local heterozygosity effects on first-year survival probability, and g2 indicated that our markers had power to detect inbreeding. We emphasize the importance of assessing how environmental stressors can influence the magnitude and direction of HFCs and of considering how parental genetic diversity can affect fitness-related traits, which could play an important role in the evolution of mate choice

Centrality categorization for R_{p(d)+A} in high-energy collisions

High-energy proton- and deuteron-nucleus collisions provide an excellent tool for studying a wide array of physics effects, including modifications of parton distribution functions in nuclei, gluon saturation, and color neutralization and hadronization in a nuclear environment, among others. All of these effects are expected to have a significant dependence on the size of the nuclear target and the impact parameter of the collision, also known as the collision centrality. In this article, we detail a method for determining centrality classes in p(d)+A collisions via cuts on the multiplicity at backward rapidity (i.e., the nucleus-going direction) and for determining systematic uncertainties in this procedure. For d+Au collisions at sqrt(s_NN) = 200 GeV we find that the connection to geometry is confirmed by measuring the fraction of events in which a neutron from the deuteron does not interact with the nucleus. As an application, we consider the nuclear modification factors R_{p(d)+A}, for which there is a potential bias in the measured centrality dependent yields due to auto-correlations between the process of interest and the backward rapidity multiplicity. We determine the bias correction factor within this framework. This method is further tested using the HIJING Monte Carlo generator. We find that for d+Au collisions at sqrt(s_NN)=200 GeV, these bias corrections are small and vary by less than 5% (10%) up to p_T = 10 (20) GeV. In contrast, for p+Pb collisions at sqrt(s_NN) = 5.02 TeV we find these bias factors are an order of magnitude larger and strongly p_T dependent, likely due to the larger effect of multi-parton interactions.Comment: 375 authors, 18 pages, 16 figures, 4 tables. Submitted to Phys. Rev. C. Plain text data tables for the points plotted in figures for this and previous PHENIX publications are (or will be) publicly available at http://www.phenix.bnl.gov/papers.htm

L\'evy-stable two-pion Bose-Einstein correlations in sNN=200\sqrt{s_{_{NN}}}=200 GeV Au++Au collisions

We present a detailed measurement of charged two-pion correlation functions in 0%-30% centrality $\sqrt{s_{_{NN}}}=200$ GeV Au$+$Au collisions by the PHENIX experiment at the Relativistic Heavy Ion Collider. The data are well described by Bose-Einstein correlation functions stemming from L\'evy-stable source distributions. Using a fine transverse momentum binning, we extract the correlation strength parameter $\lambda$, the L\'evy index of stability $\alpha$ and the L\'evy length scale parameter $R$ as a function of average transverse mass of the pair $m_T$. We find that the positively and the negatively charged pion pairs yield consistent results, and their correlation functions are represented, within uncertainties, by the same L\'evy-stable source functions. The $\lambda(m_T)$ measurements indicate a decrease of the strength of the correlations at low $m_T$. The L\'evy length scale parameter $R(m_T)$ decreases with increasing $m_T$, following a hydrodynamically predicted type of scaling behavior. The values of the L\'evy index of stability $\alpha$ are found to be significantly lower than the Gaussian case of $\alpha=2$, but also significantly larger than the conjectured value that may characterize the critical point of a second-order quark-hadron phase transition.Comment: 448 authors, 25 pages, 11 figures, 4 tables, 2010 data. v2 is version accepted for publication in Phys. Rev. C. Plain text data tables for the points plotted in figures for this and previous PHENIX publications are (or will be) publicly available at http://www.phenix.bnl.gov/papers.htm

Upsilon (1S+2S+3S) production in d+Au and p+p collisions at sqrt(s_NN)=200 GeV and cold-nuclear matter effects

The three Upsilon states, Upsilon(1S+2S+3S), are measured in d+Au and p+p collisions at sqrt(s_NN)=200 GeV and rapidities 1.2<|y|<2.2 by the PHENIX experiment at the Relativistic Heavy-Ion Collider. Cross sections for the inclusive Upsilon(1S+2S+3S) production are obtained. The inclusive yields per binary collision for d+Au collisions relative to those in p+p collisions (R_dAu) are found to be 0.62 +/- 0.26 (stat) +/- 0.13 (syst) in the gold-going direction and 0.91 +/- 0.33 (stat) +/- 0.16 (syst) in the deuteron-going direction. The measured results are compared to a nuclear-shadowing model, EPS09 [JHEP 04, 065 (2009)], combined with a final-state breakup cross section, sigma_br, and compared to lower energy p+A results. We also compare the results to the PHENIX J/psi results [Phys. Rev. Lett. 107, 142301 (2011)]. The rapidity dependence of the observed Upsilon suppression is consistent with lower energy p+A measurements.Comment: 495 authors, 11 pages, 9 figures, 5 tables. Submitted to Phys. Rev. C. Plain text data tables for the points plotted in figures for this and previous PHENIX publications are (or will be) publicly available at http://www.phenix.bnl.gov/papers.htm

Quadrupole Anisotropy in Dihadron Azimuthal Correlations in Central dd++Au Collisions at sNN\sqrt{s_{_{NN}}}=200 GeV

The PHENIX collaboration at the Relativistic Heavy Ion Collider (RHIC) reports measurements of azimuthal dihadron correlations near midrapidity in $d$$+$Au collisions at $\sqrt{s_{_{NN}}}$=200 GeV. These measurements complement recent analyses by experiments at the Large Hadron Collider (LHC) involving central $p$$+$Pb collisions at $\sqrt{s_{_{NN}}}$=5.02 TeV, which have indicated strong anisotropic long-range correlations in angular distributions of hadron pairs. The origin of these anisotropies is currently unknown. Various competing explanations include parton saturation and hydrodynamic flow. We observe qualitatively similar, but larger, anisotropies in $d$$+$Au collisions compared to those seen in $p$$+$Pb collisions at the LHC. The larger extracted $v_2$ values in $d$$+$Au collisions at RHIC are consistent with expectations from hydrodynamic calculations owing to the larger expected initial-state eccentricity compared with that from $p$$+$Pb collisions. When both are divided by an estimate of the initial-state eccentricity the scaled anisotropies follow a common trend with multiplicity that may extend to heavy ion data at RHIC and the LHC, where the anisotropies are widely thought to arise from hydrodynamic flow.Comment: 375 authors, 7 pages, 5 figures. Published in Phys. Rev. Lett. v2 has minor changes to text and figures in response to PRL referee suggestions. Plain text data tables for the points plotted in figures for this and previous PHENIX publications are (or will be) publicly available at http://www.phenix.bnl.gov/papers.htm