Hiroyasu Tajima

High-energy astrophysical neutrinos serve as crucial messengers for understanding hadronic acceleration processes and identifying the origins of cosmic rays, with blazars among the most promising neutrino sources. The KM3NeT experiment reported the detection of an ultrahigh-energy neutrino with an energy estimate of ∼220 PeV, the most energetic yet observed. The neutrino arrival direction has a 99% confidence region of 3° radius centered at R.A. 94 . ° 3, decl. −7 . ° 8 (J2000). In this work, 17 candidate blazars located within this region are identified. Comprehensive new observations and archival data analysis for these sources are presented. The study provides a complete multiwavelength coverage across radio, optical, X-ray, and gamma-ray bands, including proprietary data and dedicated follow-up observations. This study highlights flaring behavior in several candidate counterparts. One object exhibits a radio flare coinciding with the neutrino arrival time, with a pre-trial chance probability of 0.26%. Only for this object, the chance probability could be quantified in a robust way. Another candidate displays a rising trend in X-ray flux in a 1 yr window around the neutrino arrival time, while a third undergoes a gamma-ray flare during the same period. Based on the observational findings presented here, while none of these candidates can conclusively be linked to the neutrino, the implications of a possible blazar origin for the KM3NeT event are discussed.

We present very-high-energy gamma-ray observations of two binary black hole merger candidates, GW240615_113620 and GW241125_010116, performed with the Major Atmospheric Gamma Imaging Cherenkov (MAGIC) telescopes and the first Large-Sized Telescope of the Cherenkov Telescope Array Observatory’s (CTAO LST-1). GW240615_113620 was the best localized event of the fourth observing run of the LIGO-Virgo-KAGRA gravitational waves interferometers. GW241125_010116 was temporally and spatially coincident with a sub-threshold short-duration burst detected with the Swift-Burst Alert Telescope, the Swift-X-Ray Telescope and the Einstein Probe Follow-up X-ray Telescope. We observed the two events in stereoscopic mode, taking advantage of the improved sensitivity of joint MAGIC+LST-1 observations. No detection was achieved in the GeV–TeV gamma-ray band for any of the two sources. The unfavourable observing conditions of both events posed a challenge for a standard analysis and therefore a non standard analysis was necessary for both objects. Owing to the small localization area and the association with a GRB-like burst respectively, these events represented an unprecedented opportunity to study in details the electromagnetic emission from binary black holes merger events and, in particular, we discussed two theoretical models that predict a detectable gamma-ray emission and the possible future applications.

We investigate the nature of the unassociated sources detected by the Fermi-LAT close (∣ b ∣ < 10°) to the Galactic plane, representing 16% of all sources in the 4FGL-DR4 catalog. The bulk of these sources (referred to as soft Galactic unassociated sources, SGUs) exhibit properties not found in known classes of gamma-ray emitters, as confirmed by a machine-learning classification approach. In particular, these properties include a steep, curved spectrum peaking below 1 GeV and a specific Galactic-latitude distribution with both a narrow and a broad component (dubbed the spike and the shoulder, respectively). Some source clusters are highlighted. New plausible source classes are explored, but only star-forming regions are found to account for a significant fraction (at most 10%) of the unassociated population. A thorough search for counterparts to the 175 brightest sources brings out a number of plausible counterparts but does not reveal clues about the nature of the whole population. We investigate the possibility that SGUs originate from mismodeled clumps of diffuse emission. Using Monte Carlo simulations, the SGU spectra can be reproduced in this scenario under an ad hoc condition concerning the clump spatial extension. The possible connection between the SGUs and gas not accounted for by the ¹² CO tracer is explored using the ¹³ CO MOPRA data but leads to inconclusive results. The origin of SGUs being related to diffuse emission remains plausible. However, a scenario whereby SGUs represent a new class of gamma-ray emitters cannot be fully excluded.

The Cherenkov Telescope Array Observatory (CTAO) is a ground-based observatory that will improve upon the sensitivities of the current generation of very-high-energy gamma-ray instruments. The Schwarzschild-Couder Telescope (SCT) is a dual-mirror candidate design for a CTAO Medium-Sized Telescope (MST). The prototype Schwarzschild-Couder Telescope (pSCT) was inaugurated in 2019 at Fred Lawrence Whipple Observatory (FLWO) in Arizona and observed significant gamma-ray emission from the Crab Nebula with a partially populated camera. The pSCT camera is currently being upgraded to fully instrument the focal plane with 11,328 silicon photomultiplier (SiPM) pixels split between 177 camera modules. In addition, the modules will feature upgraded electronics designed to reduce electronics crosstalk and noise. A module calibration procedure has been developed using a preproduction test module. Following this calibration procedure, performance testing shows that the upgrade module has low noise, minimal electronics crosstalk, and excellent charge resolution. After calibration and optimization, the 177 production modules will be installed in the pSCT camera for commissioning. This will be followed by observations of known VHE gamma-ray sources for camera performance validation.

Flat spectrum radio quasars (FSRQs) often remain undetected at very-high-energy (VHE, E > 100 GeV) γ rays due to their typically soft spectra and large distance. In December 2023, the Large-Sized Telescope prototype (LST-1) detected for the first time VHE γ-ray emission from the FSRQ OP 313 during an extraordinarily bright flare, becoming the furthest blazar ever observed at VHE with z=0.997. We aim to characterize the γ-ray emission of OP 313 during this flare, comparing it with its average emission state in order to understand the processes leading to this detection. Its remarkable distance also enables studies on the extragalactic gackground light (EBL), with the goal of evaluating the attenuation of VHE γ-ray photons. Thanks to an intensive multiwavelength campaign, we also study in detail the broadband emission of the source. We characterize the γ-ray emission during the flare in December 2023 and the low emission state observed in January 2024 thanks to the LST-1 and MAGIC data and quasi-simultaneous -LAT observations. Fermi This dataset also enables us to evaluate the EBL attenuation by systematically exploring the EBL intensity over the γ-ray spectrum. . Finally, we study the multiwavelength emission and interpret the broadband spectral energy distribution (SED) within blazar radiative models including the thermal contributions from the accretion disk, dusty torus and broad line region. We measured significant VHE γ-ray emission during the December 2023 flare at a level of 0.3 Crab units above 100 GeV . We also characterize the flare brightness in the high-energy (HE, E>100 MeV) γ-ray band, which was found to be a factor of 50 above the average emission seen by -LAT. The HE and VHE observations allow us to set constrains on the EBL density, resulting in an upper limit at 0.6 μm of łambda I_ Fermi łambda <8.74 nW m 2 $ sr^-1 with LST-1 data, and łambda I_ łambda <14.7 nW m^2 sr^-1 with MAGIC data, including systematic uncertainties. Finally, thanks to the extensive multiwavelength campaign organized, we are able to construct and model the broadband SED of OP 313 within the framework of a two-zone leptonic model where the γ-ray emission is produced via inverse-Compton scattering of the broad line region, accretion disk, and dusty torus photon fields. However, the dominant external photon field remains unknown, as several combinations are able to successfully explain the $γ-ray emission observed .

Extremely high-peaked BL Lac objects—also named extreme blazars —are among the most energetic and persistent extragalactic accelerators in the Universe, defined by a synchrotron emission peaking above 10 ¹⁷ Hz in X-rays. Such emission is then reprocessed and produces radiation extending deeply into very-high-energy (VHE; energy E > 100 GeV) gamma-rays. Observations in this energy band—optimally investigated by the Imaging Air-Shower Cherenkov telescopes—are crucial for probing the physical processes that drive their extreme behavior. This study extends our investigation of extreme blazars in the VHE gamma-ray range, providing a second new mini-catalog of sources observed by the MAGIC telescopes. We report on the monitoring of seven targets between 2017 and 2025, including four newly observed sources and three that have been part of long-term observation campaigns, for a total of approximately 338 hr of observations. The analysis of MAGIC data reveals two new VHE detections of extreme blazars, along with three additional sources showing hints of VHE emission. Joint observations of MAGIC and the first Large-Sized Telescope (LST-1) also confirmed a new VHE extreme blazar. Our results are complemented by simultaneous multiwavelength observations in other energy bands, including optical-UV, X-rays, and high-energy gamma-rays (100 MeV < E < 100 GeV). We confirm typical behavior of extreme blazars, such as a modest variability and a “harder-when-brighter” trend in X-rays across the sample. This new set increases the population of extreme blazars, and their broadband analysis confirms the physical properties of these extreme sources.

Context. The Galactic magnetar SGR 1935+2154 has exhibited prolific high-energy (HE) bursting activity in recent years. Aims. Investigating its potential tera-electronvolt counterpart could provide insights into the underlying mechanisms of magnetar emission and very high-energy (VHE) processes in extreme astrophysical environments. We aim to search for a possible tera-electronvolt counterpart to both its persistent and sub-second-scale burst emission. Methods. We analysed over 25 hour of observations from the Large-Sized Telescope prototype (LST−1) of the Cherenkov Telescope Array Observatory (CTAO) during periods of HE activity from SGR 1935+2154 in 2021 and 2022 to search for persistent emission. For bursting emission, we selected and analysed nine 0.1 s time windows centred around known short X-ray bursts, targeting potential sub-second-scale tera-electronvolt counterparts in a low-photon-statistics regime. Results. While no persistent or bursting emission was detected in our search, we establish upper limits for the tera-electronvolt emission of a short magnetar burst simultaneous to its soft gamma-ray flux. Specifically, for the brightest burst in our sample, the ratio between tera-electronvolt and X-ray flux is ≲10 ⁻³ . Conclusions. The non-detection of either persistent or bursting tera-electronvolt emission from SGR 1935+2154 suggests that if such components exist, they may occur under specific conditions not covered by our observations. This aligns with theoretical predictions of VHE components in magnetar-powered fast radio bursts and the detection of MeV–GeV emission in giant magnetar flares. These findings underscore the potential of magnetars, fast radio bursts, and other fast transients as promising candidates for future observations in the low-photon-statistics regime with Imaging Atmospheric Cherenkov Telescopes, particularly with the CTAO.

An increasing number of pulsar wind nebulae (PWNe) are being identified in the TeV band by ground-based Imaging Air Cherenkov Telescopes such that they constitute the dominant source class of Galactic TeV emitters. However, MeV–GeV PWN counterparts are still largely lacking. To date, only a dozen PWNe are identified by the Fermi–Large Area Telescope (LAT) in the MeV–GeV band. Most PWNe are located along the Galactic plane embedded within the prominent, diffuse Galactic γ-ray emission, which makes these sources difficult to disentangle from the bright diffuse background. We present a systematic search for γ-ray counterparts to known PWNe in the 300 MeV–2 TeV energy band using the Fermi–LAT. We target the locations of previously identified PWNe that lack detected Fermi–LAT pulsars to minimize associated pulsar contamination. The sample includes six previously identified Fermi PWNe and eight Fermi–LAT sources associated with PWNe. We report the analysis of 58 regions of interest and classify Fermi–LAT detected sources as either a likely PWN or a candidate PWN counterpart based on their morphological and spectral characteristics across the broadband spectrum. There are nine unidentified Fermi–LAT sources that we consider as likely PWN counterparts, which, if confirmed to be PWNe, would greatly increase the PWN population detected by the Fermi–LAT from 12 to 21. The remaining Fermi–LAT detected sources are considered weaker PWN candidates. A second approach in the systematic search for γ-ray emitting PWNe will involve studying the off-pulse phases of Fermi–LAT detected pulsars for the presence of an obscured PWN and will be reported in a subsequent paper.

GRB 221009A is the brightest gamma-ray burst (GRB) observed to date. Extensive observations of its afterglow emission across the electromagnetic spectrum were performed, providing the first strong evidence of a jet with a nontrivial angular structure in a long GRB. We carried out an extensive observation campaign in very-high-energy (VHE) gamma rays with the first Large-Sized Telescope of the future Cherenkov Telescope Array Observatory starting on 2022 October 10, about 1 day after the burst. A dedicated analysis of the GRB 221009A data is performed to account for the different moonlight conditions under which data were recorded. We find an excess of gamma-like events with a statistical significance of 4.1σ during the observations taken 1.33 days after the burst, followed by background-compatible results for the later days. The results are compared with various models of afterglows from structured jets that are consistent with the published multiwavelength data but entail significant quantitative and qualitative differences in the VHE emission after 1 day. We disfavor models that imply VHE flux at 1 day considerably above 10⁻¹¹ erg cm⁻² s⁻¹. Our late-time VHE observations can help disentangle the degeneracy among the models and provide valuable new insight into the structure of GRB jets.