Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “methods: statistical”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 559 records · Page 31
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Correlative piezoresponse and micro-Raman imaging of CuInP 2 S 6 –In 4/3 P 2 S 6 flakes unravels phase-specific phononic fingerprint via unsupervised learning

Characterizing the novel properties of layered van der Waals materials is key for their application in functional devices. A better understanding of this type of material requires correlative imaging of diverse nanoscale material properties. Within this class of materials, CuInP 2 S 6 (CIPS) has received a significant degree of interest due to its ionically mediated room temperature ferroelectricity. Moreover, it is possible to form stable self-assembled heterostructures of ferroelectric CuInP 2 S 6 (CIPS) and non-ferroelectric (i.e., lacking Cu) In 4/3 P 2 S 6 (IPS) phases, by controlling the targeted composition and kinetics of synthesis. In this work, we present a correlative nanometric imaging study of the phononic modes and piezoelectricity of the phase-separated thin heteroepitaxial CIPS/IPS flakes. Here, we show that it is possible to isolate the different phononic modes of the two phases by spatially correlating them with their distinct ferroelectric behavior. The coupling of our experimental data with unsupervised learning statistical methods enables unraveling specific Raman peaks that are characteristic of each chemical phase (CIPS and IPS) present in the composite sample, discarding the less significant ones.

correlative microscopy↗
DOE OSTI2022

Stochastic fluctuation and transport of tokamak edge plasmas with the resonant magnetic perturbation field

Here we present that a statistical method known as the complexity–entropy analysis is useful to characterize a state of plasma turbulence and flux in the resonant magnetic perturbation (RMP) edge localized mode (ELM) control experiment. The stochastic pedestal top temperature fluctuation in the RMP ELM suppression phase is distinguished from the chaotic fluctuation in the natural ELM-free phase. It is discussed that the stochastic temperature fluctuation can be originated from the narrow layer of the field penetration on the pedestal top. The forced magnetic island can emit the resonant drift wave of comparable sizes (relatively low-k) in the RMP ELM suppression phase, and it can result in the generation of stochastic higher wavenumber fluctuations coupled to tangled fields around the island. The analysis of the ion saturation current measurement around the major outer striking point on the divertor shows that it also becomes more stochastic as the stronger plasma response to the RMP field is expected.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Temporal evolution of the light emitted by a thin, laser-ionized plasma source

We present an experimental and simulation-based investigation of the temporal evolution of light emission from a thin, laser-ionized helium plasma source. We demonstrate an analytic model to calculate the approximate scaling of the time-integrated, on-axis light emission with the initial plasma density and temperature, supported by the experiment, which enhances the understanding of plasma light measurement for plasma wakefield accelerator (PWFA) plasma sources. Our model simulates the plasma density and temperature using a split-step Fourier code and a particle-in-cell code. A fluid simulation is then used to model the plasma and neutral density, and the electron temperature as a function of time and position. We then show the numerical results of the space-and-time-resolved light emission and that collisional excitation is the dominant source of light emission. We validate our model by measuring the light emitted by a laser-ionized plasma using a novel statistical method capable of resolving the nanosecond-scale temporal dynamics of the plasma light using a cost-effective camera with microsecond-scale timing jitter. This method is ideal for deployment in the high radiation environment of a particle accelerator that precludes the use of expensive nanosecond-gated cameras. Our results show that our models can effectively simulate the dynamics of a thin, laser-ionized plasma source. In addition, this work provides a detailed understanding of the plasma light measurement, which is one of the few diagnostic signals available for the direct measurement of PWFA plasma sources.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

Bayesian Treed Calibration: An Application to Carbon Capture With AX Sorbent

In cases where field or experimental measurements are not available, computer models can model real physical or engineering systems to reproduce their outcomes. They are usually calibrated in light of experimental data to create a better representation of the real system. Statistical methods, based on Gaussian processes, for calibration and prediction have been especially important when the computer models are expensive and experimental data limited. In this paper, we develop the Bayesian treed calibration (BTC) as an extension of standard Gaussian process calibration methods to deal with non-stationarity computer models and/or their discrepancy from the field (or experimental) data. Our proposed method partitions both the calibration and observable input space, based on a binary tree partitioning, into sub-regions where existing model calibration methods can be applied to connect a computer model with the real system. The estimation of the parameters in the proposed model is carried out using Markov chain Monte Carlo (MCMC) computational techniques. Different strategies have been applied to improve mixing. We illustrate our method in two artificial examples and a real application that concerns the capture of carbon dioxide with AX amine based sorbents. The source code and the examples analyzed in this paper are available as part of the supplementary materials.

Bayesian tree↗
DOE OSTI2011

SILICON AND OXYGEN ABUNDANCES IN PLANET-HOST STARS

The positive correlation between planet detection rate and host star iron abundance lends strong support to the core accretion theory of planet formation. However, iron is not the most significant mass contributor to the cores of giant planets. Since giant planet cores are thought to grow from silicate grains with icy mantles, the likelihood of gas giant formation should depend heavily on the oxygen and silicon abundance of the planet formation environment. Here we compare the silicon and oxygen abundances of a set of 76 planet hosts and a control sample of 80 metal-rich stars without any known giant planets. Our new, independent analysis was conducted using high resolution, high signal-to-noise data obtained at McDonald Observatory. Because we do not wish to simply reproduce the known planet-metallicity correlation, we have devised a statistical method for matching the underlying [Fe/H] distributions of our two sets of stars. We find a 99% probability that planet detection rate depends on the silicon abundance of the host star, over and above the observed planet-metallicity correlation. We do not detect any such correlation for oxygen. Our results would thus seem to suggest that grain nucleation, rather than subsequent icy mantle growth, is the important limiting factor in forming giant planets via core accretion. Based on our results and interpretation, we predict that planet detection should correlate with host star abundance for refractory elements responsible for grain nucleation and that no such trends should exist for the most abundant volatile elements responsible for icy mantle growth.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

COSMOLOGICAL CONSTRAINTS FROM SUNYAEV-ZEL'DOVICH-SELECTED CLUSTERS WITH X-RAY OBSERVATIONS IN THE FIRST 178 deg{sup 2} OF THE SOUTH POLE TELESCOPE SURVEY

We use measurements from the South Pole Telescope (SPT) Sunyaev-Zel'dovich (SZ) cluster survey in combination with X-ray measurements to constrain cosmological parameters. We present a statistical method that fits for the scaling relations of the SZ and X-ray cluster observables with mass while jointly fitting for cosmology. The method is generalizable to multiple cluster observables, and self-consistently accounts for the effects of the cluster selection and uncertainties in cluster mass calibration on the derived cosmological constraints. We apply this method to a data set consisting of an SZ-selected catalog of 18 galaxy clusters at z > 0.3 from the first 178 deg{sup 2} of the 2500 deg{sup 2} SPT-SZ survey, with 14 clusters having X-ray observations from either Chandra or XMM-Newton. Assuming a spatially flat {Lambda}CDM cosmological model, we find the SPT cluster sample constrains {sigma}{sub 8}({Omega} {sub m}/0.25){sup 0.30} = 0.785 {+-} 0.037. In combination with measurements of the cosmic microwave background (CMB) power spectrum from the SPT and the seven-year Wilkinson Microwave Anisotropy Probe data, the SPT cluster sample constrains {sigma}{sub 8} = 0.795 {+-} 0.016 and {Omega} {sub m} = 0.255 {+-} 0.016, a factor of 1.5 improvement on each parameter over the CMB data alone. We consider several extensions beyond the {Lambda}CDM model by including the following as free parameters: the dark energy equation of state (w), the sum of the neutrino masses ({Sigma}m {sub {nu}}), the effective number of relativistic species (N {sub eff}), and a primordial non-Gaussianity (f {sub NL}). We find that adding the SPT cluster data significantly improves the constraints on w and {Sigma}m {sub {nu}} beyond those found when using measurements of the CMB, supernovae, baryon acoustic oscillations, and the Hubble constant. Considering each extension independently, we best constrain w = -0.973 {+-} 0.063 and the sum of neutrino masses {Sigma}m {sub {nu}} < 0.28 eV at 95% confidence, a factor of 1.25 and 1.4 improvement, respectively, over the constraints without clusters. Assuming a {Lambda}CDM model with a free N {sub eff} and {Sigma}m {sub {nu}}, we measure N {sub eff} = 3.91 {+-} 0.42 and constrain {Sigma}m {sub {nu}} < 0.63 eV at 95% confidence. We also use the SPT cluster sample to constrain f {sub NL} = -220 {+-} 317, consistent with zero primordial non-Gaussianity. Finally, we discuss the current systematic limitations due to the cluster mass calibration, and future improvements for the recently completed 2500 deg{sup 2} SPT-SZ survey. The survey has detected {approx}500 clusters with a median redshift of {approx}0.5 and a median mass of {approx}2.3 Multiplication-Sign 10{sup 14} M {sub Sun} h {sup -1} and, when combined with an improved cluster mass calibration and existing external cosmological data sets will significantly improve constraints on w.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

A statistical analysis of seeds and other high-contrast exoplanet surveys: massive planets or low-mass brown dwarfs?

We conduct a statistical analysis of a combined sample of direct imaging data, totalling nearly 250 stars. The stars cover a wide range of ages and spectral types, and include five detections (κ And b, two ∼60 M {sub J} brown dwarf companions in the Pleiades, PZ Tel B, and CD–35 2722B). For some analyses we add a currently unpublished set of SEEDS observations, including the detections GJ 504b and GJ 758B. We conduct a uniform, Bayesian analysis of all stellar ages using both membership in a kinematic moving group and activity/rotation age indicators. We then present a new statistical method for computing the likelihood of a substellar distribution function. By performing most of the integrals analytically, we achieve an enormous speedup over brute-force Monte Carlo. We use this method to place upper limits on the maximum semimajor axis of the distribution function derived from radial-velocity planets, finding model-dependent values of ∼30-100 AU. Finally, we model the entire substellar sample, from massive brown dwarfs to a theoretically motivated cutoff at ∼5 M {sub J}, with a single power-law distribution. We find that p(M, a)∝M {sup –0.65} {sup ±} {sup 0.60} a {sup –0.85} {sup ±} {sup 0.39} (1σ errors) provides an adequate fit to our data, with 1.0%-3.1% (68% confidence) of stars hosting 5-70 M {sub J} companions between 10 and 100 AU. This suggests that many of the directly imaged exoplanets known, including most (if not all) of the low-mass companions in our sample, formed by fragmentation in a cloud or disk, and represent the low-mass tail of the brown dwarfs.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

The formation and disintegration of magnetic bright points observed by sunrise/IMaX

The evolution of the physical parameters of magnetic bright points (MBPs) located in the quiet Sun (mainly in the interwork) during their lifetime is studied. First, we concentrate on the detailed description of the magnetic field evolution of three MBPs. This reveals that individual features follow different, generally complex, and rather dynamic scenarios of evolution. Next, we apply statistical methods on roughly 200 observed MBP evolutionary tracks. MBPs are found to be formed by the strengthening of an equipartition field patch, which initially exhibits a moderate downflow. During the evolution, strong downdrafts with an average velocity of 2.4 km s{sup –1} set in. These flows, taken together with the concurrent strengthening of the field, suggest that we are witnessing the occurrence of convective collapses in these features, although only 30% of them reach kG field strengths. This fraction might turn out to be larger when the new 4 m class solar telescopes are operational as observations of MBPs with current state of the art instrumentation could still be suffering from resolution limitations. Finally, when the bright point disappears (although the magnetic field often continues to exist) the magnetic field strength has dropped to the equipartition level and is generally somewhat weaker than at the beginning of the MBP's evolution. Also, only relatively weak downflows are found on average at this stage of the evolution. Only 16% of the features display upflows at the time that the field weakens, or the MBP disappears. This speaks either for a very fast evolving dynamic process at the end of the lifetime, which could not be temporally resolved, or against strong upflows as the cause of the weakening of the field of these magnetic elements, as has been proposed based on simulation results. It is noteworthy that in about 10% of the cases, we observe in the vicinity of the downflows small-scale strong (exceeding 2 km s{sup –1}) intergranular upflows related spatially and temporally to these downflows. The paper is complemented by a detailed discussion of aspects regarding the applied methods, the complementary literature, and in depth analysis of parameters like magnetic field strength and velocity distributions. An important difference to magnetic elements and associated bright structures in active region plage is that most of the quiet Sun bright points display significant downflows over a large fraction of their lifetime (i.e., in more than 46% of time instances/measurements they show downflows exceeding 1 km s{sup –1}).

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Improved estimates of the Milky Way’s stellar mass and star formation rate from hierarchical Bayesian meta-analysis

We present improved estimates of several global properties of the Milky Way, including its current star formation rate (SFR), the stellar mass contained in its disk and bulge+bar components, as well as its total stellar mass. We do so by combining previous measurements from the literature using a hierarchical Bayesian (HB) statistical method that allows us to account for the possibility that any value may be incorrect or have underestimated errors. We show that this method is robust to a wide variety of assumptions about the nature of problems in individual measurements or error estimates. Ultimately, our analysis yields an SFR for the Galaxy of M-dot {sub ⋆}=1.65±0.19 M{sub ⊙} yr{sup −1}, assuming a Kroupa initial mass function (IMF). By combining HB methods with Monte Carlo simulations that incorporate the latest estimates of the Galactocentric radius of the Sun, R{sub 0}, the exponential scale length of the disk, L{sub d}, and the local surface density of stellar mass, Σ{sub ⋆}(R{sub 0}), we show that the mass of the Galactic bulge+bar is M{sub ⋆}{sup B}=0.91±0.07×10{sup 10} M{sub ⊙}, the disk mass is M{sub ⋆}{sup D}=5.17±1.11×10{sup 10} M{sub ⊙}, and their combination yields a total stellar mass of M{sub ⋆}=6.08±1.14×10{sup 10} M{sub ⊙} (assuming a Kroupa IMF and an exponential disk profile). This analysis is based upon a new compilation of literature bulge mass estimates, normalized to common assumptions about the stellar IMF and Galactic disk properties, presented herein. We additionally find a bulge-to-total mass ratio for the Milky Way of B/T=0.150{sub −0.019}{sup +0.028} and a specific SFR of M-dot {sub ⋆}/M{sub ⋆}=2.71±0.59×10{sup −11} yr{sup −1}.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

DISENTANGLING OVERLAPPING ASTRONOMICAL SOURCES USING SPATIAL AND SPECTRAL INFORMATION

We present a powerful new algorithm that combines both spatial information (event locations and the point-spread function) and spectral information (photon energies) to separate photons from overlapping sources. We use Bayesian statistical methods to simultaneously infer the number of overlapping sources, to probabilistically separate the photons among the sources, and to fit the parameters describing the individual sources. Using the Bayesian joint posterior distribution, we are able to coherently quantify the uncertainties associated with all these parameters. The advantages of combining spatial and spectral information are demonstrated through a simulation study. The utility of the approach is then illustrated by analysis of observations of FK Aqr and FL Aqr with the XMM-Newton Observatory and the central region of the Orion Nebula Cluster with the Chandra X-ray Observatory.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Unveiling the Milky Way: A new technique for determining the optical color and luminosity of our galaxy

We demonstrate a new statistical method of determining the global photometric properties of the Milky Way (MW) to an unprecedented degree of accuracy, allowing our Galaxy to be compared directly to objects measured in extragalactic surveys. Capitalizing on the high-quality imaging and spectroscopy data set from the Sloan Digital Sky Survey (SDSS), we exploit the inherent dependence of galaxies’ luminosities and colors on their total stellar mass, M{sub ⋆}, and star formation rate (SFR), M-dot {sub ⋆}, by selecting a sample of Milky Way analog galaxies designed to reproduce the best Galactic M{sub ⋆} and M-dot {sub ⋆} measurements, including all measurement uncertainties. Making the Copernican assumption that the MW is not extraordinary among galaxies of similar stellar mass and SFR, we then analyze the photometric properties of this matched sample, constraining the characteristics of our Galaxy without suffering interference from interstellar dust. We explore a variety of potential systematic errors that could affect this method, and find that they are subdominant to random uncertainties. We present both SDSS ugriz absolute magnitudes and colors in both rest-frame z = 0 and z = 0.1 passbands for the MW, which are in agreement with previous estimates but can have up to ∼3× lower errors. We find the MW to have absolute magnitude {sup 0}M{sub r}−5logh=−21.00{sub −0.37}{sup +0.38} and integrated color {sup 0}(g−r)=0.682{sub −0.056}{sup +0.066}, indicating that it may belong to the green-valley region in color–magnitude space and ranking it among the brightest and reddest of spiral galaxies. We also present new estimates of global stellar mass-to-light ratios for our Galaxy. This work will help relate our in-depth understanding of the Galaxy to studies of more distant objects.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

The distribution of galaxies’ gravitational field stemming from their tidal interaction

We calculate the distribution function of astronomical objects’ (like galaxies and/or smooth halos of different kinds) gravitational fields due to their tidal interaction. For that we apply the statistical method of Chandrasekhar, used originally to calculate the famous Holtzmark distribution. We show that in our approach the distribution function is never Gaussian, its form being dictated by the potential of interaction between objects. This calculation permits us to perform a theoretical analysis of the relation between angular momentum and mass (richness) of the galaxy clusters. To do so, we follow the ideas of Catelan and Theuns and Heavens and Peacock. The main difference is that here we reduce the problem to a discrete many-body case, where all physical properties of the system are determined by the interaction potential V(r{sub ij}). The essence of reduction is that we use the multipole (up to quadrupole here) expansion of Newtonian potential so that all hydrodynamic, “extended” characteristics of an object, such as its density mass, are “integrated out,” leaving its “point-like” characteristics, such as mass and quadrupole moment. In that sense we do not distinguish between galaxies and smooth components such as halos. We compare our theoretical results with observational data.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Stellar velocity dispersion and anisotropy of the Milky Way inner halo

We measure the three components of velocity dispersion, σ{sub R}, σ{sub θ}, σ{sub ϕ}, for stars within 6 < R < 30 kpc of the Milky Way using a new radial velocity sample from the MMT telescope. We combine our measurements with previously published data so that we can more finely sample the stellar halo. We use a maximum likelihood statistical method for estimating mean velocities, dispersions, and covariances assuming only that velocities are normally distributed. The alignment of the velocity ellipsoid is consistent with a spherically symmetric gravitational potential. From the spherical Jeans equation, the mass of the Milky Way is M(R⩽12 kpc)=1.3×10{sup 11} M{sub ⊙} with an uncertainty of 40%. We also find a region of discontinuity, 15 ≲ R ≲ 25 kpc, where the estimated velocity dispersions and anisotropies diverge from their anticipated values, confirming the break observed by others. We argue that this break in anisotropy is physically explained by coherent stellar velocity structure in the halo, such as the Sgr stream. To significantly improve our understanding of halo kinematics will require combining radial velocities with future Gaia proper motions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

Blob-hole correlation model for edge turbulence and comparisons with NSTX gas puff imaging data [Blob-hole correlation model for edge turbulence and comparisons with NSTX GPI data]

We report that gas puff imaging (GPI) observations made in NSTX have revealed two-point spatial correlations of edge and scrape-off layer turbulence in the plane perpendicular to the magnetic field. A common feature is the occurrence of dipole-like patterns with significant regions of negative correlation. In this paper, we explore the possibility that these dipole patterns may be due to blob-hole pairs. Statistical methods are applied to determine the two-point spatial correlation that results from a model of blob-hole pair formation. It is shown that the model produces dipole correlation patterns that are qualitatively similar to the GPI data in several respects. Effects of the reference location (confined surfaces or scrape-off layer), a superimposed random background, hole velocity and lifetime, and background sheared flows are explored and discussed with respect to experimental observations. Additional analysis of the experimental GPI dataset is performed to further test this blob-hole correlation model. A time delay two-point spatial correlation study did not reveal inward propagation of the negative correlation structures that were postulated to correspond to holes in the data nor did it suggest that the negative correlation structures are due to neutral shadowing. However, tracking of the highest and lowest values (extrema) of the normalized GPI fluctuations shows strong evidence for mean inward propagation of minima and outward propagation of maxima, in qualitative agreement with theoretical expectations. Finally, other properties of the experimentally observed extrema are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2011

A Bayesian view of the current status of dark matter direct searches

Bayesian statistical methods offer a simple and consistent framework for incorporating uncertainties into a multi-parameter inference problem. In this work we apply these methods to a selection of current direct dark matter searches. We consider the simplest scenario of spin-independent elastic WIMP scattering, and infer the WIMP mass and cross-section from the experimental data with the essential systematic uncertainties folded into the analysis. We find that when uncertainties in the scintillation efficiency of XENON100 have been accounted for, the resulting exclusion limit is not sufficiently constraining to rule out the CoGeNT preferred parameter region, contrary to previous claims. In the same vein, we also investigate the impact of astrophysical uncertainties on the preferred WIMP parameters. We find that within the class of smooth and isotropic WIMP velocity distributions, it is difficult to reconcile the DAMA and the CoGeNT preferred regions by tweaking the astrophysics parameters alone. If we demand compatibility between these experiments, then the inference process naturally concludes that a high value for the sodium quenching factor for DAMA is preferred.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2012

Fermi LAT search for internal bremsstrahlung signatures from dark matter annihilation

A commonly encountered obstacle in indirect searches for galactic dark matter is how to disentangle possible signals from astrophysical backgrounds. Given that such signals are most likely subdominant, the search for pronounced spectral features plays a key role for indirect detection experiments; monochromatic gamma-ray lines or similar features related to internal bremsstrahlung, in particular, provide smoking gun signatures. We perform a dedicated search for the latter in the data taken by the Fermi gamma-ray space telescope during its first 43 months. To this end, we use a new adaptive procedure to select optimal target regions that takes into account both standard and contracted dark matter profiles. The behaviour of our statistical method is tested by a subsampling analysis of the full sky data and found to reproduce the theoretical expectations very well. The limits on the dark matter annihilation cross-section that we derive are stronger than what can be obtained from the observation of dwarf galaxies and, at least for the model considered here, collider searches. While these limits are still not quite strong enough to probe annihilation rates expected for thermally produced dark matter, future prospects to do so are very good. In fact, we already find a weak indication, with a significance of 3.1σ (4.3σ) when (not) taking into account the look-elsewhere effect, for an internal bremsstrahlung-like signal that would correspond to a dark matter mass of ∼150 GeV; the same signal is also well fitted by a gamma-ray line at around 130 GeV. Although this would be a fascinating possibility, we caution that a much more dedicated analysis and additional data will be necessary to rule out or confirm this option.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Analysis of the theoretical bias in dark matter direct detection

Fitting the model ''A'' to dark matter direct detection data, when the model that underlies the data is ''B'', introduces a theoretical bias in the fit. We perform a quantitative study of the theoretical bias in dark matter direct detection, with a focus on assumptions regarding the dark matter interactions, and velocity distribution. We address this problem within the effective theory of isoscalar dark matter-nucleon interactions mediated by a heavy spin-1 or spin-0 particle. We analyze 24 benchmark points in the parameter space of the theory, using frequentist and Bayesian statistical methods. First, we simulate the data of future direct detection experiments assuming a momentum/velocity dependent dark matter-nucleon interaction, and an anisotropic dark matter velocity distribution. Then, we fit a constant scattering cross section, and an isotropic Maxwell-Boltzmann velocity distribution to the simulated data, thereby introducing a bias in the analysis. The best fit values of the dark matter particle mass differ from their benchmark values up to 2 standard deviations. The best fit values of the dark matter-nucleon coupling constant differ from their benchmark values up to several standard deviations. We conclude that common assumptions in dark matter direct detection are a source of potentially significant bias.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Angular correlation of cosmic neutrinos with ultrahigh-energy cosmic rays and implications for their sources

Cosmic neutrino events detected by the IceCube Neutrino Observatory with energy 0∼> 3 TeV have poor angular resolutions to reveal their origin. Ultrahigh-energy cosmic rays (UHECRs), with better angular resolutions at 0>6 EeV energies, can be used to check if the same astrophysical sources are responsible for producing both neutrinos and UHECRs. We test this hypothesis, with statistical methods which emphasize invariant quantities, by using data from the Pierre Auger Observatory, Telescope Array and past cosmic-ray experiments. We find that the arrival directions of the cosmic neutrinos are correlated with 0≥ 10 EeV UHECR arrival directions at confidence level ≈ 90%. The strength of the correlation decreases with decreasing UHECR energy and no correlation exists at energy 0∼ 6 EeV . A search in astrophysical databases within 3{sup o} of the arrival directions of UHECRs with energy 0≥ 10 EeV, that are correlated with the IceCube cosmic neutrinos, resulted in 18 sources from the Swift-BAT X-ray catalog with redshift z≤ 0.06. We also found 3 objects in the Kühr catalog of radio sources using the same criteria. The sources are dominantly Seyfert galaxies with Cygnus A being the most prominent member. We calculate the required neutrino and UHECR fluxes to produce the observed correlated events, and estimate the corresponding neutrino luminosity (25 TeV–2.2 PeV) and cosmic-ray luminosity (500 TeV–180 EeV), assuming the sources are the ones we found in the Swift-BAT and Kühr catalogs. We compare these luminosities with the X-ray luminosity of the corresponding sources and discuss possibilities of accelerating protons to 0∼> 10 EeV and produce neutrinos in these sources.

79 ASTRONOMY AND ASTROPHYSICS↗