Showing 1–50 of 268 results for author: Jiang, W

Anisotropic Particle Transport from a Pulsar Wind Nebula Revealed by Einstein Probe and LHAASO

Authors: Zhen Cao, F. Aharonian, Y. X. Bai, Y. W. Bao, D. Bastieri, X. J. Bi, Y. J. Bi, W. Bian, J. Blunier, A. V. Bukevich, C. M. Cai, W. Y. Cao, Zhe Cao, J. Chang, J. F. Chang, E. S. Chen, G. H. Chen, H. K. Chen, L. F. Chen, Liang Chen, Long Chen, M. J. Chen, M. L. Chen, Q. H. Chen, S. Chen , et al. (320 additional authors not shown)

Abstract: Pulsar wind nebulae (PWNe) are major cosmic ray accelerators, yet the mechanisms transporting high-energy particles into the interstellar medium remain elusive. Building on the LHAASO discovery of an ultra-high-energy (UHE) $γ$-ray source near the bow-shock PWN powered by the pulsar PSR J1740+1000, we present a joint Einstein Probe (EP) and LHAASO study of this system. EP observations reveal an ex… ▽ More Pulsar wind nebulae (PWNe) are major cosmic ray accelerators, yet the mechanisms transporting high-energy particles into the interstellar medium remain elusive. Building on the LHAASO discovery of an ultra-high-energy (UHE) $γ$-ray source near the bow-shock PWN powered by the pulsar PSR J1740+1000, we present a joint Einstein Probe (EP) and LHAASO study of this system. EP observations reveal an extended X-ray tail far exceeding the structure previously seen by XMM-Newton. Updated LHAASO observations show that the $γ$-ray emission is elongated, with its major axis aligned with the extended X-ray tail revealed by EP. This is the first detection of an X-ray pulsar tail associated with a spatially coincident extended UHE $γ$-ray emission. The X-ray and $γ$-ray spectrum can be well explained with a single population of relativistic electrons via synchrotron and inverse Compton radiation, respectively, removing the need for particle re-acceleration during propagation. The results unambiguously show that electrons/positrons above 100 TeV are escaping from the PWN. Instead of the immediate, isotropic diffusion into ambient interstellar medium that is typically assumed, these particles are transported anisotropically over at least $\sim$10 pc, either guided by the background magnetic field or carried by an advective outflow. △ Less

Submitted 7 August, 2026; originally announced August 2026.

Comments: Accepted by Science China Physics, Mechanics, and Astronomy. Main text: 9 pages, 4 figures, 1 table; Supplementary Materials: 7 pages, 2 figures, 4 tables

Distinct Jet Properties in the X-Ray-Obscured State of GRS 1915+105

Authors: Xi Yan, Lang Cui, Sergei Trushkin, Shuangjing Xu, Wu Jiang, Zhen Yan, Sándor Frey, Timur Mufakharov, Ruchika Dhaka

Abstract: GRS 1915+105 has remained in an X-ray-obscured state since its transition from a long-lasting unobscured state in 2019. We report on 6.7-GHz East Asia VLBI Network observations of GRS 1915+105 obtained during strong radio flares detected at 2.3--11.2 GHz with the RATAN-600 radio telescope in 2025. Our images reveal two contrasting jet morphologies. The first epoch, associated with a flare evolving… ▽ More GRS 1915+105 has remained in an X-ray-obscured state since its transition from a long-lasting unobscured state in 2019. We report on 6.7-GHz East Asia VLBI Network observations of GRS 1915+105 obtained during strong radio flares detected at 2.3--11.2 GHz with the RATAN-600 radio telescope in 2025. Our images reveal two contrasting jet morphologies. The first epoch, associated with a flare evolving from an optically thick to an optically thin spectrum, shows a bright radio core accompanied by an extended jet structure. By contrast, the second epoch, observed near the peak of another flare displaying optically thin emission at lower frequencies, is dominated by two bright, symmetric, well-separated jet blobs and shows no detectable radio core. If these jets exhibited the apparent superluminal motions commonly observed prior to 2019, measurable angular shifts would be expected over the five-hour observations. However, no significant jet motion is detected. Combined with our derived jet speed of $βΓ\lesssim 0.40$, these results suggest that the jets launched during the current obscured state are slower than the relativistic jets ($βΓ\gtrsim 1$) observed earlier during the unobscured state. Together with the recently reported large variations in jet orientation, our findings in GRS 1915+105 provide robust support for the emerging paradigm that X-ray binary jets launched in obscured and unobscured states likely exhibit distinct propagation properties. △ Less

Submitted 30 July, 2026; originally announced July 2026.

Comments: This is a revised manuscript submitted to ApJ

The Extended Ultrahigh-energy Gamma-Ray Emission in the Vicinity of PSR J2238+5903

Authors: Zhen Cao, F. Aharonian, Y. X. Bai, Y. W. Bao, D. Bastieri, X. J. Bi, Y. J. Bi, W. Bian, J. Blunier, A. V. Bukevich, C. M. Cai, W. Y. Cao, Zhe Cao, J. Chang, J. F. Chang, E. S. Chen, G. H. Chen, H. K. Chen, L. F. Chen, Liang Chen, Long Chen, M. J. Chen, M. L. Chen, Q. H. Chen, S. Chen , et al. (305 additional authors not shown)

Abstract: We present a comprehensive analysis of the recently discovered TeV gamma-ray source, LHAASO J2238+5900. Based on data collected from the LHAASO, our fitting results suggest that the source is significantly extended with an angular extension of 0.54° \pm 0.01° and is spatially coincident with the pulsar PSR J2238+5903. Its spectrum is characterized by a power-law with a cutoff at 41.0\pm 3.5 TeV. A… ▽ More We present a comprehensive analysis of the recently discovered TeV gamma-ray source, LHAASO J2238+5900. Based on data collected from the LHAASO, our fitting results suggest that the source is significantly extended with an angular extension of 0.54° \pm 0.01° and is spatially coincident with the pulsar PSR J2238+5903. Its spectrum is characterized by a power-law with a cutoff at 41.0\pm 3.5 TeV. Additionally, the source exhibits a significant signal of 7.9σabove 100 TeV, implying that it is a PeVatron candidate. While the gamma-ray emission is consistent with a pulsar wind nebula (PWN) scenario, the relatively large extension size also allows for a halo interpretation, potentially caused by electron-positron pairs escaping from the PWN. △ Less

Submitted 23 July, 2026; originally announced July 2026.

Dynamics and geometry of the inner sub-parsec-scale jet in 3C 279 observed with the Event Horizon Telescope

Authors: Hendrik Mueller, Sebastiano D. von Fellenberg, Ai-Ling Zeng, Paul Tiede, Thomas P. Krichbaum, Roman Gold, Tuomas Savolainen, Jae-Young Kim, Sijia Peng, Teresa Toscano, Michael Janssen, Boris Georgiev, Dhanya G. Nair, Iniyan Natarajan, Lindy Blackburn, Kazunori Akiyama, Ezequiel Albentosa-Ruiz, Antxon Alberdi, Walter Alef, Juan Carlos Algaba, Rohan Ganesh Amanaganti, Richard Anantua, Eleni Antonopoulou, Keiichi Asada, Rebecca Azulay , et al. (253 additional authors not shown)

Abstract: The 2021 Event Horizon Telescope observations resolve the innermost jet region of the blazar 3C279 with unprecedented detail. The reconstructed images consistently reveal a compact core elongated nearly orthogonal to the large-scale jet axis. This rarely observed morphology recurs across multiple epochs and from 22-230 GHz and is therefore intrinsic rather than an imaging artifact. Geometric model… ▽ More The 2021 Event Horizon Telescope observations resolve the innermost jet region of the blazar 3C279 with unprecedented detail. The reconstructed images consistently reveal a compact core elongated nearly orthogonal to the large-scale jet axis. This rarely observed morphology recurs across multiple epochs and from 22-230 GHz and is therefore intrinsic rather than an imaging artifact. Geometric model fitting identifies several components with apparent speeds up to 10c, requiring bulk Lorentz factors greater than 10.3 and constraining viewing angles to extremely small values (smaller than one degree). Rest-frame brightness temperatures are systematically low (between 10^9 and 10^10 K), consistent with optically thin emission at 230 GHz. These results suggest that the jet bends toward the observer on sub-parsec scales, producing strong relativistic beaming. Possible drivers of the observed jet bending and temporal evolution include the jet's interaction with the interstellar medium, kink or Kelvin--Helmholtz instabilities, magnetic reconnection near the horizon, or binary-induced precession. However, the current temporal coverage of VLBI data remains insufficient to distinguish between these mechanisms. Continued multifrequency VLBI monitoring will be essential to constraining the dynamics and geometry of the jet base in 3C279. △ Less

Submitted 17 July, 2026; originally announced July 2026.

Comments: accepted for publication by A&A

Low-Frequency VLBI Network Using SKA-LOW

Authors: Hideyuki Kobayashi, Tomoaki Oyama, Kazuma Fujita, Syunsaku Suzuki, Shintaro Yoshiura, Hiroaki Misawa, Fuminori Tsuchiya, Kazuhiro Takefuji, Kenta Fujisawa, Kohtaro Niinuma, Keitaro Takahashi, Yashwant Gupta, Bhal Chandra Joshi, Vishwesh Marthi, M. A. Krishnakumar, Zhiqiang Shen, Wu Jiang, Yihua Yan, Wei Wang, Linjie Chen

Abstract: We propose the development of a low-frequency Very Long Baseline Interferometry (VLBI) network operating in the 100-350 MHz range, incorporating the Square Kilometre Array Low (SKA-LOW). SKA-LOW is expected to achieve exceptionally high sensitivity within this frequency band. By integrating SKA-LOW with other high-sensitivity radio telescopes located across the Asia-Pacific region, the proposed ne… ▽ More We propose the development of a low-frequency Very Long Baseline Interferometry (VLBI) network operating in the 100-350 MHz range, incorporating the Square Kilometre Array Low (SKA-LOW). SKA-LOW is expected to achieve exceptionally high sensitivity within this frequency band. By integrating SKA-LOW with other high-sensitivity radio telescopes located across the Asia-Pacific region, the proposed network is anticipated to deliver up to two orders of magnitude improvement in sensitivity compared to the existing VLBA. While several scientific themes utilizing low-frequency VLBI have already been proposed, we specifically advocate for astrometric studies employing existing VLBI stations to demonstrate the feasibility and scientific potential of this frequency regime. Furthermore, the combination of SKA-LOW with additional radio telescopes will enable high-fidelity imaging observations, significantly enhancing the quality and scope of low-frequency VLBI science. △ Less

Submitted 16 July, 2026; originally announced July 2026.

Comments: Published in Advancing Astrophysics with the SKA II (AASKAII), 2026 (arXiv:2606.20366). Report-no:AASKAII/Kobayashi01

Report number: AASKAII/Kobayashi01

Research Progress on Solar Small-Scale Dynamo

Authors: Wen-Jie Jiang, Lei Ni, Chun-Lan Jin, Zhi Xu, Mei Zhang

Abstract: The small-scale solar dynamo theory, as the core mechanism explaining the origin of the persistent, disperse weak magnetic field in the quiet Sun regions, has made significant progress over the past three decades in the fields of observation, theory, and simulation. Breakthrough observations from high-resolution space-based and ground-based telescopes have revealed that the quiet Sun is ubiquitous… ▽ More The small-scale solar dynamo theory, as the core mechanism explaining the origin of the persistent, disperse weak magnetic field in the quiet Sun regions, has made significant progress over the past three decades in the fields of observation, theory, and simulation. Breakthrough observations from high-resolution space-based and ground-based telescopes have revealed that the quiet Sun is ubiquitously populated by highly dynamic, mixed-polarity and possibly predominantly horizontal magnetic structures with complex topology. These observations confirm that their total magnetic flux is substantial, with a high and widely distributed magnetic flux emergence rate, strongly suggesting a local dynamo effect independent of the solar cycle and driven by intense turbulence and convection. Theoretical studies indicate that even in the challenging low magnetic Prandtl number ( Pm << 1) environment of the solar photosphere and convection zone, turbulent motions can self-excitedly convert kinetic energy into magnetic energy through the stretching, folding, and twisting of magnetic field lines. MHD simulations have successfully reproduced observed features, demonstrating that a pure small-scale dynamo can operate efficiently and sustain magnetic fields even in an open, stratified solar environment incorporating realistic physical processes. Current research strongly suggests that small-scale magnetic fields constitute the majority of the magnetic energy in the quiet Sun and also influence coronal heating, solar wind acceleration, and radiation distribution. This article conducts a literature review centered on observations, theoretical models, and numerical simulations of the small-scale dynamo, organizing and discussing the relevant research history and progress. Finally, it summarizes the content and provides an outlook on future research from multiple perspectives. △ Less

Submitted 2 July, 2026; originally announced July 2026.

Broadband multiwavelength properties of the archetypal blazar 3C 279 during the 2017 Event Horizon Telescope campaign

Authors: G. Principe, J. C. Algaba, E. Aviano, W. Y. Cheong, K. Hada, D. Haggard, A. Hahn, S. G. Jorstad, E. V. Kravchenko, Y. Kovalev, S. S. Lee, M. Lisakov, S. Markoff, A. P. Marscher, M. Sasada, P. Voitsik, Kazunori Akiyama, Ezequiel Albentosa-Ruiz, Antxon Alberdi, Walter Alef, Richard Anantua, Eleni Antonopoulou, Keiichi Asada, Rebecca Azulay, Anne-Kathrin Baczko , et al. (508 additional authors not shown)

Abstract: The archetypal blazar 3C 279 hosts a prominent relativistic jet and exhibits strong broadband variability across the electromagnetic spectrum. In April 2017, the Event Horizon Telescope (EHT) observed 3C 279, alongside one of the most extensive quasi-simultaneous multiwavelength (MWL) campaigns ever conducted. With the aim of investigating the physical processes governing 3C 279, we analyzed indiv… ▽ More The archetypal blazar 3C 279 hosts a prominent relativistic jet and exhibits strong broadband variability across the electromagnetic spectrum. In April 2017, the Event Horizon Telescope (EHT) observed 3C 279, alongside one of the most extensive quasi-simultaneous multiwavelength (MWL) campaigns ever conducted. With the aim of investigating the physical processes governing 3C 279, we analyzed individual observations and multiband light curves, and constructed a new quasi-simultaneous MWL spectrum. We also performed phenomenological modeling using the turbulent extreme multi-zone (TEMZ) model to constrain the fundamental physical properties of the source. The EHT observations reveal a clear flux increase in the innermost core between April 5 and 11, 2017. Over a broader timescale, radio measurements at longer wavelengths show concurrent enhancements in core flux and polarization around mid-April, coinciding with the ejection of a superluminal knot. Record UV-optical flares with strong polarization variability occurred in late March, followed by gamma-ray activity that declined before the end of the EHT observing period. During this interval, the source remained in a low X-ray state and showed no detectable VHE emission. The TEMZ modeling suggests that the broadband spectrum and variability of 3C 279 can be explained within a jet scenario in which turbulent plasma cells are compressed by a stationary conical shock. However, alternative interpretations, such as magnetic reconnection or a moving shock-in-jet event, remain plausible. This coordinated MWL campaign advances our understanding of the origin of jet and gamma-ray emission in 3C 279, while also providing a comprehensive publicly available dataset that will serve as a valuable reference for future studies. △ Less

Submitted 24 June, 2026; originally announced June 2026.

Comments: 29 pages, 13 figures, accepted by Astronomy & Astrophysics on June 18, 2026

Extreme PeV accelerator associated with GRS 1915+105

Authors: Zhen Cao, F. Aharonian, Y. X. Bai, Y. W. Bao, D. Bastieri, X. J. Bi, Y. J. Bi, W. Bian, J. Blunier, A. V. Bukevich, C. M. Cai, Y. Y. Cai, W. Y. Cao, Zhe Cao, J. Chang, J. F. Chang, E. S. Chen, G. H. Chen, H. K. Chen, L. F. Chen, Liang Chen, Long Chen, M. J. Chen, M. L. Chen, Q. H. Chen , et al. (304 additional authors not shown)

Abstract: Microquasars, binary systems featuring relativistic jets, have emerged as sources for particle acceleration beyond PeV energies. We present a study of the broadband $γ$-ray emission from one of the most prominent Galactic microquasars GRS 1915+105 based on data accumulated by LHAASO and Fermi-LAT over 4 and 17 years, respectively. A joint analysis of LHAASO-WCDA and LHAASO-KM2A data reveals extend… ▽ More Microquasars, binary systems featuring relativistic jets, have emerged as sources for particle acceleration beyond PeV energies. We present a study of the broadband $γ$-ray emission from one of the most prominent Galactic microquasars GRS 1915+105 based on data accumulated by LHAASO and Fermi-LAT over 4 and 17 years, respectively. A joint analysis of LHAASO-WCDA and LHAASO-KM2A data reveals extended $γ$-ray emission whose centroid appears significantly shifted, by ~ 0.13°, from the binary system and its jets. The spectral energy distribution is well described by a curved spectrum with progressive steepening that can be described by a log-parabola function with no evidence for a sharp cutoff, consistent with parent particles reaching multi-PeV energies and an extreme acceleration efficiency approaching the limit set by the available potential drop across the source. Several features, most notably the shift of the emission and single-power-law spectrum down to GeV band, favor radiation by cosmic rays accelerated in the source interacting with the dense ambient medium. Our spectral modeling implies that at least a few percent of the jet mechanical power is transferred to protons, whose maximum energy reaches beyond 5 PeV. These results strengthen the case for microquasars as exceptionally efficient accelerators in our Galaxy. △ Less

Submitted 25 June, 2026; v1 submitted 23 June, 2026; originally announced June 2026.

Comments: 10 pages, 4 figures, with supplementary material. Corrected a typo in the y-axis units of Fig. 2

Discovery of Unusual Jet Orientation Variations in the Microquasar GRS 1915+105

Authors: Xi Yan, Lang Cui, Wu Jiang, Zhen Yan, Sándor Frey, Sergei Trushkin, Timur Mufakharov, Ruchika Dhaka, Shuangjing Xu

Abstract: We report large day-timescale variations in the orientation of the southeast--northwest jet in the prototype microquasar GRS 1915+105. These results are based on three-epoch East Asia VLBI Network (EAVN) observations at 6.7 GHz, obtained during giant radio flares in 2025 detected by the RATAN-600 monitoring program. Our observations reveal the smallest position angle (PA) of… ▽ More We report large day-timescale variations in the orientation of the southeast--northwest jet in the prototype microquasar GRS 1915+105. These results are based on three-epoch East Asia VLBI Network (EAVN) observations at 6.7 GHz, obtained during giant radio flares in 2025 detected by the RATAN-600 monitoring program. Our observations reveal the smallest position angle (PA) of $118^\circ \pm 7^\circ$ ever measured for the jet in GRS 1915+105, which increases to $152^\circ \pm 2^\circ$ within 37 days. Based on the literature results, we further suggest that the jet orientation has exhibited significant variations over a PA range of $118^\circ$--$188^\circ$ since 2023. This unusual jet orientation behavior in GRS 1915+105 during its current X-ray-obscured state may arise from a warped, precessing inner accretion disk, as implied by recent X-ray spectroscopy. Notably, one image reveals a peculiar morphology in GRS 1915+105, which likely indicates lateral spreading of the approaching southeast jet. Future observations are essential to clarify the issues raised in this work. △ Less

Submitted 15 June, 2026; originally announced June 2026.

Comments: ApJL in press

Satellite gravity constraints on inner core viscosity and LLVPs density anomalies

Authors: Yachong An, Hao Ding, Fred D. Richards, Weiping Jiang, Jiancheng Li, Wenbin Shen

Abstract: Constraining the physical properties of Earth's deep interior, particularly the viscosity of the solid inner core and the density structure of large low-velocity provinces (LLVPs), remains a major challenge in geophysics. Here we develop a unified dynamical framework that combines mantle-inner core gravitational coupling (MICG) with torsional oscillations in the fluid outer core and show that thei… ▽ More Constraining the physical properties of Earth's deep interior, particularly the viscosity of the solid inner core and the density structure of large low-velocity provinces (LLVPs), remains a major challenge in geophysics. Here we develop a unified dynamical framework that combines mantle-inner core gravitational coupling (MICG) with torsional oscillations in the fluid outer core and show that their interaction can produce a distinct and testable geodetic signature. Guided by this prediction, we analyze satellite gravity observations together with independent corrections for surface mass variability. We identify a robust approximately 6-year signal in the Stokes coefficient Delta S22, while no corresponding stationary signal is detected in Delta C22. A signal with the same periodicity is independently detected in length-of-day variations (Delta LOD), and the two signals exhibit a near anti-phase relationship. Interpreting this coupled signature within the proposed framework allows us to constrain the inner core viscosity to approximately 4.6 (+/- 1.8) x 10^16 Pa s and the equatorial relief of the inner core boundary to a semi-axis difference of about 200 +/- 70 m. The inversion further indicates mean density anomalies of +5.5 (+/- 0.6) per mil at the base of LLVPs. These results indicate that satellite gravimetry provides a direct observational window into deep-Earth dynamics and the physical properties of Earth's deep interior. △ Less

Submitted 8 June, 2026; originally announced June 2026.

Comments: 53 pages, 14 figures, 1 table

A 0.03 Hz Radio Quasi-periodic Oscillation During the 2025 Flare of GRS 1915+105

Authors: Ya Xing Li, Lei Liu, Wu Jiang, Zhen Yan, Bo Xia, Zhi Qiang Shen

Abstract: Our weekly-cadence radio monitoring campaign captured a bright flare in 2025 from the microquasar GRS 1915+105, observed simultaneously in the S- and X-bands (2.25 GHz and 8.42 GHz) with a short single baseline of two radio telescopes in Shanghai. Through high time resolution analysis, we detected a significant and short-lived quasi-periodic oscillation (QPO) at $\sim$0.03 Hz and its harmonic (… ▽ More Our weekly-cadence radio monitoring campaign captured a bright flare in 2025 from the microquasar GRS 1915+105, observed simultaneously in the S- and X-bands (2.25 GHz and 8.42 GHz) with a short single baseline of two radio telescopes in Shanghai. Through high time resolution analysis, we detected a significant and short-lived quasi-periodic oscillation (QPO) at $\sim$0.03 Hz and its harmonic ($\sim$0.06 Hz) in both radio bands of two consecutive observations on MJD 60765 ($>5.9 σ$) and MJD 60772 (2.8$σ$). Crucially, the QPO frequency is identical in both radio bands and matches oscillations detected in previous years. The recurrence and wavelength independence of the QPO frequency suggest an intrinsic characteristic timescale of the accretion-jet system. △ Less

Submitted 1 June, 2026; originally announced June 2026.

Comments: 10 pages, 3 figures. Accepted for publication in ApJ

Constraints on Ultralight Scalar and Dark Photon Dark Matter from PPTA-DR3 and EPTA-DR2

Authors: Xiao-Song Hu, Siyuan Chen, Kuo Liu, Xingjiang Zhu, Shi-Yi Zhao, Wu Jiang, John Antoniadis, N. D. Ramesh Bhat, Amodio Carleo, Shi Dai, Valentina Di Marco, Huanchen Hu, Wenhua Ling, Yang Liu, Saurav Mishra, Christopher J Russell, Ryan M. Shannon, Clemente Smarra, Jingbo Wang, Lin Wang, Andrew Zic

Abstract: The cold dark matter model successfully describes the Universe on large scales, yet faces challenges at sub-galactic scales. Ultralight dark matter (ULDM), with particle masses around $10^{-22} \mathrm{eV}$, offers a promising solution to these small-scale issues. Pulsar Timing Arrays (PTAs), designed to detect nanohertz gravitational waves, can also provide a sensitive probe for ULDM signals. In… ▽ More The cold dark matter model successfully describes the Universe on large scales, yet faces challenges at sub-galactic scales. Ultralight dark matter (ULDM), with particle masses around $10^{-22} \mathrm{eV}$, offers a promising solution to these small-scale issues. Pulsar Timing Arrays (PTAs), designed to detect nanohertz gravitational waves, can also provide a sensitive probe for ULDM signals. In this work, we perform a Bayesian search for ULDM using PTA data sets, focusing on two types of signals: the oscillatory gravitational potential from scalar ULDM and the fifth-force interaction mediated by dark photon dark matter (DPDM). We incorporate pulsar distances in the analysis to better model the ULDM density. No statistically significant evidence for ULDM has been found, therefore we place 95% confidence-level upper limits on the relevant parameters. For scalar ULDM, our analysis does not exclude the scenario in which ULDM constitutes all of dark matter. The constraints from PPTA-DR3 show significant improvements over the earlier PPTA-DR2 (2018 Preview) across most of the mass range, and are consistent with the recent uncorrelated limits from other PTAs. We also present for the first time the DPDM constraints using EPTA data. The obtained bounds on the DPDM from the EPTA-DR2 and PPTA-DR3 are comparable to existing constraints. △ Less

Submitted 3 May, 2026; originally announced May 2026.

Comments: Accepted for publication in Physical Review D, 14 pages, 6 figures

A magnetotelluric image of the Curnamona Province and the adjacent Delamerian Orogen margin: new insights into the crustal architecture

Authors: Wenping Jiang, Michael Doublier, Russell Korsch, Andy Clark, Malcolm Nicoll, Adrian Hitchman, Yanbo Cheng

Abstract: We have used new magnetotelluric data collected in the Curnamona Province and the adjacent part of the Delamerian Orogen margin to image electrical conductivity structures and to inform the understanding of the crustal architecture within the regional geological context. The preferred 3D resistivity model confirms, and resolves in greater detail, crustal-scale conductive features that have been ma… ▽ More We have used new magnetotelluric data collected in the Curnamona Province and the adjacent part of the Delamerian Orogen margin to image electrical conductivity structures and to inform the understanding of the crustal architecture within the regional geological context. The preferred 3D resistivity model confirms, and resolves in greater detail, crustal-scale conductive features that have been mapped by the long-period data collected at half-degree spacing as part of the Australian Lithospheric Architecture Magnetotelluric Project (AusLAMP), that is, the prominent Curnamona Province Conductor and the two Nackara Arc conductors. The new model reveals that the eastern Nackara Arc (ENAC) conductor continues as the Broken Hill Conductor (BHC) into the Curnamona Province. Regional geological considerations suggest that its formation is possibly linked to rifting/extension in the early Cambrian. Although we recognise that the east-west trending Wilcannia Conductor could be a possible continuation of the ENAC-BHC zone, integration with recently acquired deep seismic reflection data and evaluation of the geological setting lead us to suggest that they are not genetically linked. We suggest that the Wilcannia Conductor is younger and most likely is related to late Delamerian (~500 Ma) or Siluro-Devonian magmatism. Finally, these conductivity anomalies may represent large-scale trans-crustal structures that control the emplacement of low volume alkaline ultramafic magmas, and show a spatial relationship with certain mineral deposit types, suggesting a possible control on the distribution and formation of metallogenic provinces/belts in the region. This will be further investigated in future work. △ Less

Submitted 27 April, 2026; originally announced April 2026.

Comments: this paper is currently under consideration of Earth, Planets and Space

East Asian VLBI Network astrometry toward the star-forming region G040.96+02.48 in the Extreme Outer Galaxy

Authors: Xianjin Shen, Zehao Lin, Nobuyuki Sakai, Ye Xu, Shuaibo Bian, Yuanwei Wu, Yan Sun, Dejian Liu, Jingjing Li, Bo Zhang, Shuangjing Xu, Tomoaki Oyama, Chungsik Oh, Wu Jiang, Lang Cui, Pengfei Jiang, Guanghui Li, Mareki Honma, Se-Jin Oh, Zhi-Qiang Shen, Na Wang

Abstract: Accurate astrometric measurements for star-forming regions located on the far side of the Milky Way remain scarce. In this work, we present the astrometric results for a 22\,GHz water maser associated with star-forming region G040.96+02.48 located on the far side of the Milky Way, using the East Asian VLBI Network. The target water maser's proper motion was determined to be ($μ_α\cosδ, μ_δ$) = (… ▽ More Accurate astrometric measurements for star-forming regions located on the far side of the Milky Way remain scarce. In this work, we present the astrometric results for a 22\,GHz water maser associated with star-forming region G040.96+02.48 located on the far side of the Milky Way, using the East Asian VLBI Network. The target water maser's proper motion was determined to be ($μ_α\cosδ, μ_δ$) = ($-2.06_{-0.51}^{+0.53}$, $-2.95_{-0.44}^{+0.45}$)~mas~yr$^{-1}$. The derived three-dimensional kinematic distance to the star-forming region is 20.2$\pm$3.2\,kpc, placing it slightly outside the Outer Scutum$-$Centaurus Arm. The corresponding vertical height of 872$\pm$139\,pc indicates a significant warp of the outer Galactic disk, which is in good agreement with the latest precessing warp model. Moreover, the resulting peculiar motions reveal a complex kinematic pattern, characterized by a large outward radial velocity of $-32\pm$18\,km~s$^{-1}$. Our observations substantially expand the valuable sample of star-forming regions with accurate astrometric measurements in the Extreme Outer Galaxy. △ Less

Submitted 26 April, 2026; originally announced April 2026.

Comments: 12 pages, 6 figures, 3 tables, accepted for publication in AJ

Ultra-high-energy $γ$-ray imprints from PeV particles accelerated by supernova remnants

Authors: Zhen Cao, F. Aharonian, Y. X. Bai, Y. W. Bao, D. Bastieri, X. J. Bi, Y. J. Bi, W. Bian, J. Blunier, A. V. Bukevich, C. M. Cai, Y. Y. Cai, W. Y. Cao, Zhe Cao, J. Chang, J. F. Chang, E. S. Chen, G. H. Chen, H. K. Chen, L. F. Chen, Liang Chen, Long Chen, M. J. Chen, M. L. Chen, Q. H. Chen , et al. (303 additional authors not shown)

Abstract: The quest for the origin of cosmic ray (CRs) is a fundamental issue in astrophysics. Shocks of supernova remnants (SNRs) have been considered as the dominant contributors to Galactic CRs below the spectral knee near $\sim 3$ petaelectronvolt (PeV). Whether SNRs are efficient accelerators of particles beyond PeV energies has long been debated. Here we report observations of very-high-energy $γ$-ray… ▽ More The quest for the origin of cosmic ray (CRs) is a fundamental issue in astrophysics. Shocks of supernova remnants (SNRs) have been considered as the dominant contributors to Galactic CRs below the spectral knee near $\sim 3$ petaelectronvolt (PeV). Whether SNRs are efficient accelerators of particles beyond PeV energies has long been debated. Here we report observations of very-high-energy $γ$-ray emission up to hundreds of TeV from two middle age shell-type SNRs, G150.3$+$4.5 and $γ$-Cygni, with the Large High Altitude Air Shower Observatory (LHAASO). Two (or three) distinct morphological/spectral components with convex spectral shapes are observed in both sources, with the low-energy one being more extended than the high-energy one. %Although it is possible that these high-energy components may be driven by powerful pulsars, The likely association of the high-energy component with molecular clouds at similar distances, and the weakness/absence of pulsar wind nebulae (PWNe) inside these SNRs clearly indicate for the first time that the highest energy emission is produced by collision of hadronic CRs up to PeV energies with the clouds. These results are compatible with the classic model prediction that PeV particles accelerated near the end of the free expansion phase of SNR evolution can illuminate nearby molecular clouds (MCs) to produce strong $γ$-ray emission. △ Less

Submitted 24 April, 2026; originally announced April 2026.

Comments: 31 pages, 8 figures, 3 Tables

Full-polarization millimeter wavelength variability of Sagittarius A* during the 2018 EHT campaign

Authors: Ezequiel Albentosa-Ruiz, Jasmin E. Washington, Nicola Marchili, Iván Martí-Vidal, Ciriaco Goddi, Maciek Wielgus, Alejandro Mus, Angelo Ricarte, Daniel P. Marrone, León D. S. Salas, Yuhei Iwata, Douglas F. Carlos, Alexandra J. Tetarenko, Kotaro Moriyama, Vedant Dhruv, Kazunori Akiyama, Antxon Alberdi, Walter Alef, Juan Carlos Algaba, Richard Anantua, Keiichi Asada, Rebecca Azulay, Uwe Bach, Anne-Kathrin Baczko, David Ball , et al. (250 additional authors not shown)

Abstract: Sagittarius A* (Srg A*), the supermassive black hole at the center of the Milky Way, provides a unique laboratory to study accretion dynamics and plasma processes near the event horizon. We investigated the variability and polarization properties of Srg A* using ALMA observations during the 2018 Event Horizon Telescope campaign. We analyzed high-cadence full-polarization light curves from ALMA at… ▽ More Sagittarius A* (Srg A*), the supermassive black hole at the center of the Milky Way, provides a unique laboratory to study accretion dynamics and plasma processes near the event horizon. We investigated the variability and polarization properties of Srg A* using ALMA observations during the 2018 Event Horizon Telescope campaign. We analyzed high-cadence full-polarization light curves from ALMA at millimeter wavelengths, performed time-series analysis, and investigated the temporal behavior during an X-ray flare observed by Chandra on 2018 April 24. The variability characteristics are compared with expectations from standard accretion flow models. We find low variability in total intensity ($σ/μ< 10\%$), but significantly higher variability in linear and circular polarization (~ 30% and ~ 50%, respectively). A time-series analysis reveals red-noise variability, with power spectral densities between -2 and -3 across all Stokes parameters. Polarized intensity shows stable intra-day timescales, while total intensity exhibits more variable timescales, suggesting distinct emission regions, with polarization likely arising from a coherent structure. On April 24, a statistically significant inter-band delay in polarized intensity coincides with a near-simultaneous X-ray and millimeter peak that deviates from the typical delayed flare scenario. This event also features enhanced millimeter variability and coherent polarization loop evolution. The observed simultaneity challenges standard models of transient synchrotron emission with cooling delays, favoring instead a scenario of continuous energy injection in an optically thin region. Our results offer new constraints on the physical mechanisms driving variability in Srg A*, and provide key observational input for refining theoretical models of accretion and plasma behavior in the vicinity of supermassive black holes. △ Less

Submitted 11 April, 2026; originally announced April 2026.

Journal ref: Astronomy & Astrophysics 708 (2026) A179

The Broken Similarity: Sinking and Merging of Dark Matter Subhalos Across Hierarchical Levels

Authors: Wenkang Jiang, Jiaxin Han, Kun Xu, Victor J. Forouhar Moreno, Feihong He, Zhaozhou Li, Chunyan Jiang, Yipeng Jing, Xiaohu Yang

Abstract: We investigate hierarchical mergers among subhalos within a $Λ$CDM simulation using the HBT+ subhalo finder. Unlike previous methods, HBT+ tracks subhalo evolution across hierarchy levels, identifying the coalescence of subhalo cores in phase-space as a ''sinking" event. This coalescence marks a distinct stalled phase in orbital decay, providing a physically motivated and natural definition of a r… ▽ More We investigate hierarchical mergers among subhalos within a $Λ$CDM simulation using the HBT+ subhalo finder. Unlike previous methods, HBT+ tracks subhalo evolution across hierarchy levels, identifying the coalescence of subhalo cores in phase-space as a ''sinking" event. This coalescence marks a distinct stalled phase in orbital decay, providing a physically motivated and natural definition of a resolved merger. Our main findings include: 1) Over 90% of sinking events occur between adjacent subhalo levels, while cross-level pathways arise from tidal stripping, group infall, and numerical constraints. 2) Resolved mergers are predominantly major mergers (mass ratios > 1:10), while the occurrence of minor mergers decreases with the dynamical age of the host halo. 3) Although deep-level subhalos have low mass ratios relative to the host halo, their high mass ratios relative to direct parents significantly boost merger statistics. Consequently, the satellite-satellite merger rate can rival or exceed the central-satellite rate at lower mass thresholds. 4) Satellite-satellite mergers are spatially biased toward the outer regions of the host, suggesting that the central tidal field suppresses their orbital decay. 5) A bidirectional sinking detection recovers 32% more sinking events than the original algorithm, revealing that child-dispersion-driven mergers are dominated by tidal heating at the final stage of sinking, while parent-dispersion-driven and doubly-identified events proceed primarily via orbital decay. Altogether, these results reveal a complex landscape of hierarchical satellite mergers that deviate from the self-similarity of host halo mergers, due to additional physical processes including dynamical friction and the scale-dependent halo growth history. △ Less

Submitted 8 April, 2026; v1 submitted 7 April, 2026; originally announced April 2026.

Comments: 19 pages, 16 figures

A large misalignment between continuous jet and discrete ejecta in microquasar GRS 1915+105 during its obscured phase

Authors: Wu Jiang, Xi Yan, Zhen Yan, Ya-Ping Li, Lang Cui, Zhi-Qiang Shen

Abstract: We report a large misalignment between the continuous jet and the discrete ejecta in GRS 1915+105, detected in April 2023 with the East Asian VLBI network (EAVN). Two-sided ejecta are shown at 6.7 GHz images and central continuous jets are resolved at 43 GHz by EAVN quasi-simultaneously. While the continuous jet was aligned with the long-standing jet position angle (PA) of about 147 degree, the di… ▽ More We report a large misalignment between the continuous jet and the discrete ejecta in GRS 1915+105, detected in April 2023 with the East Asian VLBI network (EAVN). Two-sided ejecta are shown at 6.7 GHz images and central continuous jets are resolved at 43 GHz by EAVN quasi-simultaneously. While the continuous jet was aligned with the long-standing jet position angle (PA) of about 147 degree, the discrete ejecta appeared at a markedly different PA about 188 degree, with the lowest intrinsic velocity about 0.35 c ever reported. A similar misalignment of PA between discrete ejecta and continuous jet was independently detected in late September of 2023 during a consecutive flare event. The pronounced and recurrent angular deviations suggest a time-variable jet launching geometry, which, in conjunction with the observed X-ray obscuration, can be attributable to a warped accretion disk. Our result could offer new insight into the fundamental differences between continuous jet and discrete ejecta, and broadly provides a clue to understand the phenomena for transient black hole X-Ray binaries and changing-look active galactic nuclei during the X-ray obscured phase. △ Less

Submitted 31 March, 2026; originally announced April 2026.

Comments: 9 pages, published in ApJL

Search for Light-Mass Fractionally Charged Particles in Space with DAMPE Experiment

Authors: F. Alemanno, Q. An, P. Azzarello, F. C. T. Barbato, P. Bernardini, X. J. Bi, H. V. Boutin, I. Cagnoli, M. S. Cai, E. Casilli, J. Chang, D. Y. Chen, J. L. Chen, Z. F. Chen, Z. X. Chen, P. Coppin, M. Y. Cui, T. S. Cui, I. De Mitri, F. de Palma, A. Di Giovanni, T. K. Dong, Z. X. Dong, G. Donvito, J. L. Duan , et al. (123 additional authors not shown)

Abstract: Free Fractionally Charged Particles (FCPs) are predicted by some theories beyond or extended to the standard model. FCPs have been widely searched for by underground and space-based experiments based on the assumption of heavy lepton-like particles. However, there is a paucity of research focusing on light-mass FCPs (LFCPs) in the sub-MeV mass range. In this work, we report the LFCPs in primary hi… ▽ More Free Fractionally Charged Particles (FCPs) are predicted by some theories beyond or extended to the standard model. FCPs have been widely searched for by underground and space-based experiments based on the assumption of heavy lepton-like particles. However, there is a paucity of research focusing on light-mass FCPs (LFCPs) in the sub-MeV mass range. In this work, we report the LFCPs in primary high energy cosmic rays, based on observational data from the Dark Matter Particle Explorer (DAMPE) satellite. This study utilized ten years on-orbit data of DAMPE to search for LFCPs with a charge of $\frac{2}{3}~e$. No LFCP candidate was observed. Upper flux limit of LFCPs with a mass of 0.511 MeV$/c^{2}$ and a charge of $\frac{2}{3}~e$ is determined to be $\rm 5.0 \times 10^{-11}\,cm^{-2}sr^{-1}s^{-1}$ at the $\rm 90\%$ confidence level. △ Less

Submitted 23 February, 2026; originally announced February 2026.

Diverse properties of electron Forbush decreases revealed by the Dark Matter Particle Explorer

Authors: F. Alemanno, Q. An, P. Azzarello, F. C. T. Barbato, P. Bernardini, X. J. Bi, H. Boutin, I. Cagnoli, M. S. Cai, E. Casilli, J. Chang, D. Y. Chen, J. L. Chen, Z. F. Chen, Z. X. Chen, P. Coppin, M. Y. Cui, T. S. Cui, I. De Mitri, F. de Palma, A. Di Giovanni, T. K. Dong, Z. X. Dong, G. Donvito, J. L. Duan , et al. (125 additional authors not shown)

Abstract: The Forbush decrease (FD) of cosmic rays is an important probe of the interplanetary environment disturbed by solar activities. In this work, we study the properties of 8 FDs electrons (including positrons) between 2 GeV and 20 GeV from January, 2016 to March, 2024, with the Dark Matter Particle Explorer. The maximum decrease amplitudes of these events are about 30% - 15%, and the amplitudes reduc… ▽ More The Forbush decrease (FD) of cosmic rays is an important probe of the interplanetary environment disturbed by solar activities. In this work, we study the properties of 8 FDs electrons (including positrons) between 2 GeV and 20 GeV from January, 2016 to March, 2024, with the Dark Matter Particle Explorer. The maximum decrease amplitudes of these events are about 30% - 15%, and the amplitudes reduce with energy. The recovery time of these events shows diverse behaviors of their energy-dependence. Some of them show strong energy-dependence, while some have a nearly constant recovery time. It has been shown that such diverse behaviors could be related with the geometry of the disturbed regions of the interplanetary space by coronal mass ejections (CME), represented by the combined effect of the CME velocity, angular spread, and ejection direction. △ Less

Submitted 21 February, 2026; originally announced February 2026.

Comments: 14 pages, 12 figures. Accepted for publication in PRL

Journal ref: Phys. Rev. Lett. 136, 101004 (2026)

LHAASO observation of Mrk 421 during 2021 March - 2024 March: a comprehensive VHE catalog of multi-timescale outbursts and its time average behavior

Authors: The LHAASO Collaboration, Zhen Cao, F. Aharonian, Y. X. Bai, Y. W. Bao, D. Bastieri, X. J. Bi, Y. J. Bi, W. Bian, J. Blunier, A. V. Bukevich, C. M. Cai, Y. Y. Cai, W. Y. Cao, Zhe Cao, J. Chang, J. F. Chang, E. S. Chen, G. H. Chen, H. K. Chen, L. F. Chen, Liang Chen, Long Chen, M. J. Chen, M. L. Chen , et al. (303 additional authors not shown)

Abstract: The Large High Altitude Air Shower Observatory (LHAASO) monitors sources within its field of view for up to 7 hours daily, achieving a duty cycle exceeding 98% and an annual point-source sensitivity of 1.5% Crab Units (CU) in the very high energy (VHE) band. This unbiased sky-survey mode facilitates systematic monitoring and investigation of outburst phenomena. In this paper, we present results fr… ▽ More The Large High Altitude Air Shower Observatory (LHAASO) monitors sources within its field of view for up to 7 hours daily, achieving a duty cycle exceeding 98% and an annual point-source sensitivity of 1.5% Crab Units (CU) in the very high energy (VHE) band. This unbiased sky-survey mode facilitates systematic monitoring and investigation of outburst phenomena. In this paper, we present results from an unprecedented three-year monitoring campaign (March 2021--March 2024) of Mrk421 using LHAASO, spanning energies from 0.4 TeV to 20 TeV. We find that the blazar stayed in a quiescent state in 2021 and became active starting in 2022 with a total of 23 VHE outburst events identified, where the highest observed daily significance reaches $20\,σ$ with a flux equivalent to approximately 3.3~CU. LHAASO's continuous monitoring suggests the flaring occupancy of Mrk~421 to be around 14%. During long-term monitoring, multiwavelength (MWL) variability and correlation analyses are conducted using complementary data from Fermi-LAT, MAXI-GSC, Swift-XRT, and ZTF. A significant correlation ($>3\,σ$) is observed between X-ray and VHE bands with no detectable time lag, while the correlation between GeV and TeV bands is weaker. The flux distribution of the TeV emission during the quiescent state is different from that in the active state, implying the existence of two modes of energy dissipation in the blazar jet. Using simultaneous MWL data, we also analyzed both the long-term and outburst-period SEDs, and discussed the possible origin of the outburst events. △ Less

Submitted 13 February, 2026; originally announced February 2026.

Comments: 34 pages, 20 figures

Locating the missing large-scale emission in the jet of M87* with short EHT baselines

Authors: Boris Georgiev, Paul Tiede, Sebastiano D. von Fellenberg, Michael Janssen, Iniyan Natarajan, Lindy Blackburn, Jongho Park, Erandi Chavez, Andrew T. West, Kotaro Moriyama, Jun Yi Koay, Hendrik Müller, Dhanya G. Nair, Avery E. Broderick, Maciek Wielgus, Kazunori Akiyama, Ezequiel Albentosa-Ruíz, Antxon Alberdi, Walter Alef, Juan Carlos Algaba, Richard Anantua, Keiichi Asada, Rebecca Azulay, Uwe Bach, Anne-Kathrin Baczko , et al. (258 additional authors not shown)

Abstract: In Very-Long Baseline Interferometric arrays, nearly co-located stations probe the largest scales and typically cannot resolve the observed source. In the absence of large-scale structure, closure phases constructed with these stations are zero and, since they are independent of station-based errors, they can be used to probe data issues. Here, we show with an expansion about co-located stations,… ▽ More In Very-Long Baseline Interferometric arrays, nearly co-located stations probe the largest scales and typically cannot resolve the observed source. In the absence of large-scale structure, closure phases constructed with these stations are zero and, since they are independent of station-based errors, they can be used to probe data issues. Here, we show with an expansion about co-located stations, how these trivial closure phases become non-zero with brightness distribution on smaller scales than their short baseline would suggest. When applied to sources that are made up of a bright compact and large-scale diffuse component, the trivial closure phases directly measure the centroid relative to the compact source and higher-order image moments. We present a technique to measure these image moments with minimal model assumptions and validate it on synthetic Event Horizon Telescope (EHT) data. We then apply this technique to 2017 and 2018 EHT observations of M87* and find a weak preference for extended emission in the direction of the large-scale jet. We also apply it to 2021 EHT data and measure the source centroid about 1 mas northwest of the compact ring, consistent with the jet observed at lower frequencies. △ Less

Submitted 19 January, 2026; originally announced January 2026.

Comments: 15 pages, 8 figures, Submitted to Astronomy & Astrophysics

Journal ref: A&A 709, A26 (2026)

Transient Large-Scale Anisotropy in TeV Cosmic Rays due to an Interplanetary Coronal Mass Ejection

Authors: Zhen Cao, F. Aharonian, Y. X. Bai, Y. W. Bao, D. Bastieri, X. J. Bi, Y. J. Bi, W. Bian, A. V. Bukevich, C. M. Cai, W. Y. Cao, Zhe Cao, J. Chang, J. F. Chang, A. M. Chen, E. S. Chen, G. H. Chen, H. X. Chen, Liang Chen, Long Chen, M. J. Chen, M. L. Chen, Q. H. Chen, S. Chen, S. H. Chen , et al. (291 additional authors not shown)

Abstract: Large- or medium-scale cosmic ray anisotropy at TeV energies has not previously been confirmed to vary with time. Transient anisotropy changes have been observed below 150 GeV, especially near the passage of an interplanetary shock and coronal mass ejection containing a magnetic flux rope ejected by a solar storm, which can trigger a geomagnetic storm with practical consequences. In such events, c… ▽ More Large- or medium-scale cosmic ray anisotropy at TeV energies has not previously been confirmed to vary with time. Transient anisotropy changes have been observed below 150 GeV, especially near the passage of an interplanetary shock and coronal mass ejection containing a magnetic flux rope ejected by a solar storm, which can trigger a geomagnetic storm with practical consequences. In such events, cosmic rays provide remote sensing of the magnetic field properties. Here we report the observation of transient large-scale anisotropy in TeV cosmic ray ions using data from the Large High Altitude Air Shower Observatory (LHAASO). We analyze hourly skymaps of the transient cosmic ray intensity excess or deficit, the gradient of which indicates the direction and magnitude of transient large-scale anisotropy across the field of view. We observe enhanced anisotropy above typical hourly fluctuations with $>$5$σ$ significance during some hours of November 4, 2021, in separate data sets for four primary cosmic ray energy ranges of median energy from $E$=0.7 to 3.1 TeV. The gradient varies with energy as $E^γ$, where $γ\approx-0.5$. At a median energy $\leq$1.0 TeV, this gradient corresponds to a dipole anisotropy of at least 1\%, or possibly a weaker anisotropy of higher order. This new type of observation opens the opportunity to study interplanetary magnetic structures using air shower arrays around the world, complementing existing in situ and remote measurements of plasma properties. △ Less

Submitted 6 January, 2026; originally announced January 2026.

Comments: 14 pages, 11 figures. Submitted to Physical Review Letters, November, 2025

Journal ref: Phys. Rev. Lett., 136, 251002 (2026)

Ring Asymmetry and Spin in M87*

Authors: Vadim Bernshteyn, Nicholas S. Conroy, Michi Bauböck, Paul Tiede, Abhishek V. Joshi, Ben S. Prather, Charles F. Gammie, the Event Horizon Telescope Collaboration, :, Kazunori Akiyama, Ezequiel Albentosa-Ruíz, Antxon Alberdi, Walter Alef, Juan Carlos Algaba, Richard Anantua, Keiichi Asada, Rebecca Azulay, Anne-Kathrin Baczko, David Ball, Bidisha Bandyopadhyay, John Barrett, Bradford A. Benson, Dan Bintley, Lindy Blackburn, Raymond Blundell , et al. (241 additional authors not shown)

Abstract: Event Horizon Telescope (EHT) images of the supermassive black hole M87* depict an asymmetric ring of emission. General relativistic magnetohydrodynamic (GRMHD) models of M87* and its accretion disk predict that the amplitude and location of the ring's peak brightness asymmetry should fluctuate due to turbulence in the source plasma. We compare the observed distribution of brightness asymmetry amp… ▽ More Event Horizon Telescope (EHT) images of the supermassive black hole M87* depict an asymmetric ring of emission. General relativistic magnetohydrodynamic (GRMHD) models of M87* and its accretion disk predict that the amplitude and location of the ring's peak brightness asymmetry should fluctuate due to turbulence in the source plasma. We compare the observed distribution of brightness asymmetry amplitudes to the simulated distribution in GRMHD models, across varying black hole spin $a_{*}$. We show that, for strongly magnetized (MAD) models, three epochs of EHT data marginally disfavor $|a_{*}| \lesssim 0.2$. This is consistent with the Blandford-Znajek model for M87's jet, which predicts that M87* should have nonzero spin. We show quantitatively how future observations could improve spin constraints, and discuss how improved spin constraints could distinguish between differing jet-launching mechanisms and black hole growth scenarios. △ Less

Submitted 1 April, 2026; v1 submitted 1 January, 2026; originally announced January 2026.

Comments: Published by ApJ. 10 figures, 3 tables

Journal ref: ApJ 1000 231 (2026)

Observations of the Fermi bubbles and the Galactic center excess with the DArk Matter Particle Explorer

Authors: F. Alemanno, Q. An, P. Azzarello, F. C. T. Barbato, P. Bernardini, X. J. Bi, H. Boutin, I. Cagnoli, M. S. Cai, E. Casilli, J. Chang, D. Y. Chen, J. L. Chen, Z. F. Chen, Z. X. Chen, P. Coppin, M. Y. Cui, T. S. Cui, I. De Mitri, F. de Palma, A. Di Giovanni, T. K. Dong, Z. X. Dong, G. Donvito, J. L. Duan , et al. (123 additional authors not shown)

Abstract: The DArk Matter Particle Explorer (DAMPE) is a space-borne high-energy particle detector that surveys the $γ$-ray sky above$\sim 2~\rm GeV$ with a peak acceptance of $\sim 0.2~\rm m^2\,sr$. With the 102 months of data collected by DAMPE, we show that the Fermi bubbles are detected at a significance of $\sim 26σ$ and identify a GeV excess in the direction of Galactic center at $\sim 7 σ$ confidence… ▽ More The DArk Matter Particle Explorer (DAMPE) is a space-borne high-energy particle detector that surveys the $γ$-ray sky above$\sim 2~\rm GeV$ with a peak acceptance of $\sim 0.2~\rm m^2\,sr$. With the 102 months of data collected by DAMPE, we show that the Fermi bubbles are detected at a significance of $\sim 26σ$ and identify a GeV excess in the direction of Galactic center at $\sim 7 σ$ confidence. Both spectra and morphology are consistent with those observed by Fermi-LAT and the GeV excess component can be interpreted by the dark matter annihilation with a mass of $\sim 50$ GeV and a velocity-averaged cross section of $\sim 10^{-26}~{\rm cm^{3}~s^{-1}}$ for the $χχ\rightarrow b\bar{b}$ channel. Our results thus provide the first independent detection of these two intriguing diffuse gamma-ray sources besides Fermi-LAT. △ Less

Submitted 30 March, 2026; v1 submitted 29 December, 2025; originally announced December 2025.

Comments: 24 pages, 18 figures and 6 tables, accepted for publication in ApJS

Journal ref: ApJS 284:22 (2026)

Energy-Dependent Shifts of Medium-Scale Anisotropies in Very-High-Energy Cosmic Rays Observed by LHAASO-KM2A

Authors: The LHAASO collabration, Zhen Cao, F. Aharonian, Y. X. Bai, Y. W. Bao, D. Bastieri, X. J. Bi, Y. J. Bi, W. Bian, A. V. Bukevich, C. M. Cai, W. Y. Cao, Zhe Cao, J. Chang, J. F. Chang, A. M. Chen, E. S. Chen, G. H. Chen, H. X. Chen, Liang Chen, Long Chen, M. J. Chen, M. L. Chen, Q. H. Chen, S. Chen , et al. (292 additional authors not shown)

Abstract: Small deviations from isotropy in the arrival directions of Galactic cosmic rays serve as a unique probe of the local magnetic environment. In this Letter, we report observations of medium-scale anisotropies (MSA) at energies above 10 TeV using the LHAASO-KM2A array. Our analysis identifies four regions of excess and four regions of deficit, each spanning angular scales of approximately ten degree… ▽ More Small deviations from isotropy in the arrival directions of Galactic cosmic rays serve as a unique probe of the local magnetic environment. In this Letter, we report observations of medium-scale anisotropies (MSA) at energies above 10 TeV using the LHAASO-KM2A array. Our analysis identifies four regions of excess and four regions of deficit, each spanning angular scales of approximately ten degrees. Crucially, we detect significant energy-dependent shifts in the centroids of two excess regions: Region B and the newly identified Region $\mathrm{\widetilde{D}}$. We also characterize the energy evolution of the fractional relative intensity across both excess and deficit regions. These findings imply that the observed anisotropies are shaped by the specific realization of the local turbulent magnetic field within the cosmic ray scattering length. Such energy-dependent behaviors impose strict constraints on local turbulence models and cosmic ray propagation theories. △ Less

Submitted 1 July, 2026; v1 submitted 20 December, 2025; originally announced December 2025.

Cygnus X-3: A variable petaelectronvolt gamma-ray source

Authors: The LHAASO Collaboration, Zhen Cao, F. Aharonian, Y. X. Bai, Y. W. Bao, D. Bastieri, X. J. Bi, Y. J. Bi, W. Bian, J. Blunier, A. V. Bukevich, C. M. Cai, Y. Y. Cai, W. Y. Cao, Zhe Cao, J. Chang, J. F. Chang, E. S. Chen, G. H. Chen, H. K. Chen, L. F. Chen, Liang Chen, Long Chen, M. J. Chen, M. L. Chen , et al. (306 additional authors not shown)

Abstract: We report the discovery of variable $γ$-rays up to petaelectronvolt from Cygnus X-3, an iconic X-ray binary. The $γ$-ray signal was detected with a statistical significance of approximately 10 $σ$ by the Large High Altitude Air Shower Observatory (LHAASO). Its intrinsic spectral energy distribution (SED), extending from 0.06 to 3.7 PeV, shows a pronounced rise toward 1 PeV after accounting for abs… ▽ More We report the discovery of variable $γ$-rays up to petaelectronvolt from Cygnus X-3, an iconic X-ray binary. The $γ$-ray signal was detected with a statistical significance of approximately 10 $σ$ by the Large High Altitude Air Shower Observatory (LHAASO). Its intrinsic spectral energy distribution (SED), extending from 0.06 to 3.7 PeV, shows a pronounced rise toward 1 PeV after accounting for absorption by the cosmic microwave background radiation. We find variability on month-long timescales at a significance of $8.6 σ$, coinciding with a high state of the GeV gamma-ray flux detected by the Fermi-LAT. This,together with a 3.2$σ$ evidence for orbital modulation, suggests that the PeV $γ$-rays originate within, or in close proximity to, the binary system itself. The observed energy spectrum and temporal modulation can be naturally explained by $γ$-ray production through photomeson processes in the innermost region of the relativistic jet, where protons need to be accelerated to tens of PeV energies. △ Less

Submitted 12 April, 2026; v1 submitted 18 December, 2025; originally announced December 2025.

Comments: Submitted to NSR

Measurement of the cosmic ray nickel energy spectrum from 10 GeV/n to 2 TeV/n with the DAMPE

Authors: F. Alemanno, Q. An, P. Azzarello, F. C. T. Barbato, P. Bernardini, X. J. Bi, H. V. Boutin, I. Cagnoli, M. S. Cai, E. Casilli, J. Chang, D. Y. Chen, J. L. Chen, Z. F. Chen, Z. X. Chen, P. Coppin, M. Y. Cui, T. S. Cui, I. De Mitri, F. de Palma, A. Di Giovanni, T. K. Dong, Z. X. Dong, G. Donvito, J. L. Duan , et al. (123 additional authors not shown)

Abstract: Nickel, one of the most tightly bound nuclei alongside iron, is the most abundant heavy element beyond iron in cosmic rays. With DAMPE's excellent charge resolution and broad energy range, a high-precision energy spectrum provides valuable insights into the acceleration sources of heavy nuclei and their propagation through the interstellar medium. In this analysis, we report the direct measurement… ▽ More Nickel, one of the most tightly bound nuclei alongside iron, is the most abundant heavy element beyond iron in cosmic rays. With DAMPE's excellent charge resolution and broad energy range, a high-precision energy spectrum provides valuable insights into the acceleration sources of heavy nuclei and their propagation through the interstellar medium. In this analysis, we report the direct measurement of cosmic-ray nickel spectrum from 10 GeV/n to 2 TeV/n with nine years of flight data. The nickel spectrum is consistent with a single power law with spectral index -2.60 +/- 0.03 from 40 GeV/n to 1 TeV/n. This work provides an accurate measurement of differential flux of nickel with kinetic energy extending to TeV/n for the first time. △ Less

Submitted 1 August, 2026; v1 submitted 12 December, 2025; originally announced December 2025.

Probing jet base emission of M87* with the 2021 Event Horizon Telescope observations

Authors: Saurabh, Hendrik Müller, Sebastiano D. von Fellenberg, Paul Tiede, Michael Janssen, Lindy Blackburn, Avery E. Broderick, Erandi Chavez, Boris Georgiev, Thomas P. Krichbaum, Kotaro Moriyama, Dhanya G. Nair, Iniyan Natarajan, Jongho Park, Andrew Thomas West, Maciek Wielgus, Kazunori Akiyama, Ezequiel Albentosa-Ruíz, Antxon Alberdi, Walter Alef, Juan Carlos Algaba, Richard Anantua, Keiichi Asada, Rebecca Azulay, Uwe Bach , et al. (260 additional authors not shown)

Abstract: We investigate the presence and spatial characteristics of the jet base emission in M87* at 230 GHz, enabled by the enhanced uv coverage in the 2021 Event Horizon Telescope (EHT) observations. The addition of the 12-m Kitt Peak Telescope and NOEMA provides two key intermediate-length baselines to SMT and the IRAM 30-m, giving sensitivity to emission structures at scales of $\sim250~μ$as and… ▽ More We investigate the presence and spatial characteristics of the jet base emission in M87* at 230 GHz, enabled by the enhanced uv coverage in the 2021 Event Horizon Telescope (EHT) observations. The addition of the 12-m Kitt Peak Telescope and NOEMA provides two key intermediate-length baselines to SMT and the IRAM 30-m, giving sensitivity to emission structures at scales of $\sim250~μ$as and $\sim2500~μ$as (0.02 pc and 0.2 pc). Without these baselines, earlier EHT observations lacked the capability to constrain emission on large scales, where a "missing flux" of order $\sim1$ Jy is expected. To probe these scales, we analyzed closure phases, robust against station-based gain errors, and modeled the jet base emission using a simple Gaussian offset from the compact ring emission at separations $>100~μ$as. Our analysis reveals a Gaussian feature centered at ($Δ$RA $\approx320~μ$as, $Δ$Dec $\approx60~μ$as), a projected separation of $\approx5500$ AU, with a flux density of only $\sim60$ mJy, implying that most of the missing flux in previous studies must arise from larger scales. Brighter emission at these scales is ruled out, and the data do not favor more complex models. This component aligns with the inferred direction of the large-scale jet and is consistent with emission from the jet base. While our findings indicate detectable jet base emission at 230 GHz, coverage from only two intermediate baselines limits reconstruction of its morphology. We therefore treat the recovered Gaussian as an upper limit on the jet base flux density. Future EHT observations with expanded intermediate-baseline coverage will be essential to constrain the structure and nature of this component. △ Less

Submitted 1 December, 2025; originally announced December 2025.

Comments: 19 pages, 18 figures. Abstract shortened with respect to the manuscript. Accepted in Astronomy & Astrophysics

Journal ref: A&A, 706, A27 (2026)

Charge-dependent spectral softenings of primary cosmic-rays below the knee

Authors: DAMPE Collaboration, Francesca Alemanno, Qi An, Philipp Azzarello, Felicia-Carla-Tiziana Barbato, Paolo Bernardini, Xiao-Jun Bi, Hugo Valentin Boutin, Irene Cagnoli, Ming-Sheng Cai, Elisabetta Casilli, Jin Chang, Deng-Yi Chen, Jun-Ling Chen, Zhan-Fang Chen, Zi-Xuan Chen, Paul Coppin, Ming-Yang Cui, Tian-Shu Cui, Ivan De Mitri, Francesco de Palma, Adriano Di Giovanni, Tie-Kuang Dong, Zhen-Xing Dong, Giacinto Donvito , et al. (124 additional authors not shown)

Abstract: In most particle acceleration or propagation theories, the characteristic features of the cosmic ray spectra due to acceleration limits or propagation phase changes are charge dependent. Alternatively, the interaction scenario would expect mass dependent spectral features in general. The observational verification of which relation takes effect in nature is still lack due to the difficulty of meas… ▽ More In most particle acceleration or propagation theories, the characteristic features of the cosmic ray spectra due to acceleration limits or propagation phase changes are charge dependent. Alternatively, the interaction scenario would expect mass dependent spectral features in general. The observational verification of which relation takes effect in nature is still lack due to the difficulty of measuring the spectra of individual particles up to very high energies. Here we report direct measurements of the carbon, oxygen, and iron spectra from ~20 gigavolts to ~100 teravolts (~60 teravolts for iron) with 9 years of on-orbit data collected by the Dark Matter Particle Explorer. Distinct spectral softenings have been directly detected in these spectra for the first time. Combined with the updated proton and helium spectra, the spectral softening appears universally at a rigidity of ~15 teravolts. A nuclei mass dependent softening is rejected at a confidence level of >99.999%. Possible interpretations of these results, including a nearby cosmic ray source and other models such as the propagation effect, are discussed. △ Less

Submitted 30 April, 2026; v1 submitted 7 November, 2025; originally announced November 2025.

Comments: 36 pages, 8 figures, published online in Nature

LHAASO Detection of Ultra-High-Energy Gamma-Ray Emission toward the Giant Molecular Clouds

Authors: Zhen Cao, F. Aharonian, Y. X. Bai, Y. W. Bao, D. Bastieri, X. J. Bi, Y. J. Bi, W. Bian, A. V. Bukevich, C. M. Cai, W. Y. Cao, Zhe Cao, J. Chang, J. F. Chang, A. M. Chen, E. S. Chen, G. H. Chen, H. X. Chen, Liang Chen, Long Chen, M. J. Chen, M. L. Chen, Q. H. Chen, S. Chen, S. H. Chen , et al. (292 additional authors not shown)

Abstract: The $γ$-ray from Giant molecular clouds (GMCs) is regarded as the most ideal tool to perform in-situ measurement of cosmic ray (CR) density and spectra in our Galaxy. We report the first detection of $γ$-ray emissions in the very-high-energy (VHE) domain from the five nearby GMCs with a stacking analysis based on a 4.5-year $γ$-ray observation with the Large High Altitude Air Shower Observatory (L… ▽ More The $γ$-ray from Giant molecular clouds (GMCs) is regarded as the most ideal tool to perform in-situ measurement of cosmic ray (CR) density and spectra in our Galaxy. We report the first detection of $γ$-ray emissions in the very-high-energy (VHE) domain from the five nearby GMCs with a stacking analysis based on a 4.5-year $γ$-ray observation with the Large High Altitude Air Shower Observatory (LHAASO) experiment. The spectral energy distributions derived from the GMCs are consistent with the expected $γ$-ray flux produced via CR interacting with the ISM in the energy interval 1 - 100 $~\rm$ TeV. In addition, we investigate the presence of the CR spectral `knee' by introducing a spectral break in the $γ$-ray data. While no significant evidence for the CR knee is found, the current KM2A measurements from GMCs strongly favor a proton CR knee located above 0.9$~\rm$ PeV, which is consistent with the latest measurement of the CR spectrum by ground-based experiments. △ Less

Submitted 5 January, 2026; v1 submitted 7 November, 2025; originally announced November 2025.

Precise Measurement of the Cosmic Ray Helium Spectrum above 0.1 PeV

Authors: LHAASO Collaboration, Zhen Cao, F. Aharonian, Y. X. Bai, Y. W. Bao, D. Bastieri, X. J. Bi, Y. J. Bi, W. Bian, J. Blunier, A. V. Bukevich, C. M. Cai, Y. Y. Cai, W. Y. Cao, Zhe Cao, J. Chang, J. F. Chang, E. S. Chen, G. H. Chen, H. K. Chen, L. F. Chen, Liang Chen, Long Chen, M. J. Chen, M. L. Chen , et al. (303 additional authors not shown)

Abstract: We report a measurement of the cosmic ray helium energy spectrum in the energy interval 0.16 -- 13~PeV, derived by subtracting the proton spectrum from the light component~(proton and helium) spectrum obtained with observations made by the Large High Altitude Air Shower Observatory~(LHAASO) under a consistent energy scale. The helium spectrum shows a significant hardening centered at $E \simeq$ 1.… ▽ More We report a measurement of the cosmic ray helium energy spectrum in the energy interval 0.16 -- 13~PeV, derived by subtracting the proton spectrum from the light component~(proton and helium) spectrum obtained with observations made by the Large High Altitude Air Shower Observatory~(LHAASO) under a consistent energy scale. The helium spectrum shows a significant hardening centered at $E \simeq$ 1.1~PeV, followed by a softening at $\sim$ 7 PeV, indicating the appearance of a helium `knee'. Comparing the proton and helium spectra in the LHAASO energy range reveals some remarkable facts. In the lower part of this range, in contrast to the behavior at lower energies, the helium spectrum is significantly softer than the proton spectrum. This results in protons overtaking helium nuclei and becoming the largest cosmic ray component at $E \simeq$ 0.7 PeV. A second crossing of the two spectra is observed at $E \simeq$ 5 PeV, above the proton knee, when helium nuclei overtake protons to become the largest cosmic ray component again. These results have important implications for our understanding of the Galactic cosmic ray sources. △ Less

Submitted 1 April, 2026; v1 submitted 7 November, 2025; originally announced November 2025.

A semi-analytical mock galaxy catalog for the CSST extragalactic surveys from the Jiutian simulations

Authors: Zhenlin Tan, Lizhi Xie, Jiaxin Han, Yisheng Qiu, Fabio Fontanot, Gabriella De Lucia, Qi Guo, Qingyang Li, Jiale Zhou, Wenkang Jiang, Xin Wang, Feihong He, Chichuan Jin, Yipeng Jing, Ming Li, Xiaodong Li, Wenxiang Pei, Wenting Wang, Xiaohu Yang, Yu Yu

Abstract: We introduce a mock galaxy catalog built for the CSST extragalactic surveys using the primary runs of the Jiutian $N$-body simulation suites. The catalogs are built by coupling the GAlaxy Evolution and Assembly (GAEA) semi-analytical model of galaxy formation with merger trees extracted from the simulations using the Hierarchical Bound-Tracing (HBT+) algorithm. The spectral energy distributions (S… ▽ More We introduce a mock galaxy catalog built for the CSST extragalactic surveys using the primary runs of the Jiutian $N$-body simulation suites. The catalogs are built by coupling the GAlaxy Evolution and Assembly (GAEA) semi-analytical model of galaxy formation with merger trees extracted from the simulations using the Hierarchical Bound-Tracing (HBT+) algorithm. The spectral energy distributions (SEDs) and broadband magnitudes are computed using the neural-network-based stellar population synthesizer StarDuster, which is trained on radiative transfer simulations to account for detailed galaxy geometry in modeling dust obscuration. Galaxy light-cones up to $z=5$ are subsequently generated with the BLiC light-cone builder which interpolates the properties of galaxies over time using an optimized interpolation scheme. The resulting catalogs exhibit good convergence in many statistical properties of the galaxy population produced from two different resolution simulations. The catalogs reproduce a number of observed galaxy properties across a range of galaxy mass and redshift, including the stellar mass functions, the luminosity function, gas mass fraction, galaxy size-mass relation and galaxy clustering. We also present the photometric and redshift distributions of galaxies expected to be observed in the CSST surveys. △ Less

Submitted 9 November, 2025; v1 submitted 5 November, 2025; originally announced November 2025.

Comments: accepted by SCPMA, data products will be released at https://jiutian.sjtu.edu.cn

Evidence of cosmic-ray acceleration up to sub-PeV energies in the supernova remnant IC 443

Authors: Zhen Cao, F. Aharonian, Y. X. Bai, Y. W. Bao, D. Bastieri, X. J. Bi, Y. J. Bi, W. Bian, A. V. Bukevich, C. M. Cai, W. Y. Cao, Zhe Cao, J. Chang, J. F. Chang, A. M. Chen, E. S. Chen, G. H. Chen, H. X. Chen, Liang Chen, Long Chen, M. J. Chen, M. L. Chen, Q. H. Chen, S. Chen, S. H. Chen , et al. (291 additional authors not shown)

Abstract: Supernova remnants (SNRs) have been considered as the primary contributors to cosmic rays (CRs) in our Galaxy. However, the maximum energy of particles that can be accelerated by shocks of SNRs is uncertain, and SNRs' contribution to CRs around PeV energies is unclear. In this study, we present observations of high-energy $γ$-ray emission from the SNR IC 443 using the Large High Altitude Air Showe… ▽ More Supernova remnants (SNRs) have been considered as the primary contributors to cosmic rays (CRs) in our Galaxy. However, the maximum energy of particles that can be accelerated by shocks of SNRs is uncertain, and SNRs' contribution to CRs around PeV energies is unclear. In this study, we present observations of high-energy $γ$-ray emission from the SNR IC 443 using the Large High Altitude Air Shower Observatory (LHAASO). The morphological analysis reveals a pointlike source whose location and spectrum are consistent with those of the Fermi-LAT-detected compact source with $π^0$-decay signature, and a more extended source that is associated with a newly discovered Fermi source. The spectrum of the point source can be described by a power-law function with an index of $\sim3.0$, extending beyond $\sim 30$ TeV without apparent cutoff. Assuming a hadronic origin of the $γ$-ray emission, the $95\%$ lower limit of accelerated protons reaches about 300 TeV. The extended source might be associated with IC 443, SNR G189.6+3.3 or the putative pulsar wind nebula CXOU J061705.3+222127, and can be explained by either a hadronic or a leptonic model with particles reaching hundreds of TeV. These LHAASO results provide compelling evidence that sub-PeV CRs can be accelerated efficiently by shocks of SNRs. △ Less

Submitted 5 March, 2026; v1 submitted 29 October, 2025; originally announced October 2025.

First detection of ultra-high energy emission from gamma-ray binary LS I +61 303

Authors: Zhen Cao, F. Aharonian, Y. X. Bai, Y. W. Bao, D. Bastieri, X. J. Bi, Y. J. Bi, W. Bian, J. Blunier, A. V. Bukevich, C. M. Cai, Y. Y. Cai, W. Y. Cao, Zhe Cao, J. Chang, J. F. Chang, E. S. Chen, G. H. Chen, H. K. Chen, L. F. Chen, Liang Chen, Long Chen, M. J. Chen, M. L. Chen, Q. H. Chen , et al. (302 additional authors not shown)

Abstract: We report the first detection of gamma-ray emission up to ultra-high-energy (UHE; $>$100 TeV) emission from the prototypical gamma-ray binary system LS I +61 303 using data from the Large High Altitude Air Shower Observatory (LHAASO). It is detected with significances of 9.2$σ$ in WCDA (1.4--30.5 TeV) and 6.2$σ$ in KM2A (25--267 TeV); in KM2A alone we identify 16 photon-like events above 100 TeV a… ▽ More We report the first detection of gamma-ray emission up to ultra-high-energy (UHE; $>$100 TeV) emission from the prototypical gamma-ray binary system LS I +61 303 using data from the Large High Altitude Air Shower Observatory (LHAASO). It is detected with significances of 9.2$σ$ in WCDA (1.4--30.5 TeV) and 6.2$σ$ in KM2A (25--267 TeV); in KM2A alone we identify 16 photon-like events above 100 TeV against an estimated 5.1 background events, corresponding to a 3.8$σ$ detection. These results provide compelling evidence of extreme particle acceleration in LS I +61 303. Furthermore, we observe orbital modulation at 3.9$σ$ confidence level, between 25 and 100 TeV, with a hint that the orbital modulation is energy-dependent. These features can be understood in a composite scenario in which leptonic and hadronic processes jointly contribute. △ Less

Submitted 10 March, 2026; v1 submitted 27 October, 2025; originally announced October 2025.

Comments: 14 pages, 7 figures, 1 table, accepted for publication in PRL

Energy calibration of LHAASO-KM2A using the cosmic ray Moon shadow

Authors: Zhen Cao, F. Aharonian, Y. X. Bai, Y. W. Bao, D. Bastieri, X. J. Bi, Y. J. Bi, W. Bian, J. Blunier, A. V. Bukevich, C. M. Cai, Y. Y. Cai, W. Y. Cao, Zhe Cao, J. Chang, J. F. Chang, E. S. Chen, G. H. Chen, H. K. Chen, L. F. Chen, Liang Chen, Long Chen, M. J. Chen, M. L. Chen, Q. H. Chen , et al. (305 additional authors not shown)

Abstract: We present a precise measurement of the westward, rigidity-dependent shift of the Moon's shadow using three and a half years of cosmic-ray data collected by the Kilometer Square Array (KM2A) of the Large High Altitude Air Shower Observatory (LHAASO). These measurements enable us to calibrate the detector energy response in the range 20-260 TeV, with results showing excellent agreement with the res… ▽ More We present a precise measurement of the westward, rigidity-dependent shift of the Moon's shadow using three and a half years of cosmic-ray data collected by the Kilometer Square Array (KM2A) of the Large High Altitude Air Shower Observatory (LHAASO). These measurements enable us to calibrate the detector energy response in the range 20-260 TeV, with results showing excellent agreement with the response derived from Monte Carlo (MC) simulations of the KM2A detector. We also measure a best-fit parameter $ε= 0.015 \pm 0.08$, corresponding to a 95% confidence interval of [-14%, +17%] for the energy-scale estimation. This result establishes the exceptional accuracy of the KM2A-MC in simulating the detector's response within this energy range. △ Less

Submitted 7 January, 2026; v1 submitted 14 October, 2025; originally announced October 2025.

Comments: 12 pages, 13 figures, 2 tables

Journal ref: Physical Review D 112, 123051 (2025)

Towards the Giant Radio Array for Neutrino Detection (GRAND): the GRANDProto300 and GRAND@Auger prototypes

Authors: GRAND Collaboration, Jaime Álvarez-Muniz, Rafael Alves Batista, Aurélien Benoit-Lévy, Teresa Bister, Martina Bohacova, Mauricio Bustamante, Washington Carvalho, Yiren Chen, LingMei Cheng, Simon Chiche, Jean-Marc Colley, Pablo Correa, Nicoleta Cucu Laurenciu, Zigao Dai, Rogerio M. de Almeida, Beatriz de Errico, João R. T. de Mello Neto, Krijn D. de Vries, Valentin Decoene, Peter B. Denton, Bohao Duan, Kaikai Duan, Ralph Engel, William Erba , et al. (96 additional authors not shown)

Abstract: The Giant Radio Array for Neutrino Detection (GRAND) is a proposed multi-messenger observatory of Ultra-High-Energy (UHE) particles of cosmic origin. Its main goal is to find the long-sought origin of UHE cosmic rays by detecting large numbers of them and the secondary particles created by their interactions like gamma rays and neutrinos. The GRAND Collaboration plans to achieve this using large a… ▽ More The Giant Radio Array for Neutrino Detection (GRAND) is a proposed multi-messenger observatory of Ultra-High-Energy (UHE) particles of cosmic origin. Its main goal is to find the long-sought origin of UHE cosmic rays by detecting large numbers of them and the secondary particles created by their interactions like gamma rays and neutrinos. The GRAND Collaboration plans to achieve this using large arrays of radio antennas that look for the radio signals emitted by the air showers initiated by the interactions of the UHE particles in the atmosphere. Since 2023, three small-scale prototype GRAND arrays have been in operation: GRAND@Nançay in France, GRAND@Auger in Argentina, and GRANDProto300 in China. Together, their goal is to validate the detection principle of GRAND under prolonged field conditions, achieving efficient, autonomous radio-detection of air showers. We describe the hardware, software, layout, and operation of the GRAND prototypes. Using their data, we show a first characterization of the local electromagnetic environment of each site and a measurement of the Galactic synchrotron emission. Despite challenges, the successful operation of the prototypes confirms that the GRAND instrumentation is apt to address the goals of the experiment and lays the groundwork for its ensuing stages. △ Less

Submitted 24 February, 2026; v1 submitted 25 September, 2025; originally announced September 2025.

Comments: 30 pages, 16 figures, 2 tables

Journal ref: IOP Publishing, 21 (2026) 02

Investigation of hadronic cross sections of cosmic ray carbon and oxygen on BGO from 200 GeV to 10 TeV energy at the DAMPE experiment

Authors: F. Alemanno, Q. An, P. Azzarello, F. C. T. Barbato, P. Bernardini, X. J. Bi, H. Boutin, I. Cagnoli, M. S. Cai, E. Casilli, E. Catanzani, J. Chang, D. Y. Chen, J. L. Chen, Z. F. Chen, Z. X. Chen, P. Coppin, M. Y. Cui, T. S. Cui, Y. X. Cui, I. De Mitri, F. de Palma, A. Di Giovanni, T. K. Dong, Z. X. Dong , et al. (122 additional authors not shown)

Abstract: The Dark Matter Particle Explorer (DAMPE) has made significant progress in measuring the fluxes of cosmic rays. These new measurements are pivotal in advancing our understanding of the origins and propagation mechanisms of cosmic rays. The bismuth germanium oxide (BGO) calorimeter plays a crucial role in these measurements, particularly in the precise determination of cosmic ray fluxes. However, f… ▽ More The Dark Matter Particle Explorer (DAMPE) has made significant progress in measuring the fluxes of cosmic rays. These new measurements are pivotal in advancing our understanding of the origins and propagation mechanisms of cosmic rays. The bismuth germanium oxide (BGO) calorimeter plays a crucial role in these measurements, particularly in the precise determination of cosmic ray fluxes. However, for a calorimetric experiment like DAMPE, uncertainties in hadronic models persist as a major barrier in achieving more accurate measurements of fluxes of cosmic ray nuclei. This study centers on the measurement of the inelastic hadronic cross sections of carbon and oxygen nuclei interacting with BGO crystals target over an extensive energy range, spanning from 200 GeV to 10 TeV. For carbon nuclei interacting with the BGO target, the measurements of the cross sections have achieved a total relative uncertainty of less than 10% below 8 TeV for carbon, and below 3 TeV for oxygen. For oxygen nuclei, the same level of precision was attained below 3 TeV. Additionally, we compare the experimental results with Geant4 and FLUKA simulations to validate the accuracy and consistency of these simulation tools. Through comprehensive analysis of the inelastic hadronic interaction cross sections, this research provides validation for the hadronic interaction models used in DAMPE's cosmic-ray flux measurements. △ Less

Submitted 21 September, 2025; originally announced September 2025.

Millisecond Pulsars in M2: New discoveries and a detailed timing analysis

Authors: Baoda Li, Kuo Liu, Lin Wang, Shunyi Lan, P. C. C. Freire, Pinsong Zhao, Liyun Zhang, Zhengwei Liu, Lei Qian, Wu Jiang, Dejiang Yin, Yaowei Li, Yinfeng Dai, Yang Liu, Xiangcun Meng, Zhichen Pan

Abstract: Globular clusters (GCs) offer a unique environment for discovering and studying millisecond pulsars. In this paper, we present a multi-epoch search and detailed timing analysis of millisecond pulsars in the GC M2, using the Five-hundred-meter Aperture Spherical Telescope. We have discovered two new binary millisecond pulsars in M2, designated M2F and M2G, respectively. We provide measurements of t… ▽ More Globular clusters (GCs) offer a unique environment for discovering and studying millisecond pulsars. In this paper, we present a multi-epoch search and detailed timing analysis of millisecond pulsars in the GC M2, using the Five-hundred-meter Aperture Spherical Telescope. We have discovered two new binary millisecond pulsars in M2, designated M2F and M2G, respectively. We provide measurements of the emission properties of all known pulsars in M2, including their polarization profiles, rotation measures, flux densities, scintillation characteristics, and so forth. In particular, we report the first rotation measure at the distance and direction of this cluster. Additionally, we report the first phase-coherent timing solutions for the M2 pulsars. From our Bayesian timing analysis, we have measured their spin and orbital parameters with high precision, including the advance of periastron for M2A and M2E indicating total system masses of 1.75(13) and 1.80(5) solar masses respectively. Using archival data from the Hubble Space Telescope, we have identified an optical counterpart of M2C, which is likely the white dwarf companion of the pulsar. By combining results from optical and radio observations, we have reconstructed the binary evolution track of this system and estimated the cooling age of the companion to be approximately 10\,Myr, making it the youngest white dwarf in any known GC binary pulsars. Furthermore, using the spin period derivatives of M2 pulsars, we have investigated the gravitational potential of the cluster and found that our results strongly support the latest central-stellar-velocity dispersion measurement in M2. △ Less

Submitted 8 August, 2025; originally announced August 2025.

Comments: 12 figures, 6 tables, accepted for publication in ApJ

Jet Collimation Profile of Low-Luminosity AGN M84: Insight into the Jet Formation in the Low Accretion Regime

Authors: Elika Prameswari Fariyanto, Kazuhiro Hada, Yuzhu Cui, Mareki Honma, Masanori Nakamura, Keiichi Asada, Xuezheng Wang, Wu Jiang

Abstract: Recent advancements in high-resolution Very Long Baseline Interferometry (VLBI) have significantly improved our understanding of jet collimation near supermassive black holes in active galactic nuclei (AGNs), particularly in high-power systems. However, the collimation properties of jets in low-luminosity AGNs (LLAGNs) remain poorly explored. In this study, we investigate the jet structure of M84,… ▽ More Recent advancements in high-resolution Very Long Baseline Interferometry (VLBI) have significantly improved our understanding of jet collimation near supermassive black holes in active galactic nuclei (AGNs), particularly in high-power systems. However, the collimation properties of jets in low-luminosity AGNs (LLAGNs) remain poorly explored. In this study, we investigate the jet structure of M84, a nearby radio galaxy and a representative LLAGN, to probe jet collimation properties in a low-accretion regime. Utilizing astrometric phase-referencing observations from the Very Long Baseline Array (VLBA), supplemented by archival Very Large Array (VLA) data, we trace the jet geometry of M84 over a broad range of scales, from approximately 10^2 to 10^7 Schwarzschild radii (rs). Our analysis reveals a well-defined transition from a semi-parabolic profile, W(r) proportional to r^0.71, to a conical shape, W(r) proportional to r^1.16, occurring at approximately 1.67 x 10^4 rs. This indicates that the M84 jet is notably less collimated than those in other known LLAGN sources. Our findings provide new insights into the relationship between jet collimation and accretion rate, offering crucial constraints for jet formation models in LLAGNs. △ Less

Submitted 28 July, 2025; originally announced July 2025.

The Giant Radio Array for Neutrino Detection (GRAND) Collaboration -- Contributions to the 39th International Cosmic Ray Conference (ICRC 2025)

Authors: Jaime Álvarez-Muñiz, Rafael Alves Batista, Aurélien Benoit-Lévy, Teresa Bister, Martina Bohacova, Mauricio Bustamante, Washington Carvalho Jr., Yiren Chen, LingMei Cheng, Simon Chiche, Jean-Marc Colley, Pablo Correa, Nicoleta Cucu Laurenciu, Zigao Dai, Rogerio M. de Almeida, Beatriz de Errico, João R. T. de Mello Neto, Krijn D. de Vries, Valentin Decoene, Peter B. Denton, Bohao Duan, Kaikai Duan, Ralph Engel, William Erba, Yizhong Fan , et al. (113 additional authors not shown)

Abstract: The Giant Radio Array for Neutrino Detection (GRAND) is an envisioned observatory of ultra-high-energy particles of cosmic origin, with energies in excess of 100 PeV. GRAND uses large surface arrays of antennas to look for the radio emission from extensive air showers that are triggered by the interaction of ultra-high-energy cosmic rays, gamma rays, and neutrinos in the atmosphere or underground.… ▽ More The Giant Radio Array for Neutrino Detection (GRAND) is an envisioned observatory of ultra-high-energy particles of cosmic origin, with energies in excess of 100 PeV. GRAND uses large surface arrays of antennas to look for the radio emission from extensive air showers that are triggered by the interaction of ultra-high-energy cosmic rays, gamma rays, and neutrinos in the atmosphere or underground. In particular, for ultra-high-energy neutrinos, the future final phase of GRAND aims to be sensitive enough to detect them in spite of their plausibly tiny flux. Three prototype GRAND radio arrays have been in operation since 2023: GRANDProto300, in China, GRAND@Auger, in Argentina, and GRAND@Nançay, in France. Their goals are to field-test the GRAND detection units, understand the radio background to which they are exposed, and develop tools for diagnostic, data gathering, and data analysis. This list of contributions to the 39th International Cosmic Ray Conference (ICRC 2025) presents an overview of GRAND, in its present and future incarnations, and a first look at data collected by GRANDProto300 and GRAND@Auger, including the first cosmic-ray candidates detected by them. △ Less

Submitted 13 July, 2025; originally announced July 2025.

Comments: Note: To access the list of contributions, please follow the "HTML" link that can be found on the arXiv page

Langmuir Wave Excitation in Solar-wind Magnetic Holes

Authors: Jingting Liu, Daniel Verscharen, Jesse Coburn, Georgios Nicolaou, Xiangyu Wu, Wence Jiang, Oreste Pezzi, Francesco Pucci, Matteo Zuin, Christopher J. Owen, Hamish Reid

Abstract: Magnetic holes are structures commonly observed in various space plasma environments throughout the solar system, including the solar wind. These structures are characterized by a localized decrease in magnetic field strength, coincident with an increase in plasma density. Previous observational studies in the solar wind link the presence of Langmuir waves to magnetic holes, suggesting a strong co… ▽ More Magnetic holes are structures commonly observed in various space plasma environments throughout the solar system, including the solar wind. These structures are characterized by a localized decrease in magnetic field strength, coincident with an increase in plasma density. Previous observational studies in the solar wind link the presence of Langmuir waves to magnetic holes, suggesting a strong correlation between these phenomena. We develop a model based on magnetic-moment conservation and its violation to explain the excitation of Langmuir waves in magnetic holes. Our model illustrates that magnetic holes induce changes in the electron velocity distribution function that emit electrostatic Langmuir waves due to the bump-on-tail instability. Using data from the Solar Orbiter spacecraft, we provide a comprehensive analysis of this process and test our predictions with observations. The consistency between the model and observations indicates that our proposed process is a viable mechanism for producing Langmuir waves in magnetic holes in the solar wind. △ Less

Submitted 2 July, 2025; originally announced July 2025.

The enhanced X-ray Timing and Polarimetry mission -- eXTP for launch in 2030

Authors: Shuang-Nan Zhang, Andrea Santangelo, Yupeng Xu, Hua Feng, Fangjun Lu, Yong Chen, Mingyu Ge, Kirpal Nandra, Xin Wu, Marco Feroci, Margarita Hernanz, Congzhan Liu, Huilin He, Yusa Wang, Weichun Jiang, Weiwei Cui, Yanji Yang, Juan Wang, Wei Li, Xiaohua Liu, Bin Meng, Xiangyang Wen, Aimei Zhang, Jia Ma, Maoshun Li , et al. (136 additional authors not shown)

Abstract: In this paper we present the current status of the enhanced X-ray Timing and Polarimetry mission, which has been fully approved for launch in 2030. eXTP is a space science mission designed to study fundamental physics under extreme conditions of matter density, gravity, and magnetism. The mission aims at determining the equation of state of matter at supra-nuclear density, measuring the effects of… ▽ More In this paper we present the current status of the enhanced X-ray Timing and Polarimetry mission, which has been fully approved for launch in 2030. eXTP is a space science mission designed to study fundamental physics under extreme conditions of matter density, gravity, and magnetism. The mission aims at determining the equation of state of matter at supra-nuclear density, measuring the effects of quantum electro-dynamics, and understanding the dynamics of matter in strong-field gravity. In addition to investigating fundamental physics, the eXTP mission is poised to become a leading observatory for time-domain and multi-messenger astronomy in the 2030's, as well as providing observations of unprecedented quality on a variety of galactic and extragalactic objects. After briefly introducing the history and a summary of the scientific objectives of the eXTP mission, this paper presents a comprehensive overview of: 1) the cutting-edge technology, technical specifications, and anticipated performance of the mission's scientific instruments; 2) the full mission profile, encompassing spacecraft design, operational capabilities, and ground segment infrastructure. △ Less

Submitted 8 September, 2025; v1 submitted 9 June, 2025; originally announced June 2025.

Comments: accepted for publication in the SCIENCE CHINA Physics, Mechanics & Astronomy

All-sky search for individual Primordial Black Hole bursts with LHAASO

Authors: Zhen Cao, F. Aharonian, Y. X. Bai, Y. W. Bao, D. Bastieri, X. J. Bi, Y. J. Bi, W. Bian, A. V. Bukevich, C. M. Cai, W. Y. Cao, Zhe Cao, J. Chang, J. F. Chang, A. M. Chen, E. S. Chen, G. H. Chen, H. X. Chen, Liang Chen, Long Chen, M. J. Chen, M. L. Chen, Q. H. Chen, S. Chen, S. H. Chen , et al. (293 additional authors not shown)

Abstract: Primordial Black Holes~(PBHs) are hypothetical black holes with a wide range of masses that formed in the early universe. As a result, they may play an important cosmological role and provide a unique probe of the early universe. A PBH with an initial mass of approximately $10^{15}$~g is expected to explode today in a final burst of Hawking radiation. In this work, we conduct an all-sky search for… ▽ More Primordial Black Holes~(PBHs) are hypothetical black holes with a wide range of masses that formed in the early universe. As a result, they may play an important cosmological role and provide a unique probe of the early universe. A PBH with an initial mass of approximately $10^{15}$~g is expected to explode today in a final burst of Hawking radiation. In this work, we conduct an all-sky search for individual PBH burst events using the data collected from March 2021 to July 2024 by the Water Cherenkov Detector Array of the Large High Altitude Air Shower Observatory (LHAASO). Three PBH burst durations, 10~s, 20~s, and 100~s, are searched, with no significant PBH bursts observed. The upper limit on the local PBH burst rate density is set to be as low as 181~pc$^{-3}$~yr$^{-1}$ at 99$\%$ confidence level, representing the most stringent limit achieved to date. △ Less

Submitted 2 November, 2025; v1 submitted 30 May, 2025; originally announced May 2025.

Comments: 8 pages, 2 figures

Precise measurements of the cosmic ray proton energy spectrum in the "knee'' region

Authors: The LHAASO Collaboration, Zhen Cao, F. Aharonian, Y. X. Bai, Y. W. Bao, D. Bastieri, X. J. Bi, Y. J. Bi, W. Bian, A. V. Bukevich, C. M. Cai, W. Y. Cao, Zhe Cao, J. Chang, J. F. Chang, A. M. Chen, E. S. Chen, G. H. Chen, H. X. Chen, Liang Chen, Long Chen, M. J. Chen, M. L. Chen, Q. H. Chen, S. Chen , et al. (292 additional authors not shown)

Abstract: We report the high-purity identification of cosmic-ray (CR) protons and a precise measurement of their energy spectrum from 0.15 to 12 PeV using the Large High Altitude Air Shower Observatory (LHAASO). Abundant event statistics, combined with the simultaneous detection of electrons/photons, muons, and Cherenkov light in air showers, enable spectroscopic measurements with statistical and systematic… ▽ More We report the high-purity identification of cosmic-ray (CR) protons and a precise measurement of their energy spectrum from 0.15 to 12 PeV using the Large High Altitude Air Shower Observatory (LHAASO). Abundant event statistics, combined with the simultaneous detection of electrons/photons, muons, and Cherenkov light in air showers, enable spectroscopic measurements with statistical and systematic precision comparable to satellite data at lower energies. The proton spectrum shows significant hardening relative to low-energy extrapolations, culminating at 3 PeV, followed by sharp softening. This distinct spectral structure closely aligned with the knee in the all-particle spectrum points to the emergence of a new CR component at PeV energies that might be linked to the dozens of PeVatrons recently discovered by LHAASO, and offers crucial clues to the origin of Galactic cosmic rays. △ Less

Submitted 24 December, 2025; v1 submitted 20 May, 2025; originally announced May 2025.

Journal ref: Science Bulletin 70(2025)4173-4180

Origin of the ring ellipticity in the black hole images of M87*

Authors: Rohan Dahale, Ilje Cho, Kotaro Moriyama, Kaj Wiik, Paul Tiede, José L. Gómez, Chi-kwan Chan, Roman Gold, Vadim Y. Bernshteyn, Marianna Foschi, Britton Jeter, Hung-Yi Pu, Boris Georgiev, Abhishek V. Joshi, Alejandro Cruz-Osorio, Iniyan Natarajan, Avery E. Broderick, León D. S. Salas, Koushik Chatterjee, Kazunori Akiyama, Ezequiel Albentosa-Ruíz, Antxon Alberdi, Walter Alef, Juan Carlos Algaba, Richard Anantua , et al. (251 additional authors not shown)

Abstract: We investigate the origin of the elliptical ring structure observed in the images of the supermassive black hole M87*, aiming to disentangle contributions from gravitational, astrophysical, and imaging effects. Leveraging the enhanced capabilities of the Event Horizon Telescope (EHT) 2018 array, including improved $(u,v)$-coverage from the Greenland Telescope, we measure the ring's ellipticity usi… ▽ More We investigate the origin of the elliptical ring structure observed in the images of the supermassive black hole M87*, aiming to disentangle contributions from gravitational, astrophysical, and imaging effects. Leveraging the enhanced capabilities of the Event Horizon Telescope (EHT) 2018 array, including improved $(u,v)$-coverage from the Greenland Telescope, we measure the ring's ellipticity using five independent imaging methods, obtaining a consistent average value of $τ= 0.08_{-0.02}^{+0.03}$ with a position angle $ξ= 50.1_{-7.6}^{+6.2}$ degrees. To interpret this measurement, we compare against General Relativistic Magnetohydrodynamic (GRMHD) simulations spanning a wide range of physical parameters including thermal or non-thermal electron distribution function, spins, and ion-to-electron temperature ratios in both low and high-density regions. We find no statistically significant correlation between spin and ellipticity in GRMHD images. Instead, we identify a correlation between ellipticity and the fraction of non-ring emission, particularly in non-thermal models and models with higher jet emission. These results indicate that the ellipticity measured from the \m87 emission structure is consistent with that expected from simulations of turbulent accretion flows around black holes, where it is dominated by astrophysical effects rather than gravitational ones. Future high-resolution imaging, including space very long baseline interferometry and long-term monitoring, will be essential to isolate gravitational signatures from astrophysical effects. △ Less

Submitted 15 May, 2025; originally announced May 2025.

Comments: 18 pages, 13 figures

Journal ref: A&A 699, A279 (2025)

Measurement of separate electron and positron spectra from 10 GeV to 20GeV with the geomagnetic field on DAMPE

Authors: DAMPE Collaboration, F. Alemanno, Q. An, P. Azzarello, F. C. T. Barbato, P. Bernardini, X. J. Bi, H. Boutin, I. Cagnoli, M. S. Cai, E. Casilli, E. Catanzani, J. Chang, D. Y. Chen, J. L. Chen, Z. F. Chen, Z. X. Chen, P. Coppin, M. Y. Cui, T. S. Cui, Y. X. Cui, I. DeMitri, F. dePalma, A. DiGiovanni, T. K. Dong , et al. (127 additional authors not shown)

Abstract: The cosmic-ray (CR) electrons and positrons in space are of great significance for studying the origin and propagation of cosmic-rays. The satellite-borne experiment DArk Matter Particle Explorer (DAMPE) has been used to measure the separate electron and positron spectra, as well as the positron fraction. In this work, the Earth's magnetic field is used to distinguish CR electrons and positrons, a… ▽ More The cosmic-ray (CR) electrons and positrons in space are of great significance for studying the origin and propagation of cosmic-rays. The satellite-borne experiment DArk Matter Particle Explorer (DAMPE) has been used to measure the separate electron and positron spectra, as well as the positron fraction. In this work, the Earth's magnetic field is used to distinguish CR electrons and positrons, as the DAMPE detector does not carry an onboard magnet. The energy range for the measurements is from 10 to 20 GeV, being currently limited at high energy by the zenith pointing orientation of DAMPE. The results are consistent with previous measurements based on the magnetic spectrometer by AMS-02 and PAMELA, while the results of Fermi-LAT seem then to be systematically shifted to larger values. △ Less

Submitted 21 August, 2025; v1 submitted 9 May, 2025; originally announced May 2025.

Comments: Accepted for publication in Chinese Physics C

Study of Ultra-High-Energy Gamma-Ray Source 1LHAASO J0056+6346u and Its Possible Origins

Authors: Zhen Cao, F. Aharonian, Y. X. Bai, Y. W. Bao, D. Bastieri, X. J. Bi, Y. J. Bi, W. Bian, J. Blunier, A. V. Bukevich, C. M. Cai, Y. Y. Cai, W. Y. Cao, Zhe Cao, J. Chang, J. F. Chang, E. S. Chen, G. H. Chen, H. K. Chen, L. F. Chen, Liang Chen, Long Chen, M. J. Chen, M. L. Chen, Q. H. Chen , et al. (304 additional authors not shown)

Abstract: We report a dedicated study of the newly discovered extended UHE $γ$-ray source 1LHAASO J0056+6346u. Analyzing 979 days of LHAASO-WCDA data and 1389 days of LHAASO-KM2A data, we observed a significant excess of $γ$-ray events with both WCDA and KM2A. Assuming a point power-law source with a fixed spectral index, the significance maps reveal excesses of ${\sim}12.65\,σ$, ${\sim}22.18\,σ$, and… ▽ More We report a dedicated study of the newly discovered extended UHE $γ$-ray source 1LHAASO J0056+6346u. Analyzing 979 days of LHAASO-WCDA data and 1389 days of LHAASO-KM2A data, we observed a significant excess of $γ$-ray events with both WCDA and KM2A. Assuming a point power-law source with a fixed spectral index, the significance maps reveal excesses of ${\sim}12.65\,σ$, ${\sim}22.18\,σ$, and ${\sim}10.24\,σ$ in the energy ranges of 1--25 TeV, 25--100 TeV, and $> 100$ TeV, respectively. We use a 3D likelihood algorithm to derive the morphological and spectral parameters, and the source is detected with significances of $12.65\,σ$ by WCDA and $25.27\,σ$ by KM2A. The best-fit positions derived from WCDA and KM2A data are (R.A. = $13.96^\circ\pm0.09^\circ$, Decl. = $63.92^\circ\pm0.05^\circ$) and (R.A. = $14.00^\circ\pm0.05^\circ$, Decl. = $63.79^\circ\pm0.02^\circ$), respectively. The angular size ($r_{39}$) of 1LHAASO J0056+6346u is $0.34^\circ\pm0.04^\circ$ at 1--25 TeV and $0.24^\circ\pm0.02^\circ$ at $> 25$ TeV. The differential flux of this UHE $γ$-ray source can be described by an exponential cutoff power-law function: $(2.67\pm0.25) \times 10^{-15} (E/20\,\text{TeV})^{-1.97\pm0.10} e^{-E/(55.1\pm7.2)\,\text{TeV}} \,\text{TeV}^{-1}\,\text{cm}^{-2}\,\text{s}^{-1}$. To explore potential sources of $γ$-ray emission, we investigated the gas distribution around 1LHAASO J0056+6346u. 1LHAASO J0056+6346u is likely to be a TeV PWN powered by an unknown pulsar, which would naturally explain both its spatial and spectral properties. Another explanation is that this UHE $γ$-ray source might be associated with gas content illuminated by a nearby CR accelerator, possibly the SNR candidate G124.0+1.4. △ Less

Submitted 30 December, 2025; v1 submitted 1 April, 2025; originally announced April 2025.

Comments: 16 pages, 8 figures, Accepted for publication in SCPMA

Journal ref: Sci. China-Phys. Mech. Astron. 68, 279505 (2025)

The Jiutian simulations for the CSST extra-galactic surveys

Authors: Jiaxin Han, Ming Li, Wenkang Jiang, Zhao Chen, Huiyuan Wang, Chengliang Wei, Feihong He, Jianhua He, Jiajun Zhang, Yu Liu, Weiguang Cui, Yizhou Gu, Qi Guo, Yipeng Jing, Xi Kang, Guoliang Li, Xiong Luo, Yu Luo, Wenxiang Pei, Yisheng Qiu, Zhenlin Tan, Lizhi Xie, Xiaohu Yang, Hao-Ran Yu, Yu Yu , et al. (1 additional authors not shown)

Abstract: We provide an overview of the Jiutian simulations, a hybrid simulation suite for the China Space Survey Telescope (CSST) extragalactic surveys. It consists of four complementary modules: the primary runs with high resolutions with the fiducial concordance cosmology, the emulator runs exploring the parameter uncertainties around the fiducial cosmology, the reconstruction runs intended for recoverin… ▽ More We provide an overview of the Jiutian simulations, a hybrid simulation suite for the China Space Survey Telescope (CSST) extragalactic surveys. It consists of four complementary modules: the primary runs with high resolutions with the fiducial concordance cosmology, the emulator runs exploring the parameter uncertainties around the fiducial cosmology, the reconstruction runs intended for recovering the observed Universe position by position, and the extension runs employing extended cosmologies beyond the standard model. For the primary runs, two independent pipelines are adopted to construct subhaloes and merger trees. On top of them, four sets of mock galaxy light-cone catalogs are produced from semi-analytical models and subhalo abundance matching, providing a variety of observational properties including galaxy SED, emission lines, lensing distortions, and mock images. The 129 emulator runs are used to train the CSST emulator, achieving one percent accuracy in predicting the matter power spectrum over $k\leq 10h{\rm Mpc}^{-1}$ and $z\leq 2$. The reconstruction runs employ a number of subgrid baryonic models to predict the evolution and galaxy population resembling certain regions in the real Universe with constrained initial conditions, enabling controlled investigation of galaxy formation on top of structure formation. The extension runs cover models with warm dark matter, $f(R)$ gravity, interacting dark energy, and nonzero neutrino masses, revealing differences in the cosmic structure under alternative cosmological models. We introduce the specifications for each run, the data products derived from them, the corresponding pipeline developments, and present some main tests. Using the primary runs, we also show that the subhalo peak mass functions of different levels are approximately universal. These simulations form a comprehensive and open library for CSST surveys and beyond. △ Less

Submitted 2 June, 2025; v1 submitted 27 March, 2025; originally announced March 2025.

Comments: minor text revisions. simulation data products will be gradually released at https://jiutian.sjtu.edu.cn

Self-similar decomposition of the hierarchical merger tree of dark matter halos

Authors: Wenkang Jiang, Jiaxin Han, Fuyu Dong, Feihong He

Abstract: In the $Λ$CDM universe, structure formation is generally not a self-similar process, while some self-similarity remains in certain statistics which can greatly simplify our description and understanding of the cosmic structures. In this work, we show that the merger tree of dark matter halos is approximately self-similar by investigating the universality of the subhalo peak mass function (PMF) des… ▽ More In the $Λ$CDM universe, structure formation is generally not a self-similar process, while some self-similarity remains in certain statistics which can greatly simplify our description and understanding of the cosmic structures. In this work, we show that the merger tree of dark matter halos is approximately self-similar by investigating the universality of the subhalo peak mass function (PMF) describing the mass distribution of progenitor halos. Using a set of cosmological simulations and identifying subhalos of different merger levels with HBT+, we verify that the level-1 subhalo PMF is close to universal across halo mass, redshift, and cosmology. This approximate self-similarity allows us to analytically derive the subhalo PMF for subhalos accreted at any level (i.e., for sub-sub...halos) through self-convolutions of the level-1 PMF, and the resulting model shows good agreement with simulation measurements. We further derive a number of analytical properties on the hierarchical origin of subhalos. We show that higher-level subhalos dominate at progressively lower peak mass in the PMF and are more likely to originate from major mergers than lower-level ones. At a given merger mass ratio, the subhalo accretion rates at each level track the growth rate of the host halo. At a fixed final mass ratio, however, subhalos of higher-level, higher-mass-ratio and in more massive haloes tend to be accreted more recently. Matching subhalo peak mass to galaxy mass, these results have direct implications on the hierarchical origin of satellite galaxies. △ Less

Submitted 11 June, 2025; v1 submitted 27 February, 2025; originally announced February 2025.

Comments: 20 pages, 13 figures