Source author record

Ke Fang

Ke Fang appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

28works
9topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

28 published item(s)

preprint2026arXiv

MOSAIC: Module Discovery via Sparse Additive Identifiable Causal Learning for Scientific Time Series

Causal representation learning (CRL) seeks to recover latent variables with identifiability guarantees, typically up to permutation and component-wise reparameterization under appropriate assumptions. However, identifiability does not imply interpretability: latent semantics are typically assigned post hoc by alignment with known ground-truth factors. This limitation is particularly acute in scientific time series, where underlying mechanisms are unknown and discovering interpretable structure is a primary goal. In contrast, scientific observations (such as residue-pair distances, climate indices, or process sensors) are inherently semantic, as they correspond to named physical quantities. This raises a key question: can the interpretability of observations be transferred to the identifiable latent space? We propose MOSAIC (Module discovery via Sparse Additive Identifiable Causal learning), a sparse temporal VAE that integrates temporal CRL identifiability with support recovery over observed variables. MOSAIC identifies latent variables via regime-conditioned temporal variation, and recovers for each latent a sparse set of associated observations through an additive decoder, yielding module-level interpretability. We show that ANOVA main-effect supports are identifiable under general smooth mixing functions, and provide finite-sample recovery guarantees for a tractable sparse-additive variant. Empirically, MOSAIC recovers domain-consistent variable groups across RNA molecular dynamics, solar wind, ENSO climate, the Tennessee Eastman process, and a synthetic tokamak benchmark, enabling interpretable discovery of latent mechanisms in scientific time series.

preprint2024arXiv

High-Energy Neutrinos from Gamma-Ray-Faint Accretion-Powered Hypernebulae

Hypernebulae are inflated by accretion-powered winds accompanying hyper-Eddington mass transfer from an evolved post-main sequence star onto a black hole or neutron star companion. The ions accelerated at the termination shock -- where the collimated fast disk winds/jet collide with the slower, wide-angled wind-fed shell -- can generate high-energy neutrinos via hadronic ($pp$) reactions, and photohadronic ($pγ$) interactions with the disk thermal and Comptonized nonthermal background photons. It has been suggested that some fast radio bursts (FRBs) may be powered by such short-lived jetted hyper-accreting engines. Although neutrino emission associated with the ms-duration bursts themselves is challenging to detect, the persistent radio counterparts of some FRB sources -- if associated with hypernebulae -- could contribute to the high energy neutrino diffuse background flux. If the hypernebula birth rate follows that of steller-merger transients and common envelope events, we find that their volume-integrated neutrino emission -- depending on the population-averaged mass-transfer rates -- could explain up to $\sim25\%$ of the high-energy diffuse neutrino flux observed by the IceCube Observatory and the Baikal-Gigaton Volume Detector (GVD) Telescope. The time-averaged neutrino spectrum from hypernebula -- depending on the population parameters -- can also reproduce the observed diffuse neutrino spectrum. The neutrino emission could in some cases furthermore extend to >100 PeV, detectable by future ultra-high-energy neutrino observatories. The large optical depth through the nebula to Breit-Wheeler ($γγ$) interaction attenuates the escape of GeV-PeV gamma-rays co-produced with the neutrinos, rendering these gamma-ray-faint neutrino sources, consistent with the \textit{Fermi} observations of the isotropic gamma-ray background.

preprint2023arXiv

Neutrinos from the Brightest Gamma-Ray Burst?

We discuss implications that can be obtained by searches for neutrinos from the brightest gamma-ray burst, GRB 221009A. We derive constraints on GRB model parameters such as the cosmic-ray loading factor and dissipation radius, taking into account both neutrino spectra and effective areas. The results are strong enough to constrain proton acceleration near the photosphere, and we find that the single burst limits are comparable to those from stacking analysis. Quasithermal neutrinos from subphotospheres and ultrahigh-energy neutrinos from external shocks are not yet constrained. We show that GeV-TeV neutrinos originating from neutron collisions are detectable, and urge dedicated analysis on these neutrinos with DeepCore and IceCube as well as ORCA and KM3NeT.

preprint2022arXiv

Evidence for PeV Proton Acceleration from Fermi-LAT Observations of SNR G106.3+2.7

The existence of a "knee" at energy ~1 PeV in the cosmic-ray spectrum suggests the presence of Galactic PeV proton accelerators called "PeVatrons". Supernova Remnant (SNR) G106.3+2.7 is a prime candidate for one of these. The recent detection of very high energy (0.1-100 TeV) gamma rays from G106.3+2.7 may be explained either by the decay of neutral pions or inverse Compton scattering by relativistic electrons. We report an analysis of 12 years of Fermi-LAT gamma-ray data which shows that the GeV-TeV gamma-ray spectrum is much harder and requires a different total electron energy than the radio and X-ray spectra, suggesting it has a distinct, hadronic origin. The non-detection of gamma rays below 10 GeV implies additional constraints on the relativistic electron spectrum. A hadronic interpretation of the observed gamma rays is strongly supported. This observation confirms the long-sought connection between Galactic PeVatrons and SNRs. Moreover, it suggests that G106.3+2.7 could be the brightest member of a new population of SNRs whose gamma-ray energy flux peaks at TeV energies. Such a population may contribute to the cosmic-ray knee and be revealed by future very high energy gamma-ray detectors.

preprint2022arXiv

Hard X-ray emission from the eastern jet of SS 433 powering the W50 `Manatee' nebula: Evidence for particle re-acceleration

We present a broadband X-ray study of W50 (`the Manatee nebula'), the complex region powered by the microquasar SS 433, that provides a test-bed for several important astrophysical processes. The W50 nebula, a Galactic PeVatron candidate, is classified as a supernova remnant but has an unusual double-lobed morphology likely associated with the jets from SS 433. Using NuSTAR, XMM-Newton, and Chandra observations of the inner eastern lobe of W50, we have detected hard non-thermal X-ray emission up to $\sim$30 keV, originating from a few-arcminute size knotty region (`Head') located $\lesssim$ 18$^{\prime}$ (29 pc for a distance of 5.5 kpc) east of SS 433, and constrain its photon index to 1.58$\pm$0.05 (0.5-30 keV band). The index gradually steepens eastward out to the radio `ear' where thermal soft X-ray emission with a temperature $kT$$\sim$0.2 keV dominates. The hard X-ray knots mark the location of acceleration sites within the jet and require an equipartition magnetic field of the order of $\gtrsim$12$μ$G. The unusually hard spectral index from the `Head' region challenges classical particle acceleration processes and points to particle injection and re-acceleration in the sub-relativistic SS 433 jet, as seen in blazars and pulsar wind nebulae.

preprint2022arXiv

The TeV Diffuse Cosmic Neutrino Spectrum and the Nature of Astrophysical Neutrino Sources

The diffuse flux of cosmic neutrinos has been measured by the IceCube Observatory from TeV to PeV energies. We show that an improved characterization of this flux at the lower energies, TeV and sub-TeV, reveals important information on the nature of the astrophysical neutrino sources in a model-independent way. Most significantly, it could confirm the present indications that neutrinos originate in cosmic environments that are optically thick to GeV-TeV $γ$-rays. This conclusion will become inevitable if a steeper or even uninterrupted neutrino power law is observed in the TeV region. In such $γ$-ray-obscured sources, the $γ$-rays that inevitably accompany cosmic neutrinos will cascade down to MeV-GeV energies. The requirement that the cascaded $γ$-ray flux accompanying cosmic neutrinos should not exceed the observed diffuse $γ$-ray background, puts constraints on the peak energy and density of the radiation fields in the sources. Our calculations inspired by the existing data suggest that a fraction of the observed diffuse MeV-GeV $γ$-ray background may be contributed by neutrino sources with intense radiation fields that obscure the high-energy $γ$-ray emission accompanying the neutrinos.

preprint2021arXiv

High-Energy Neutrinos and Gamma-Rays from Non-Relativistic Shock-Powered Transients

Shock interaction has been argued to play a role in powering a range of optical transients, including supernovae (particularly the superluminous class), classical novae, stellar mergers, tidal disruption events, and fast blue optical transients. These same shocks can accelerate relativistic ions, generating high-energy neutrino and gamma-ray emission via hadronic pion production. The recent discovery of time-correlated optical and gamma-ray emission in classical novae has revealed the important role of radiative shocks in powering these events, enabling an unprecedented view of the properties of ion acceleration, including its efficiency and energy spectrum, under similar physical conditions to shocks in extragalactic transients. Here we introduce a model for connecting the radiated optical fluence of non-relativistic transients to their maximal neutrino and gamma-ray fluence. We apply this technique to a wide range of extragalactic transient classes in order to place limits on their contributions to the cosmological high-energy gamma-ray and neutrino backgrounds. Based on a simple model for diffusive shock acceleration at radiative shocks, calibrated to novae, we demonstrate that several of the most luminous transients can accelerate protons up to energies $E_{\rm max} \gtrsim 10^{16}$ eV, sufficient to contribute to the IceCube astrophysical background. Furthermore, several of the considered sources$-$particularly hydrogen-poor supernovae$-$may serve as "hidden" gamma-ray sources due to the high gamma-ray opacity of their ejecta, evading constraints imposed by the non-blazar Fermi-LAT background. However, adopting an ion acceleration efficiency $\sim$ 0.3-1$\%$ motivated by nova observations, we find that currently known classes of non-relativistic, potentially shock-powered transients contribute at most a few percent of the total IceCube background.

preprint2020arXiv

A Cross-Correlation Study of High-energy Neutrinos and Tracers of Large-Scale Structure

The origin of the bulk of the astrophysical neutrinos detected by the IceCube Observatory remains a mystery. Previous source-finding analyses compare the directions of IceCube events and individual sources in astrophysical catalogs. The source association method is technically challenging when the number of source candidates is much larger than the number of the observed astrophysical neutrinos. We show that in this large source number regime, a cross-correlation analysis of neutrino data and source catalog can instead be used to constrain potential source populations for the high-energy astrophysical neutrinos, and provide spatial evidence for the existence of astrophysical neutrinos. We present an analysis of the cross-correlation of the IceCube 2010-2012 point-source data and a WISE-2MASS galaxy sample. While we find no significant detection of cross-correlation with the publicly available neutrino dataset, we show that, when applied to the full IceCube data, which has a longer observation time and higher astrophysical neutrino purity, our method has sufficient statistical power to detect a cross-correlation signal if the neutrino sources trace the Large Scale Structure of the Universe.

preprint2020arXiv

GeV-TeV Counterparts of SS 433/W50 from Fermi-LAT and HAWC Observations

The extended jets of the microquasar SS 433 have been observed in optical, radio, X-ray, and recently very-high-energy (VHE) $γ$-rays by HAWC. The detection of HAWC $γ$-rays with energies as great as 25 TeV motivates searches for high-energy $γ$-ray counterparts in the Fermi-LAT data in the 100 MeV--300 GeV band. In this paper, we report on the first-ever joint analysis of Fermi-LAT and HAWC observations to study the spectrum and location of $γ$-ray emission from SS~433. Our analysis finds common emission sites of GeV-to-TeV $γ$-rays inside the eastern and western lobes of SS 433. The total flux above 1 GeV is $\sim 1\times10^{-10}\,\rm cm^{-2}\,s^{-1}$ in both lobes. The $γ$-ray spectrum in the eastern lobe is consistent with inverse-Compton emission by an electron population that is accelerated by jets. To explain both the GeV and TeV flux, the electrons need to have a soft intrinsic energy spectrum, or undergo a quick cooling process due to synchrotron radiation in a magnetized environment.

preprint2020arXiv

High-energy neutrinos from fallback accretion of binary neutron star merger remnants

Following the coalescence of binary neutron stars, debris from the merger which remains marginally bound to the central compact remnant will fallback at late times, feeding a sustained accretion flow. Unbound winds or a wide-angle jet from this radiatively-inefficient disk may collide with the comparatively slow dense kilonova ejecta released from an earlier phase. Under the assumption that such interaction accelerate cosmic rays to ultra-high energies, we numerically simulate their propagation and interactions through the dynamical ejecta. The hadronuclear and photo-hadronic processes experienced by particles produce isotropic high-energy neutrino fluxes, peaking at times $10^{3-4}\,$s, which we calculate for two sets of parameters. A first set is inspired by the observations of GW170817. In the second scenario, which we call optimistic, parameters are chosen so as to optimize the neutrino flux, within the range allowed by observation and theory. We find that single sources can only be detected with IceCube-Gen2 for optimistic scenarios and if located within $\sim 4\,$Mpc. The cumulative flux could contribute to $\sim 0.5-10\%$ of the diffuse flux observed by the IceCube Observatory, depending on the fall-back power and the cosmic ray composition. The neutrino emission powered by fallback is nearly isotropic, and can be used for future correlation studies with gravitational wave signals.

preprint2020arXiv

Selective Transfer with Reinforced Transfer Network for Partial Domain Adaptation

One crucial aspect of partial domain adaptation (PDA) is how to select the relevant source samples in the shared classes for knowledge transfer. Previous PDA methods tackle this problem by re-weighting the source samples based on their high-level information (deep features). However, since the domain shift between source and target domains, only using the deep features for sample selection is defective. We argue that it is more reasonable to additionally exploit the pixel-level information for PDA problem, as the appearance difference between outlier source classes and target classes is significantly large. In this paper, we propose a reinforced transfer network (RTNet), which utilizes both high-level and pixel-level information for PDA problem. Our RTNet is composed of a reinforced data selector (RDS) based on reinforcement learning (RL), which filters out the outlier source samples, and a domain adaptation model which minimizes the domain discrepancy in the shared label space. Specifically, in the RDS, we design a novel reward based on the reconstruct errors of selected source samples on the target generator, which introduces the pixel-level information to guide the learning of RDS. Besides, we develope a state containing high-level information, which used by the RDS for sample selection. The proposed RDS is a general module, which can be easily integrated into existing DA models to make them fit the PDA situation. Extensive experiments indicate that RTNet can achieve state-of-the-art performance for PDA tasks on several benchmark datasets.

preprint2019arXiv

Ultra-High Energy Cosmic Rays and Neutrinos from Tidal Disruptions by Massive Black Holes

Tidal disruptions are extremely powerful phenomena that have been designated as candidate sources of ultra-high-energy cosmic rays. The disruption of a star by a black hole can naturally provide protons and heavier nuclei, which can be injected and accelerated to ultra-high energies within a jet. Inside the jet, accelerated nuclei are likely to interact with a dense photon field, leading to a significant production of neutrinos and secondary particles. We model numerically the propagation and interactions of high-energy nuclei in jetted tidal disruption events in order to evaluate consistently their signatures in cosmic rays and neutrinos. We propose a simple model of the light curve of tidal disruption events, consisting of two stages: a high state with bright luminosity and short duration and a medium state, less bright and longer lasting. These two states have different impacts on the production of cosmic rays and neutrinos. In order to calculate the diffuse fluxes of cosmic rays and neutrinos, we model the luminosity function and redshift evolution of jetted tidal disruption events. We find that we can fit the latest ultra-high-energy cosmic-ray spectrum and composition results of the Auger experiment for a range of reasonable parameters. The diffuse neutrino flux associated with this scenario is found to be subdominant, but nearby events can be detected by IceCube or next-generation detectors such as IceCube-Gen2.

preprint2016arXiv

A New Method for Finding Point Sources in High-energy Neutrino Data

The IceCube collaboration has reported the first detection of high-energy astrophysical neutrinos including $\sim 50$ high-energy starting events, but no individual sources have been identified. It is therefore important to develop the most sensitive and efficient possible algorithms to identify point sources of these neutrinos. The most popular current method works by exploring a dense grid of possible directions to individual sources, and identifying the single direction with the maximum probability of having produced multiple detected neutrinos. This method has numerous strengths, but it is computationally intensive and, because it focuses on the single best location for a point source, additional point sources are not included in the evidence. We propose a new maximum likelihood method that uses the angular separations between all pairs of neutrinos in the data. Unlike existing autocorrelation methods for this type of analysis, which also use angular separations between neutrino pairs, our method incorporates information about the point spread function and can identify individual point sources. We find that if the angular resolution is a few degrees or better, then this approach reduces both false positive and false negative errors compared to the current method, and is also more computationally efficient up to, potentially, hundreds of thousands of detected neutrinos.

preprint2016arXiv

Cluster Mergers and the Origin of the ARCADE-2 Excess

Radio observations at multiple frequencies have detected a significant isotropic emission component between 22~MHz and 10~GHz, commonly termed the ARCADE-2 Excess. The origin of this radio emission is unknown, as the intensity, spectrum and isotropy of the signal are difficult to model with either traditional astrophysical mechanisms or novel physics such as dark matter annihilation. We posit a new model capable of explaining the key components of the excess radio emission. Specifically, we show that the re-acceleration of non-thermal electrons via turbulence in merging galaxy clusters are capable of explaining the intensity, spectrum, and isotropy of the ARCADE-2 data. We examine the parameter spaces of cluster re-acceleration, magnetic field, and merger rate, finding that the radio excess can be reproduced assuming reasonable assumptions for each. Finally, we point out that future observations will definitively confirm or rule-out the contribution of cluster mergers to the isotropic radio background.

preprint2016arXiv

Energy Estimation of Cosmic Rays with the Engineering Radio Array of the Pierre Auger Observatory

The Auger Engineering Radio Array (AERA) is part of the Pierre Auger Observatory and is used to detect the radio emission of cosmic-ray air showers. These observations are compared to the data of the surface detector stations of the Observatory, which provide well-calibrated information on the cosmic-ray energies and arrival directions. The response of the radio stations in the 30 to 80 MHz regime has been thoroughly calibrated to enable the reconstruction of the incoming electric field. For the latter, the energy deposit per area is determined from the radio pulses at each observer position and is interpolated using a two-dimensional function that takes into account signal asymmetries due to interference between the geomagnetic and charge-excess emission components. The spatial integral over the signal distribution gives a direct measurement of the energy transferred from the primary cosmic ray into radio emission in the AERA frequency range. We measure 15.8 MeV of radiation energy for a 1 EeV air shower arriving perpendicularly to the geomagnetic field. This radiation energy -- corrected for geometrical effects -- is used as a cosmic-ray energy estimator. Performing an absolute energy calibration against the surface-detector information, we observe that this radio-energy estimator scales quadratically with the cosmic-ray energy as expected for coherent emission. We find an energy resolution of the radio reconstruction of 22% for the data set and 17% for a high-quality subset containing only events with at least five radio stations with signal.

preprint2016arXiv

High-energy neutrinos from sources in clusters of galaxies

High-energy cosmic rays can be accelerated in clusters of galaxies, by mega-parsec scale shocks induced by accretion of gas during the formation of large-scale structure, or by powerful sources harbored in clusters. Once accelerated, the highest energy particles leave the cluster via almost rectilinear trajectories, while lower energy ones can be confined by the cluster magnetic field up to cosmological time and interact with the intracluster gas. Using a realistic model of the baryon distribution and the turbulent magnetic field in clusters, we studied the propagation and hadronic interaction of high-energy protons in the intracluster medium. We report the cumulative cosmic ray and neutrino spectra generated by galaxy clusters including embedded sources, and demonstrate that clusters can contribute a significant fraction of the observed IceCube neutrinos above 30 TeV while remaining undetected in high-energy cosmic rays and $γ$ rays for reasonable choices of parameters and source scenarios.

preprint2016arXiv

IceCube Constraints on Fast-Spinning Pulsars as High-Energy Neutrino Sources

Relativistic winds of fast-spinning pulsars have been proposed as a potential site for cosmic-ray acceleration from very high energies (VHE) to ultrahigh energies (UHE). We re-examine conditions for high-energy neutrino production, considering the interaction of accelerated particles with baryons of the expanding supernova ejecta and the radiation fields in the wind nebula. We make use of the current IceCube sensitivity in diffusive high-energy neutrino background, in order to constrain the parameter space of the most extreme neutron stars as sources of VHE and UHE cosmic rays. We demonstrate that the current non-observation of $10^{18}$ eV neutrinos put stringent constraints on the pulsar scenario. For a given model, birthrates, ejecta mass and acceleration efficiency of the magnetar sources can be constrained. When we assume a proton cosmic ray composition and spherical supernovae ejecta, we find that the IceCube limits almost exclude their significant contribution to the observed UHE cosmic-ray flux. Furthermore, we consider scenarios where a fraction of cosmic rays can escape from jet-like structures piercing the ejecta, without significant interactions. Such scenarios would enable the production of UHE cosmic rays and help remove the tension between their EeV neutrino production and the observational data.

preprint2016arXiv

Measurement of the Radiation Energy in the Radio Signal of Extensive Air Showers as a Universal Estimator of Cosmic-Ray Energy

We measure the energy emitted by extensive air showers in the form of radio emission in the frequency range from 30 to 80 MHz. Exploiting the accurate energy scale of the Pierre Auger Observatory, we obtain a radiation energy of 15.8 \pm 0.7 (stat) \pm 6.7 (sys) MeV for cosmic rays with an energy of 1 EeV arriving perpendicularly to a geomagnetic field of 0.24 G, scaling quadratically with the cosmic-ray energy. A comparison with predictions from state-of-the-art first-principle calculations shows agreement with our measurement. The radiation energy provides direct access to the calorimetric energy in the electromagnetic cascade of extensive air showers. Comparison with our result thus allows the direct calibration of any cosmic-ray radio detector against the well-established energy scale of the Pierre Auger Observatory.

preprint2016arXiv

Testing the Newborn Pulsar Origin of Ultrahigh Energy Cosmic Rays with EeV Neutrinos

Fast-spinning newborn pulsars are intriguing candidate sources of ultrahigh energy cosmic rays (UHECRs). The acceleration of particles with a given composition in a fraction of the extragalactic pulsar population can give a consistent explanation for the measurements of the Auger Observatory. We calculate the associated diffuse neutrino flux produced while particles cross the supernova ejecta surrounding the stars. We show that in the minimal pulsar scenario that are compatible with the UHECR data, the effective optical depth to hadronuclear interactions is larger than unity at ultrahigh energies. Thus, even in the most pessimistic case, one expects energy fluxes of $\sim 0.1 - 1 \,\rm EeV$ neutrinos that should be detectable with IceCube or Askaryan Radio Array within a decade.

preprint2016arXiv

The Highest-Energy Cosmic Rays Cannot be Dominantly Protons from Steady Sources

The bulk of observed ultrahigh energy cosmic rays could be light or heavier elements, and originate from an either steady or transient population of sources. This leaves us with four general categories of sources. Energetic requirements set a lower limit on single source luminosities, while the distribution of particle arrival directions in the sky sets a lower limit on the source number density. The latter constraint depends on the angular smearing in the skymap due to the magnetic deflections of the charged particles during their propagation from the source to the Earth. We contrast these limits with the luminosity functions from surveys of existing luminous steady objects in the nearby universe, and strongly constrain one of the four categories of source models, namely, steady proton sources. The possibility that cosmic rays with energy $> 8\times 10^{19}\,$eV are dominantly pure protons coming from steady sources is excluded at 95\% confidence level, under the safe assumption that protons experience less than $30^\circ$ magnetic deflection on flight.

preprint2015arXiv

An improved limit to the diffuse flux of ultra-high energy neutrinos from the Pierre Auger Observatory

Neutrinos in the cosmic ray flux with energies near 1 EeV and above are detectable with the Surface Detector array of the Pierre Auger Observatory. We report here on searches through Auger data from 1 January 2004 until 20 June 2013. No neutrino candidates were found, yielding a limit to the diffuse flux of ultra-high energy neutrinos that challenges the Waxman-Bahcall bound predictions. Neutrino identification is attempted using the broad time-structure of the signals expected in the SD stations, and is efficiently done for neutrinos of all flavors interacting in the atmosphere at large zenith angles, as well as for "Earth-skimming" neutrino interactions in the case of tau neutrinos. In this paper the searches for downward-going neutrinos in the zenith angle bins $60^\circ-75^\circ$ and $75^\circ-90^\circ$ as well as for upward-going neutrinos, are combined to give a single limit. The $90\%$ C.L. single-flavor limit to the diffuse flux of ultra-high energy neutrinos with an $E^{-2}$ spectrum in the energy range $1.0 \times 10^{17}$ eV - $2.5 \times 10^{19}$ eV is $E_ν^2 dN_ν/dE_ν< 6.4 \times 10^{-9}~ {\rm GeV~ cm^{-2}~ s^{-1}~ sr^{-1}}$.

preprint2015arXiv

High-Energy Neutrino Signatures of Newborn Pulsars In the Local Universe

Charged particles can be accelerated to higher than PeV energies in the electro-magnetic wind of a fast-spinning newborn pulsar to produce high-energy neutrinos, through hadronuclear interactions in the supernova remnant. Here we explore the detectability and observational signatures of these high-energy neutrinos. We show that their spectral index varies approximately from 1.5 to 2, depending on the relevant pulsar properties and observation time. We also apply the scenario to existing young pulsars in the local universe and find the corresponding neutrino flux well below current detection limits. Finally, we estimate the birth rate of fast-spinning pulsars in the local Universe that can be observed by the IceCube observatory to be 0.07 per year, with an upper limit of 0.29 per year.

preprint2015arXiv

Measurement of the cosmic ray spectrum above $4{\times}10^{18}$ eV using inclined events detected with the Pierre Auger Observatory

A measurement of the cosmic-ray spectrum for energies exceeding $4{\times}10^{18}$ eV is presented, which is based on the analysis of showers with zenith angles greater than $60^{\circ}$ detected with the Pierre Auger Observatory between 1 January 2004 and 31 December 2013. The measured spectrum confirms a flux suppression at the highest energies. Above $5.3{\times}10^{18}$ eV, the "ankle", the flux can be described by a power law $E^{-γ}$ with index $γ=2.70 \pm 0.02 \,\text{(stat)} \pm 0.1\,\text{(sys)}$ followed by a smooth suppression region. For the energy ($E_\text{s}$) at which the spectral flux has fallen to one-half of its extrapolated value in the absence of suppression, we find $E_\text{s}=(5.12\pm0.25\,\text{(stat)}^{+1.0}_{-1.2}\,\text{(sys)}){\times}10^{19}$ eV.

preprint2015arXiv

The Giant Radio Array for Neutrino Detection

High-energy neutrino astronomy will probe the working of the most violent phenomena in the Universe. The Giant Radio Array for Neutrino Detection (GRAND) project consists of an array of $\sim10^5$ radio antennas deployed over $\sim$200000km$^2$ in a mountainous site. It aims at detecting high-energy neutrinos via the measurement of air showers induced by the decay in the atmosphere of $τ$ leptons produced by the interaction of the cosmic neutrinos under the Earth surface. Our objective with GRAND is to reach a neutrino sensitivity of $3\times10^{-11}E^{-2}$GeV$^{-1}$cm$^{-2}$s$^{-1}$sr$^{-1}$ above $3 \times10^{16}$eV. This sensitivity ensures the detection of cosmogenic neutrinos in the most pessimistic source models, and about 100 events per year are expected for the standard models. GRAND would also probe the neutrino signals produced at the potential sources of UHECRs. We show how our preliminary design should enable us to reach our sensitivity goals, and present the experimental characteristics. We assess the possibility to adapt GRAND to other astrophysical radio measurements. We discuss in this token the technological options for the detector and the steps to be taken to achieve the GRAND project.

preprint2014arXiv

Is the Ultra-High Energy Cosmic-Ray Excess Observed by the Telescope Array Correlated with IceCube Neutrinos?

The Telescope Array (TA) has observed a statistically significant excess in cosmic-rays with energies above 57 EeV in a region of approximately 1150 square degrees centered on coordinates (R.A. = 146.7, Dec. = 43.2). We note that the location of this excess correlates with two of the 28 extraterrestrial neutrinos recently observed by IceCube. The overlap between the two IceCube neutrinos and the TA excess is statistically significant at the 2$σ$ level. Furthermore, the spectrum and intensity of the IceCube neutrinos is consistent with a single source which would also produce the TA excess. Finally, we discuss possible source classes with the correct characteristics to explain the cosmic-ray and neutrino fluxes with a single source.

preprint2014arXiv

The Anisotropy of the Extragalactic Radio Background from Dark Matter Annihilation

Observations of the extragalactic radio background have uncovered a significant isotropic emission across multiple frequencies spanning from 22 MHz to 10 GHz. The intensity of this non-thermal emission component significantly exceeds the expected contribution from known astrophysical sources. Interestingly, models have indicated that the annihilation of dark matter particles may reproduce both the flux and spectrum of the excess. However, the lack of a measurable anisotropy in the residual emission remains challenging for both dark matter and standard astrophysical interpretations of the ARCADE-2 data. We calculate the expected synchrotron anisotropy from dark matter annihilation and show that these models can produce very small anisotropies, though this requires galaxy clusters to have large substructure contributions and strong magnetic fields. We show that this constraint can be significantly relaxed, however, in scenarios where electrons produced via dark matter annihilation can be efficiently reaccelerated by Alfvén waves in the intra-Cluster medium. Our analysis indicates that any source capable of explaining the intensity and isotropy of the extragalactic radio excess must have a spatial extension far larger than typical for baryons in galaxies, suggesting a novel physics interpretation.

preprint2013arXiv

Ultrahigh Energy Cosmic Ray Nuclei from Extragalactic Pulsars and the effect of their Galactic counterparts

The acceleration of ultrahigh energy nuclei in fast spinning newborn pulsars can explain the observed spectrum of ultrahigh energy cosmic rays and the trend towards heavier nuclei for energies above $10^{19}\,$eV as reported by the Auger Observatory. Pulsar acceleration implies a hard injection spectrum ($\sim E^{-1}$) due to pulsar spin down and a maximum energy $E_{\rm max} \sim Z \, 10^{19}$ eV due to the limit on the spin rate of neutron stars. We have previously shown that the escape through the young supernova remnant softens the spectrum, decreases slightly the maximum energy, and generates secondary nuclei. Here we show that the distribution of pulsar birth periods and the effect of propagation in the interstellar and intergalactic media modifies the combined spectrum of all pulsars. By assuming a normal distribution of pulsar birth periods centered at 300 ms, we show that the contribution of extragalactic pulsar births to the ultrahigh energy cosmic ray spectrum naturally gives rise to a contribution to very high energy cosmic rays (VHECRs, between $10^{16}$ and $10^{18}$ eV) by Galactic pulsar births. The required injected composition to fit the observed spectrum depends on the absolute energy scale, which is uncertain, differing between Auger Observatory and Telescope Array. The contribution of Galactic pulsar births can also bridge the gap between predictions for cosmic ray acceleration in supernova remnants and the observed spectrum just below the ankle, depending on the composition of the cosmic rays that escape the supernova remnant and the diffusion behavior of VHECRs in the Galaxy.

preprint2012arXiv

Newly-born pulsars as sources of ultrahigh energy cosmic rays

Newly-born pulsars offer favorable sites for the injection of heavy nuclei, and for their further acceleration to ultrahigh energies. Once accelerated in the pulsar wind, nuclei have to escape from the surrounding supernova envelope. We examine this escape analytically and numerically, and discuss the pulsar source scenario in light of the latest ultrahigh energy cosmic ray (UHECR) data. Our calculations show that, at early times, when protons can be accelerated to energies E>10^20 eV, the young supernova shell tends to prevent their escape. In contrast, because of their higher charge, iron-peaked nuclei are still accelerated to the highest observed energies at later times, when the envelope has become thin enough to allow their escape. Ultrahigh energy iron nuclei escape newly-born pulsars with millisecond periods and dipole magnetic fields of ~10^(12-13) G, embedded in core-collapse supernovae. Due to the production of secondary nucleons, the envelope crossing leads to a transition of composition from light to heavy elements at a few EeV, as observed by the Auger Observatory. The escape also results in a softer spectral slope than that initially injected via unipolar induction, which allows for a good fit to the observed UHECR spectrum. We conclude that the acceleration of iron-peaked elements in a reasonably small fraction (< 0.01%) of extragalactic rotation-powered young pulsars would reproduce satisfactorily the current UHECR data. Possible signatures of this scenario are also discussed.