Researcher profile

Shunsaku Horiuchi

Shunsaku Horiuchi contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
20works
0followers
9topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

20 published item(s)

preprint2026arXiv

The Cosmic Neutrino Background is within Reach of Future Neutrino Telescopes

The cosmic neutrino background (C$ν$B) can be boosted to high energies due to scatterings with energetic cosmic rays (CRs) across cosmological scales. Previous calculations focused on neutral current incoherent and coherent elastic scatterings of cosmic-ray protons off relic neutrinos. However, charged current interactions and deep inelastic scatterings are also expected to occur, which enhances the boosted relic neutrino fluxes on Earth. Here, we compute the \textit{total} diffuse boosted cosmic neutrino background (DBC$ν$B) arising from CRs at all redshifts in the Universe, accounting for neutral current and charged current elastic and deep inelastic scatterings. We find that IceCube already places an upper limit on the cosmic neutrino background overdensity in cosmological scales of ~$\mathcal{O}(100-1000)$ at $E_ν=10^{10}$ GeV, for a lightest neutrino mass of $m_ν \gtrsim 0.1$ eV. We further show that IceCube-Gen2 could test $\mathcal{O}(1-10)$ C$ν$B overdensities, and the combination of $10$ future neutrino telescopes with similar sensitivity would allow us to test the $Λ$CDM expected C$ν$B density for a lightest neutrino mass compatible with the KATRIN bound.

preprint2023arXiv

Long-Exposure NuSTAR Constraints on Decaying Dark Matter in the Galactic Halo

We present two complementary NuSTAR x-ray searches for keV-scale dark matter decaying to mono-energetic photons in the Milky Way halo. In the first, we utilize the known intensity pattern of unfocused stray light across the detector planes -- the dominant source of photons from diffuse sources -- to separate astrophysical emission from internal instrument backgrounds using ${\sim}$7-Ms/detector deep blank-sky exposures. In the second, we present an updated parametric model of the full NuSTAR instrument background, allowing us to leverage the statistical power of an independent ${\sim}$20-Ms/detector stacked exposures spread across the sky. Finding no evidence of anomalous x-ray lines using either method, we set limits on the active-sterile mixing angle $\sin^2(2θ)$ for sterile-neutrino masses 6--40 keV. The first key result is that we strongly disfavor a ${\sim}$7-keV sterile neutrino decaying into a 3.5-keV photon. The second is that we derive leading limits on sterile neutrinos with masses ${\sim}$15--18 keV and ${\sim}$25--40 keV, reaching or extending below the Big Bang Nucleosynthesis limit. In combination with previous results, the parameter space for the Neutrino Minimal Standard Model ($ν$MSM) is now nearly closed.

preprint2022arXiv

Gamma-ray emission from the Sagittarius Dwarf Spheroidal galaxy due to millisecond pulsars

The Fermi Bubbles are giant, gamma-ray emitting lobes emanating from the nucleus of the Milky Way discovered in ~1-100 GeV data collected by the Large Area Telescope on board the Fermi Gamma-Ray Space Telescope. Previous work has revealed substructure within the Fermi Bubbles that has been interpreted as a signature of collimated outflows from the Galaxy's super-massive black hole. Here we show via a spatial template analysis that much of the gamma-ray emission associated to the brightest region of substructure -- the so-called cocoon -- is likely due to the Sagittarius dwarf spheroidal (Sgr dSph) galaxy. This large Milky Way satellite is viewed through the Fermi Bubbles from the position of the Solar System. As a tidally and ram-pressure stripped remnant, the Sgr dSph has no on-going star formation, but we nevertheless demonstrate that the dwarf's millisecond pulsar (MSP) population can plausibly supply the gamma-ray signal that our analysis associates to its stellar template. The measured spectrum is naturally explained by inverse Compton scattering of cosmic microwave background photons by high-energy electron-positron pairs injected by MSPs belonging to the Sgr dSph, combined with these objects' magnetospheric emission. This finding plausibly suggests that MSPs produce significant gamma-ray emission amongst old stellar populations, potentially confounding indirect dark matter searches in regions such as the Galactic Centre, the Andromeda galaxy, and other massive Milky Way dwarf spheroidals.

preprint2022arXiv

Shock Revival in Core-collapse Supernovae Assisted by Heavy Axion-like Particles

Axion-like particles (ALPs) are a class of hypothetical pseudoscalar particles which feebly interact with ordinary matter. The hot plasma of core-collapse supernovae is a possible laboratory to explore physics beyond the standard model including ALPs. Once produced, some of the ALPs can be absorbed by the supernova matter and affect energy transfer. In this study, we calculate the ALP emission in core-collapse supernovae and the backreaction on supernova dynamics consistently. It is found that the stalled bounce shock can be revived if the coupling between ALPs and photons is as high as g_{ag}~10^{-9} GeV^{-1} and the ALP mass is 40-400 MeV. Most of the models result in more energetic explosions than the average observed supernova. While this can be used to place constraints on those ALPs, long-term simulations across multiple progenitors need to be further investigated to place robust limits.

preprint2022arXiv

Snowmass2021 Cosmic Frontier White Paper: Puzzling Excesses in Dark Matter Searches and How to Resolve Them

Intriguing signals with excesses over expected backgrounds have been observed in many astrophysical and terrestrial settings, which could potentially have a dark matter origin. Astrophysical excesses include the Galactic Center GeV gamma-ray excess detected by the Fermi Gamma-Ray Space Telescope, the AMS antiproton and positron excesses, and the 511 and 3.5 keV X-ray lines. Direct detection excesses include the DAMA/LIBRA annual modulation signal, the XENON1T excess, and low-threshold excesses in solid state detectors. We discuss avenues to resolve these excesses, with actions the field can take over the next several years.

preprint2022arXiv

Snowmass2021 Cosmic Frontier: Synergies between dark matter searches and multiwavelength/multimessenger astrophysics

This whitepaper focuses on the astrophysical systematics which are encountered in dark matter searches. Oftentimes in indirect and also in direct dark matter searches, astrophysical systematics are a major limiting factor to sensitivity to dark matter. Just as there are many forms of dark matter searches, there are many forms of backgrounds. We attempt to cover the major systematics arising in dark matter searches using photons -- radio and gamma rays -- to cosmic rays, neutrinos and gravitational waves. Examples include astrophysical sources of cosmic messengers and their interactions which can mimic dark matter signatures. In turn, these depend on commensurate studies in understanding the cosmic environment -- gas distributions, magnetic field configurations -- as well as relevant nuclear astrophysics. We also cover the astrophysics governing celestial bodies and galaxies used to probe dark matter, from black holes to dwarf galaxies. Finally, we cover astrophysical backgrounds related to probing the dark matter distribution and kinematics, which impact a wide range of dark matter studies. In the future, the rise of multi-messenger astronomy, and novel analysis methods to exploit it for dark matter, will offer various strategic ways to continue to enhance our understanding of astrophysical backgrounds to deliver improved sensitivity to dark matter.

preprint2022arXiv

Systematic exploration of heavy element nucleosynthesis in protomagnetar outflows

We study the nucleosynthesis products in neutrino-driven winds from rapidly rotating, highly magnetised and misaligned protomagnetars using the nuclear reaction network SkyNet. We adopt a semi-analytic parametrized model for the protomagnetar and systematically study the capabilities of its neutrino-driven wind for synthesizing nuclei and eventually producing ultra-high energy cosmic rays (UHECRs). We find that for neutron-rich outflows ($Y_e<0.5$), synthesis of heavy elements ($\overline{A}\sim 20-65$) is possible during the first $\sim 10$ seconds of the outflow, but these nuclei are subjected to composition-altering photodisintegration during the epoch of particle acceleration at the dissipation radii. However, after the first $\sim 10$ seconds of the outflow, nucleosynthesis reaches lighter elements ($\overline{A}\sim 10-50$) that are not subjected to subsequent photodisintegration. For proton-rich ($Y_e \geq 0.5$) outflows, synthesis is more limited ($\overline{A}\sim 4-15$). These suggest that while protomagnetars typically do not synthesize nuclei heavier than second r-process peak elements, they are intriguing sources of intermediate/heavy mass UHECRs. For all configurations, the most rapidly rotating protomagnetars are more conducive for nucleosynthesis with a weaker dependence on the magnetic field strength.

preprint2022arXiv

The Andromeda Gamma-Ray Excess: Background Systematics of the Millisecond Pulsars and Dark Matter Interpretations

Since the discovery of an excess in gamma rays in the direction of M31, its cause has been unclear. Published interpretations focus on a dark matter or stellar related origin. Studies of a similar excess in the Milky Way center motivate a correlation of the spatial morphology of the signal with the distribution of stellar mass in M31. However, a robust determination of the best theory for the observed excess emission is very challenging due to large uncertainties in the astrophysical gamma-ray foreground model. Here we perform a spectro-morphological analysis of the M31 gamma-ray excess using state-of-the-art templates for the distribution of stellar mass in M31 and novel astrophysical foreground models for its sky region. We construct maps for the old stellar populations of M31 based on observational data from the PAndAS survey and carefully remove the foreground stars. We also produce improved astrophysical foreground models by using novel image inpainting techniques based on machine learning methods. We find that our stellar maps, taken as a proxy for the location of a putative population of millisecond pulsars in the bulge of M31, reach a statistical significance of $5.4σ$, making them as strongly favoured as the simple phenomenological models usually considered in the literature, e.g., a disk-like template with uniform brightness. Our detection of the stellar templates is robust to generous variations of the astrophysical foreground model. Once the stellar templates are included in the astrophysical model, we show that the dark matter annihilation interpretation of the signal is unwarranted. Using the results of a binary population synthesis model we demonstrate that a population of about one million unresolved MSPs could naturally explain the observed gamma-ray luminosity per stellar mass, energy spectrum, and stellar bulge-to-disk flux ratio.

preprint2021arXiv

Impact of binary interactions on the diffuse supernova neutrino background

Binary interactions, especially mass transfer and mergers, can strongly influence the evolution of massive stars and change their final properties and the occurrence of supernovae. Here, we investigate how binary interactions affect predictions of the diffuse flux of neutrinos. By performing stellar population syntheses including prescriptions for binary interactions, we show that the resulting detection rates of the diffuse supernova neutrino background is enhanced by 15%-20% compared to estimates without binary considerations. A source of significant uncertainty arises due to the presently sparse knowledge of the evolution of rapidly rotating carbon-oxygen cores, especially those created as a result of mergers near the white dwarf to core collapse boundary. The enhancement effect may be as small as a few percent if the effects of rotation in postmerger systems are neglected, or as large as 75% if trends are extrapolated. Our estimates serve to highlight that binary effects can be important.

preprint2021arXiv

Strong constraints on thermal relic dark matter from Fermi-LAT observations of the Galactic Center

The extended excess toward the Galactic Center (GC) in gamma rays inferred from Fermi-LAT observations has been interpreted as being due to dark matter (DM) annihilation. Here, we perform new likelihood analyses of the GC and show that, when including templates for the stellar galactic and nuclear bulges, the GC shows no significant detection of a DM annihilation template, even after generous variations in the Galactic diffuse emission models and a wide range of DM halo profiles. We include Galactic diffuse emission models with combinations of three-dimensional inverse Compton maps, variations of interstellar gas maps, and a central source of electrons. For the DM profile, we include both spherical and ellipsoidal DM morphologies and a range of radial profiles from steep cusps to kiloparsec-sized cores, motivated in part by hydrodynamical simulations. Our derived upper limits on the dark matter annihilation flux place strong constraints on DM properties. In the case of the pure $b$-quark annihilation channel, our limits on the annihilation cross section are more stringent than those from the Milky Way dwarfs up to DM masses of approximately TeV and rule out the thermal relic cross section up to approximately 300 GeV. Better understanding of the DM profile, as well as the Fermi-LAT data at its highest energies, would further improve the sensitivity to DM properties.

preprint2020arXiv

Cross Correlation of the Extragalactic Gamma-ray Background with Thermal Sunyaev-Zel&#39;dovich Effect in the Cosmic Microwave Background

Cosmic rays in galaxy clusters are unique probes of energetic processes operating with large-scale structures in the Universe. Precise measurements of cosmic rays in galaxy clusters are essential for improving our understanding of non-thermal components in the intracluster-medium (ICM) as well as the accuracy of cluster mass estimates in cosmological analyses. In this paper, we perform a cross-correlation analysis with the extragalactic gamma-ray background and the thermal Sunyaev-Zel&#39;dovich (tSZ) effect in the cosmic microwave background. The expected cross-correlation signal would contain rich information about the cosmic-ray-induced gamma-ray emission in the most massive galaxy clusters at $z\sim0.1-0.2$. We analyze the gamma-ray background map with 8 years of data taken by the Large Area Telescope onboard Fermi satellite and the publicly available tSZ map by Planck. We confirm that the measured cross-correlation is consistent with a null detection, and thus it enables us to put the tightest constraint on the acceleration efficiency of cosmic ray protons at shocks in and around galaxy clusters. We find the acceleration efficiency must be below 5\% with a $2σ$ confidence level when the hydrostatic mass bias of clusters is assumed to be 30\%, and our result is not significantly affected by the assumed value of the hydrostatic mass bias. Our constraint implies that the non-thermal cosmic-ray pressure in the ICM can introduce only a $\le 3\%$ level of the hydrostatic mass bias, highlighting that cosmic rays alone do not account for the mass bias inferred by the Planck analyses. Finally, we discuss future detectability prospects of cosmic-ray-induced gamma rays from the Perseus cluster for the Cherenkov Telescope Array.

preprint2020arXiv

Detectability of Collective Neutrino Oscillation Signatures in the Supernova Explosion of a 8.8 $M_\odot$ star

In order to investigate the impact of collective neutrino oscillations (CNO) on the neutrino signal from a nearby supernova, we perform 3-flavor neutrino oscillation simulations employing the multiangle effect. The background hydrodynamic model is based on the neutrino hydrodynamic simulation of a 8.8 \Msun progenitor star. We find that CNO commences after some 100 ms post bounce. Before this, CNO is suppressed by matter-induced decoherence. In the inverted mass hierarchy, the spectrum of $\barν_e$ becomes softer after the onset of CNO. To evaluate the detectability of this modification, we define a hardness ratio between the number of high energy neutrino events and low energy neutrino events adopting a fixed critical energy. We show that Hyper-Kamiokande (HK) can distinguish the effect of CNO for supernova distances out to $\sim 10$ kpc. On the other hand, for the normal mass hierarchy, the spectrum of $ν_e$ becomes softer after the onset of CNO, and we show that DUNE can distinguish this feature for supernova distances out to $\sim 10$ kpc. More work is necessary to optimize the best value of critical energy for maximum sensitivity. We also show that if the spectrum of $\barν_e$ in HK becomes softer due to CNO, the spectrum of $ν_e$ in DUNE becomes harder, and vice versa. This synergistic observations in $\barν_e$ and $ν_e$, by HK and DUNE respectively, will be an intriguing opportunity to test the occurrence of CNO.

preprint2020arXiv

Neutrino halo effect on collective neutrino oscillation in iron core-collapse supernova model of a 9.6 $M_{\odot}$ star

We extend the multi-angle computational framework and investigate the time evolution of the neutrino halo on collective neutrino oscillation in the core collapse of an iron core progenitor. We find that in the case of the $9.6\, \rm M_\odot$ progenitor adopted in this work, there are windows of time when the effects of neutrino halo and collective neutrino oscillation are not simultaneously large. Inside the shock, the impact of the inward-scattered halo neutrino cannot in general be neglected compared to the outward-propagating neutrino flux. However, during early epochs, collective neutrino oscillation is effectively shut down by multi-angle matter suppression. During the intermediate epoch, collective neutrino oscillation is not suppressed, but its onset radius is beyond the still relatively small explosion shock front where the halo is prominent. We also find in the case of the $9.6\, \rm M_\odot$ progenitor the halo neutrinos induce a delay in the onset of collective neutrino oscillations. This causes novel flavor conversions which sharpen collective neutrino oscillation spectral features. We predict that the inclusion of neutrino halo effects makes neutrino signals that are more clearly distinct from thermal emission that when halo neutrinos are omitted.

preprint2020arXiv

NuSTAR Tests of Sterile-Neutrino Dark Matter: New Galactic Bulge Observations and Combined Impact

We analyze two dedicated NuSTAR observations with exposure ${\sim}190$ ks located ${\sim}10^\circ$ from the Galactic plane, one above and the other below, to search for x-ray lines from the radiative decay of sterile-neutrino dark matter. These fields were chosen to minimize astrophysical x-ray backgrounds while remaining near the densest region of the dark matter halo. We find no evidence of anomalous x-ray lines in the energy range 5--20 keV, corresponding to sterile neutrino masses 10--40 keV. Interpreted in the context of sterile neutrinos produced via neutrino mixing, these observations provide the leading constraints in the mass range 10--12 keV, improving upon previous constraints in this range by a factor ${\sim}2$. We also compare our results to Monte Carlo simulations, showing that the fluctuations in our derived limit are not dominated by systematic effects. An updated model of the instrumental background, which is currently under development, will improve NuSTAR&#39;s sensitivity to anomalous x-ray lines, particularly for energies 3--5 keV.

preprint2020arXiv

The Predicted Properties of Helium-Enriched Globular Cluster Progenitors at High Redshift

Globular cluster progenitors may have been detected by \textit{HST}, and are predicted to be observable with \textit{JWST} and ground-based extremely-large telescopes with adaptive optics. This has the potential to elucidate the issue of globular cluster formation and the origins of significantly helium-enriched subpopulations, a problem in Galactic astronomy with no satisfactory theoretical solution. Given this context, we use model stellar tracks and isochrones to investigate the predicted observational properties of helium-enriched stellar populations in globular cluster progenitors. We find that, relative to helium-normal populations, helium-enriched ($ΔY=+0.12$) stellar populations similar to those inferred in the most massive globular clusters, are expected, modulo some rapid fluctuations in the first $\sim$30 Myr, to be brighter and redder in the rest frame. At fixed age, stellar mass, and metallicity, a helium-enriched population is predicted to converge to being $\sim$0.40 mag brighter at $λ\approx 2.0\, μm$, and to be 0.30 mag redder in the \textit{JWST}-NIRCam colour $(F070W-F200W)$, and to actually be fainter for $λ\lesssim 0.50 \, μm$. Separately, we find that the time-integrated shift in ionizing radiation is a negligible $\sim 5\%$, though we show that the Lyman-$α$ escape fraction could end up higher for helium-enriched stars.

preprint2020arXiv

Time of Flight and Supernova Progenitor Effects on the Neutrino Halo

We argue that the neutrino halo, a population of neutrinos that have undergone direction-changing scattering in the stellar envelope of a core-collapse supernova (CCSNe), is sensitive to neutrino emission history through time of flight. We show that the constant time approximation, commonly used in calculating the neutrino halo, does not capture the spatiotemporal evolution of the halo neutrino population and that correcting for time of flight can produce conditions which may trigger fast neutrino flavor conversion. We also find that there exists a window of time early in all CCSNe where the neutrino halo population is sufficiently small that it may be negligible. This suggests that collective neutrino oscillation calculations which neglect the Halo may be well founded at sufficiently early times.

preprint2011arXiv

The Cosmic Core-collapse Supernova Rate does not match the Massive-Star Formation Rate

We identify a &#34;supernova rate problem&#34;: the measured cosmic core-collapse supernova rate is a factor of ~ 2 smaller (with significance ~ 2 sigma) than that predicted from the measured cosmic massive-star formation rate. The comparison is critical for topics from galaxy evolution and enrichment to the abundance of neutron stars and black holes. We systematically explore possible resolutions. The accuracy and precision of the star formation rate data and conversion to the supernova rate are well supported, and proposed changes would have far-reaching consequences. The dominant effect is likely that many supernovae are missed because they are either optically dim (low-luminosity) or dark, whether intrinsically or due to obscuration. We investigate supernovae too dim to have been discovered in cosmic surveys by a detailed study of all supernova discoveries in the local volume. If possible supernova impostors are included, then dim supernovae are common enough by fraction to solve the supernova rate problem. If they are not included, then the rate of dark core collapses is likely substantial. Other alternatives are that there are surprising changes in our understanding of star formation or supernova rates, including that supernovae form differently in small galaxies than in normal galaxies. These possibilities can be distinguished by upcoming supernova surveys, star formation measurements, searches for disappearing massive stars, and measurements of supernova neutrinos.

preprint2011arXiv

The Production of Ultra High Energy Cosmic Rays during the Early Epochs of Radio-loud AGN

Powerful radio-loud active galactic nuclei (AGN) with large Mpc-scale jets have been theoretically motivated as emitters of high-energy cosmic rays. Recent radio observations have established a populous class of young radio-loud galaxies with compact ($< 1$ kpc) symmetric jets that are morphologically similar to large-scale AGNs. We show that these compact AGNs, so-called compact symmetric objects (CSOs), can accelerate protons up to $10^{20}$ eV at their hot spots via a Fermi type mechanism on the assumption of efficient acceleration. The required magnetic field strengths are comparable to those derived from the minimum energy condition. We further show that the accelerated protons can escape through the photon fields of the cocoon without significant energy loss. However, the local number density of powerful CSOs is insufficient for CSOs to power the entire observed flux of ultra-high-energy cosmic rays, providing maximally only a few percent. A heavy composition of UHECRs allows more CSOs to accelerate particles to UHECR energies, but escaping the cocoon is difficult. We comment on a method that may test CSOs as UHECR sources.

preprint2010arXiv

Revealing Type Ia supernova physics with cosmic rates and nuclear gamma rays

Type Ia supernovae (SNIa) remain mysterious despite their central importance in cosmology and their rapidly increasing discovery rate. The progenitors of SNIa can be probed by the delay time between progenitor birth and explosion as SNIa. The explosions and progenitors of SNIa can be probed by MeV nuclear gamma rays emitted in the decays of radioactive nickel and cobalt into iron. We compare the cosmic star formation and SNIa rates, finding that their different redshift evolution requires a large fraction of SNIa to have large delay times. A delay time distribution of the form t^{-1.0 +/- 0.3} provides a good fit, implying 50% of SNIa explode more than ~ 1 Gyr after progenitor birth. The extrapolation of the cosmic SNIa rate to z = 0 agrees with the rate we deduce from catalogs of local SNIa. We investigate prospects for gamma-ray telescopes to exploit the facts that escaping gamma rays directly reveal the power source of SNIa and uniquely provide tomography of the expanding ejecta. We find large improvements relative to earlier studies by Gehrels et al. in 1987 and Timmes & Woosley in 1997 due to larger and more certain SNIa rates and advances in gamma-ray detectors. The proposed Advanced Compton Telescope, with a narrow-line sensitivity ~ 60 times better than that of current satellites, would, on an annual basis, detect up to ~ 100 SNIa (3 sigma) and provide revolutionary model discrimination for SNIa within 20 Mpc, with gamma-ray light curves measured with ~ 10 sigma significance daily for ~ 100 days. Even more modest improvements in detector sensitivity would open a new and invaluable astronomy with frequent SNIa gamma-ray detections.