Showing 1–50 of 91 results for author: Fang, J

Constraints on Ultra Heavy Dark Matter Properties from Dwarf Spheroidal Galaxies with LHAASO Observations

Authors: Zhen Cao, F. Aharonian, Q. An, Axikegu, Y. X. Bai, Y. W. Bao, D. Bastieri, X. J. Bi, Y. J. Bi, J. T. Cai, Q. Cao, W. Y. Cao, Zhe Cao, J. Chang, J. F. Chang, A. M. Chen, E. S. Chen, Liang Chen, Lin Chen, Long Chen, M. J. Chen, M. L. Chen, Q. H. Chen, S. H. Chen, S. Z. Chen , et al. (255 additional authors not shown)

Abstract: In this work we try to search for signals generated by ultra-heavy dark matter at the Large High Altitude Air Shower Observatory (LHAASO) data. We look for possible gamma-ray by dark matter annihilation or decay from 16 dwarf spheroidal galaxies in the field of view of LHAASO. Dwarf spheroidal galaxies are among the most promising targets for indirect detection of dark matter which have low fluxes… ▽ More In this work we try to search for signals generated by ultra-heavy dark matter at the Large High Altitude Air Shower Observatory (LHAASO) data. We look for possible gamma-ray by dark matter annihilation or decay from 16 dwarf spheroidal galaxies in the field of view of LHAASO. Dwarf spheroidal galaxies are among the most promising targets for indirect detection of dark matter which have low fluxes of astrophysical $γ$-ray background while large amount of dark matter. By analyzing more than 700 days observational data at LHAASO, no significant dark matter signal from 1 TeV to 1 EeV is detected. Accordingly we derive the most stringent constraints on the ultra-heavy dark matter annihilation cross-section up to EeV. The constraints on the lifetime of dark matter in decay mode are also derived. △ Less

Submitted 12 June, 2024; originally announced June 2024.

Comments: 17 pages, 12 figures, accepted by PRL

Nonthermal Acceleration of Electrons, Positrons and Protons at a Nonrelativistic Quasiparallel Collisionless Shock

Authors: Huan Yu, Qi Xia, Jun Fang

Abstract: Energetic positrons have been observed in the interstellar medium, and high-energy positrons with relativistic energies up to approximately 1 TeV have been detected in Galactic cosmic rays. We conducted a study on the acceleration of particles, specifically positrons, in a nonrelativistic quasiparallel collisionless shock induced by a plasma consisting of protons, electrons, and positrons. The pos… ▽ More Energetic positrons have been observed in the interstellar medium, and high-energy positrons with relativistic energies up to approximately 1 TeV have been detected in Galactic cosmic rays. We conducted a study on the acceleration of particles, specifically positrons, in a nonrelativistic quasiparallel collisionless shock induced by a plasma consisting of protons, electrons, and positrons. The positron-to-proton number density ratio in the plasma is 0.1. We focused on a representative shock with a sonic Mach number of 17.1 and an Alfvénic Mach number of 16.8 in the rest frame of the shock. To investigate the acceleration mechanisms of particles including positrons in the shock, we utilized one-dimensional particle-in-cell (PIC) simulations. It was found that all three species of particles in the shock can be accelerated and exhibit power law spectra. At the shock front, a significant portion of incoming upstream particles are reflected and undergo significant energy increase, and these reflected particles can be efficiently injected into the process of diffusive shock acceleration (DSA). Moveover, the reflected positrons can be further accelerated by an electric field parallel to the magnetic field when they move along the magnetic field upstream of the shock. As a result, positrons can be preferentially accelerated to be injected in the DSA process compared to electrons. △ Less

Submitted 29 May, 2024; originally announced May 2024.

Comments: 7 pages, 7 figures, Accepted for publication in ApJ

Data quality control system and long-term performance monitor of the LHAASO-KM2A

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

Abstract: The KM2A is the largest sub-array of the Large High Altitude Air Shower Observatory (LHAASO). It consists of 5216 electromagnetic particle detectors (EDs) and 1188 muon detectors (MDs). The data recorded by the EDs and MDs are used to reconstruct primary information of cosmic ray and gamma-ray showers. This information is used for physical analysis in gamma-ray astronomy and cosmic ray physics. To… ▽ More The KM2A is the largest sub-array of the Large High Altitude Air Shower Observatory (LHAASO). It consists of 5216 electromagnetic particle detectors (EDs) and 1188 muon detectors (MDs). The data recorded by the EDs and MDs are used to reconstruct primary information of cosmic ray and gamma-ray showers. This information is used for physical analysis in gamma-ray astronomy and cosmic ray physics. To ensure the reliability of the LHAASO-KM2A data, a three-level quality control system has been established. It is used to monitor the status of detector units, stability of reconstructed parameters and the performance of the array based on observations of the Crab Nebula and Moon shadow. This paper will introduce the control system and its application on the LHAASO-KM2A data collected from August 2021 to July 2023. During this period, the pointing and angular resolution of the array were stable. From the observations of the Moon shadow and Crab Nebula, the results achieved using the two methods are consistent with each other. According to the observation of the Crab Nebula at energies from 25 TeV to 100 TeV, the time averaged pointing errors are estimated to be $-0.003^{\circ} \pm 0.005^{\circ}$ and $0.001^{\circ} \pm 0.006^{\circ}$ in the R.A. and Dec directions, respectively. △ Less

Submitted 13 June, 2024; v1 submitted 20 May, 2024; originally announced May 2024.

Comments: 15 pages, 9 figures

Discovery of Very-high-energy Gamma-ray Emissions from the Low Luminosity AGN NGC 4278 by LHAASO

Authors: Zhen Cao, F. Aharonian, Q. An, Axikegu, Y. X. Bai, Y. W. Bao, D. Bastieri, X. J. Bi, Y. J. Bi, J. T. Cai, Q. Cao, W. Y. Cao, Zhe Cao, J. Chang, J. F. Chang, A. M. Chen, E. S. Chen, Liang Chen, Lin Chen, Long Chen, M. J. Chen, M. L. Chen, Q. H. Chen, S. H. Chen, S. Z. Chen , et al. (255 additional authors not shown)

Abstract: The first source catalog of Large High Altitude Air Shower Observatory reported the detection of a very-high-energy gamma ray source, 1LHAASO J1219+2915. In this paper a further detailed study of the spectral and temporal behavior of this point-like source have been carried. The best-fit position of the TeV source ($\rm{RA}=185.05^{\circ}\pm0.04^{\circ}$, $\rm{Dec}=29.25^{\circ}\pm0.03^{\circ}$) i… ▽ More The first source catalog of Large High Altitude Air Shower Observatory reported the detection of a very-high-energy gamma ray source, 1LHAASO J1219+2915. In this paper a further detailed study of the spectral and temporal behavior of this point-like source have been carried. The best-fit position of the TeV source ($\rm{RA}=185.05^{\circ}\pm0.04^{\circ}$, $\rm{Dec}=29.25^{\circ}\pm0.03^{\circ}$) is compatible with NGC 4278 within $\sim0.03$ degree. Variation analysis shows an indication of the variability at a few months level in the TeV band, which is consistent with low frequency observations. Based on these observations, we report the detection of TeV $γ$-ray emissions from this low-luminosity AGN NGC 4278. The observations by LHAASO-WCDA during active period has a significance level of 8.8\,$σ$ with best-fit photon spectral index $\varGamma=2.56\pm0.14$ and a flux $f_{1-10\,\rm{TeV}}=(7.0\pm1.1_{\rm{sta}}\pm0.35_{\rm{syst}})\times10^{-13}\,\rm{photons\,cm^{-2}\,s^{-1}}$, or approximately $5\%$ of the Crab Nebula. The discovery of VHE from NGC 4278 indicates that the compact, weak radio jet can efficiently accelerate particles and emit TeV photons. △ Less

Submitted 13 May, 2024; originally announced May 2024.

Comments: 11 pages, 5 figures

LHAASO-KM2A detector simulation using Geant4

Authors: Zhen Cao, F. Aharonian, Q. An, Axikegu, Y. X. Bai, Y. W. Bao, D. Bastieri, X. J. Bi, Y. J. Bi, J. T. Cai, Q. Cao, W. Y. Cao, Zhe Cao, J. Chang, J. F. Chang, A. M. Chen, E. S. Chen, Liang Chen, Lin Chen, Long Chen, M. J. Chen, M. L. Chen, Q. H. Chen, S. H. Chen, S. Z. Chen , et al. (254 additional authors not shown)

Abstract: KM2A is one of the main sub-arrays of LHAASO, working on gamma ray astronomy and cosmic ray physics at energies above 10 TeV. Detector simulation is the important foundation for estimating detector performance and data analysis. It is a big challenge to simulate the KM2A detector in the framework of Geant4 due to the need to track numerous photons from a large number of detector units (>6000) with… ▽ More KM2A is one of the main sub-arrays of LHAASO, working on gamma ray astronomy and cosmic ray physics at energies above 10 TeV. Detector simulation is the important foundation for estimating detector performance and data analysis. It is a big challenge to simulate the KM2A detector in the framework of Geant4 due to the need to track numerous photons from a large number of detector units (>6000) with large altitude difference (30 m) and huge coverage (1.3 km^2). In this paper, the design of the KM2A simulation code G4KM2A based on Geant4 is introduced. The process of G4KM2A is optimized mainly in memory consumption to avoid memory overffow. Some simpliffcations are used to signiffcantly speed up the execution of G4KM2A. The running time is reduced by at least 30 times compared to full detector simulation. The particle distributions and the core/angle resolution comparison between simulation and experimental data of the full KM2A array are also presented, which show good agreement. △ Less

Submitted 7 April, 2024; originally announced April 2024.

Measurements of All-Particle Energy Spectrum and Mean Logarithmic Mass of Cosmic Rays from 0.3 to 30 PeV with LHAASO-KM2A

Authors: The LHAASO Collaboration, Zhen Cao, F. Aharonian, Q. An, A. Axikegu, Y. X. Bai, Y. W. Bao, D. Bastieri, X. J. Bi, Y. J. Bi, J. T. Cai, Q. Cao, W. Y. Cao, Zhe Cao, J. Chang, J. F. Chang, A. M. Chen, E. S. Chen, Liang Chen, Lin Chen, Long Chen, M. J. Chen, M. L. Chen, Q. H. Chen, S. H. Chen , et al. (256 additional authors not shown)

Abstract: We present the measurements of all-particle energy spectrum and mean logarithmic mass of cosmic rays in the energy range of 0.3-30 PeV using data collected from LHAASO-KM2A between September 2021 and December 2022, which is based on a nearly composition-independent energy reconstruction method, achieving unprecedented accuracy. Our analysis reveals the position of the knee at… ▽ More We present the measurements of all-particle energy spectrum and mean logarithmic mass of cosmic rays in the energy range of 0.3-30 PeV using data collected from LHAASO-KM2A between September 2021 and December 2022, which is based on a nearly composition-independent energy reconstruction method, achieving unprecedented accuracy. Our analysis reveals the position of the knee at $3.67 \pm 0.05 \pm 0.15$ PeV. Below the knee, the spectral index is found to be -$2.7413 \pm 0.0004 \pm 0.0050$, while above the knee, it is -$3.128 \pm 0.005 \pm 0.027$, with the sharpness of the transition measured with a statistical error of 2%. The mean logarithmic mass of cosmic rays is almost heavier than helium in the whole measured energy range. It decreases from 1.7 at 0.3 PeV to 1.3 at 3 PeV, representing a 24% decline following a power law with an index of -$0.1200 \pm 0.0003 \pm 0.0341$. This is equivalent to an increase in abundance of light components. Above the knee, the mean logarithmic mass exhibits a power law trend towards heavier components, which is reversal to the behavior observed in the all-particle energy spectrum. Additionally, the knee position and the change in power-law index are approximately the same. These findings suggest that the knee observed in the all-particle spectrum corresponds to the knee of the light component, rather than the medium-heavy components. △ Less

Submitted 26 March, 2024; v1 submitted 15 March, 2024; originally announced March 2024.

Comments: 8 pages, 3 figures

Journal ref: Physical Review Letters 132, 131002 (2024)

Possible GeV gamma-ray emission from the pulsar wind nebula in CTA 1

Authors: Liancheng Zhou, Keyao Wu, Yunlu Gong, Jun Fang

Abstract: We report a detection of GeV $γ$-ray emission potentially originating from the pulsar wind nebula in CTA 1 by analyzing about 15 yr of Fermi Large Area Telescope data. By selecting an energy range from 50 GeV to 1 TeV to remove contamination from the $γ$-ray pulsar PSR J0007+7303, we have discovered an extended $γ$-ray source with a TS value of $\sim$ 44.94 in the region of CTA 1. The obtained flu… ▽ More We report a detection of GeV $γ$-ray emission potentially originating from the pulsar wind nebula in CTA 1 by analyzing about 15 yr of Fermi Large Area Telescope data. By selecting an energy range from 50 GeV to 1 TeV to remove contamination from the $γ$-ray pulsar PSR J0007+7303, we have discovered an extended $γ$-ray source with a TS value of $\sim$ 44.94 in the region of CTA 1. The obtained flux is measured to be 6.71 $\pm$ 2.60 $\times$ $10^{-12}$ erg $\mathrm{cm}^{-2}$ $\mathrm{s}^{-1}$ with a spectral index of 1.61 $\pm$ 0.36, which allows for a smooth connection with the flux in the TeV band. CTA 1 is also considered to be associated with 1LHAASO J0007+7303u, which is an Ultra-High-Energy source listed in the recently published catalog of the Large High Altitude Air Shower Observatory. We assume that the radiation originates from the pulsar wind nebula and that its multi-wavelength spectral energy distribution can be explained well with a time-dependent one-zone model. △ Less

Submitted 10 March, 2024; originally announced March 2024.

Comments: 7 pages, 5 figures, MNRAS accepted

A discovery of Two Slow Pulsars with FAST: "Ronin" from the Globular Cluster M15

Authors: Dengke Zhou, Pei Wang, Di Li, Jianhua Fang, Chenchen Miao, Paulo C. C. Freire, Lei Zhang, Dandan Zhang, Huaxi Chen, Yi Feng, Yifan Xiao, Jintao Xie, Xu Zhang, Chenwu Jin, Han Wang, Yinan Ke, Xuerong Guo, Rushuang Zhao, Chenhui Niu, Weiwei Zhu, Mengyao Xue, Yabiao Wang, Jiafu Wu, Zhenye Gan, Zhongyi Sun , et al. (4 additional authors not shown)

Abstract: Globular clusters harbor numerous millisecond pulsars, but long-period pulsars ($P \gtrsim 100$ ms) are rarely found. In this study, we employed a fast folding algorithm to analyze observational data from multiple globular clusters obtained by the Five-hundred-meter Aperture Spherical radio Telescope (FAST), aiming to detect the existence of long-period pulsars. We estimated the impact of the medi… ▽ More Globular clusters harbor numerous millisecond pulsars, but long-period pulsars ($P \gtrsim 100$ ms) are rarely found. In this study, we employed a fast folding algorithm to analyze observational data from multiple globular clusters obtained by the Five-hundred-meter Aperture Spherical radio Telescope (FAST), aiming to detect the existence of long-period pulsars. We estimated the impact of the median filtering algorithm in eliminating red noise on the minimum detectable flux density ($S_{\rm min}$) of pulsars. Subsequently, we successfully discovered two isolated long-period pulsars in M15 with periods approximately equal to 1.928451 seconds and 3.960716 seconds, respectively. On the $P-\dot{P}$ diagram, both pulsars are positioned below the spin-up line, suggesting a possible history of partial recycling in X-ray binary systems disrupted by dynamical encounters later on. According to timing results, these two pulsars exhibit remarkably strong magnetic fields. If the magnetic fields were weakened during the accretion process, then a short duration of accretion might explain the strong magnetic fields of these pulsars. △ Less

Submitted 18 April, 2024; v1 submitted 10 December, 2023; originally announced December 2023.

Comments: Accepted by SCIENCE CHINA Physics, Mechanics & Astronomy

Journal ref: Sci. China-Phys. Mech. Astron. 67, 269512 (2024)

Reconstruction of a Long-term spatially Contiguous Solar-Induced Fluorescence (LCSIF) over 1982-2022

Authors: Jianing Fang, Xu Lian, Youngryel Ryu, Sungchan Jeong, Chongya Jiang, Pierre Gentine

Abstract: Satellite-observed solar-induced chlorophyll fluorescence (SIF) is a powerful proxy for diagnosing the photosynthetic characteristics of terrestrial ecosystems. Despite the increasing spatial and temporal resolutions of these satellite retrievals, records of SIF are primarily limited to the recent decade, impeding their application in detecting long-term dynamics of ecosystem function and structur… ▽ More Satellite-observed solar-induced chlorophyll fluorescence (SIF) is a powerful proxy for diagnosing the photosynthetic characteristics of terrestrial ecosystems. Despite the increasing spatial and temporal resolutions of these satellite retrievals, records of SIF are primarily limited to the recent decade, impeding their application in detecting long-term dynamics of ecosystem function and structure. In this study, we leverage the two surface reflectance bands (red and near-infrared) available both from Advanced Very High-Resolution Radiometer (AVHRR, 1982-2022) and MODerate-resolution Imaging Spectroradiometer (MODIS, 2001-2022). Importantly, we calibrate and orbit-correct the AVHRR bands against their MODIS counterparts during their overlapping period. Using the long-term bias-corrected reflectance data, a neural network is then built to reproduce the Orbiting Carbon Observatory-2 SIF using AVHRR and MODIS, and used to map SIF globally over the entire 1982-2022 period. Compared with the previous MODIS-based CSIF product relying on four reflectance bands, our two-band-based product has similar skill but can be advantageously extended to the bias-corrected AVHRR period. Further comparison with three widely used vegetation indices (NDVI, kNDVI, NIRv; all based empirically on red and near-infrared bands) shows a higher or comparable correlation of LCSIF with satellite SIF and site-level GPP estimates across vegetation types, ensuring a greater capacity of LCSIF for representing terrestrial photosynthesis. Globally, LCSIF-AVHRR shows an accelerating upward trend since 1982, with an average rate of 0.0025 mW m-2 nm-1 sr-1 per decade during 1982-2000 and 0.0038 mW m-2 nm-1 sr-1 per decade during 2001-2022. Our LCSIF data provide opportunities to better understand the long-term dynamics of ecosystem photosynthesis and their underlying driving processes. △ Less

Submitted 19 June, 2024; v1 submitted 25 November, 2023; originally announced November 2023.

Does or did the supernova remnant Cassiopeia A operate as a PeVatron?

Authors: Zhen Cao, F. Aharonian, Q. An, Axikegu, Y. X. Bai, Y. W. Bao, D. Bastieri, X. J. Bi, Y. J. Bi, J. T. Cai, Q. Cao, W. Y. Cao, Zhe Cao, J. Chang, J. F. Chang, A. M. Chen, E. S. Chen, Liang Chen, Lin Chen, Long Chen, M. J. Chen, M. L. Chen, Q. H. Chen, S. H. Chen, S. Z. Chen , et al. (255 additional authors not shown)

Abstract: For decades, supernova remnants (SNRs) have been considered the prime sources of Galactic Cosmic rays (CRs). But whether SNRs can accelerate CR protons to PeV energies and thus dominate CR flux up to the knee is currently under intensive theoretical and phenomenological debate. The direct test of the ability of SNRs to operate as CR PeVatrons can be provided by ultrahigh-energy (UHE;… ▽ More For decades, supernova remnants (SNRs) have been considered the prime sources of Galactic Cosmic rays (CRs). But whether SNRs can accelerate CR protons to PeV energies and thus dominate CR flux up to the knee is currently under intensive theoretical and phenomenological debate. The direct test of the ability of SNRs to operate as CR PeVatrons can be provided by ultrahigh-energy (UHE; $E_γ\geq 100$~TeV) $γ$-rays. In this context, the historical SNR Cassiopeia A (Cas A) is considered one of the most promising target for UHE observations. This paper presents the observation of Cas A and its vicinity by the LHAASO KM2A detector. The exceptional sensitivity of LHAASO KM2A in the UHE band, combined with the young age of Cas A, enabled us to derive stringent model-independent limits on the energy budget of UHE protons and nuclei accelerated by Cas A at any epoch after the explosion. The results challenge the prevailing paradigm that Cas A-type SNRs are major suppliers of PeV CRs in the Milky Way. △ Less

Submitted 25 October, 2023; originally announced October 2023.

Comments: 11 pages, 3 figures, Accepted by the APJL

Very high energy gamma-ray emission beyond 10 TeV from GRB 221009A

Authors: Zhen Cao, F. Aharonian, Q. An, A. Axikegu, Y. X. Bai, Y. W. Bao, D. Bastieri, X. J. Bi, Y. J. Bi, J. T. Cai, Q. Cao, W. Y. Cao, Zhe Cao, J. Chang, J. F. Chang, A. M. Chen, E. S. Chen, Liang Chen, Lin Chen, Long Chen, M. J. Chen, M. L. Chen, Q. H. Chen, S. H. Chen, S. Z. Chen , et al. (255 additional authors not shown)

Abstract: The highest energy gamma-rays from gamma-ray bursts (GRBs) have important implications for their radiation mechanism. Here we report for the first time the detection of gamma-rays up to 13 TeV from the brightest GRB 221009A by the Large High Altitude Air-shower Observatory (LHAASO). The LHAASO-KM2A detector registered more than 140 gamma-rays with energies above 3 TeV during 230$-$900s after the t… ▽ More The highest energy gamma-rays from gamma-ray bursts (GRBs) have important implications for their radiation mechanism. Here we report for the first time the detection of gamma-rays up to 13 TeV from the brightest GRB 221009A by the Large High Altitude Air-shower Observatory (LHAASO). The LHAASO-KM2A detector registered more than 140 gamma-rays with energies above 3 TeV during 230$-$900s after the trigger. The intrinsic energy spectrum of gamma-rays can be described by a power-law after correcting for extragalactic background light (EBL) absorption. Such a hard spectrum challenges the synchrotron self-Compton (SSC) scenario of relativistic electrons for the afterglow emission above several TeV. Observations of gamma-rays up to 13 TeV from a source with a measured redshift of z=0.151 hints more transparency in intergalactic space than previously expected. Alternatively, one may invoke new physics such as Lorentz Invariance Violation (LIV) or an axion origin of very high energy (VHE) signals. △ Less

Submitted 22 November, 2023; v1 submitted 13 October, 2023; originally announced October 2023.

Comments: 49pages, 11figures

Journal ref: Science Advances, 9, eadj2778 (2023) 15 November 2023

Real-time Monitoring for the Next Core-Collapse Supernova in JUNO

Authors: Angel Abusleme, Thomas Adam, Shakeel Ahmad, Rizwan Ahmed, Sebastiano Aiello, Muhammad Akram, Abid Aleem, Fengpeng An, Qi An, Giuseppe Andronico, Nikolay Anfimov, Vito Antonelli, Tatiana Antoshkina, Burin Asavapibhop, João Pedro Athayde Marcondes de André, Didier Auguste, Weidong Bai, Nikita Balashov, Wander Baldini, Andrea Barresi, Davide Basilico, Eric Baussan, Marco Bellato, Marco Beretta, Antonio Bergnoli , et al. (606 additional authors not shown)

Abstract: The core-collapse supernova (CCSN) is considered one of the most energetic astrophysical events in the universe. The early and prompt detection of neutrinos before (pre-SN) and during the supernova (SN) burst presents a unique opportunity for multi-messenger observations of CCSN events. In this study, we describe the monitoring concept and present the sensitivity of the system to pre-SN and SN neu… ▽ More The core-collapse supernova (CCSN) is considered one of the most energetic astrophysical events in the universe. The early and prompt detection of neutrinos before (pre-SN) and during the supernova (SN) burst presents a unique opportunity for multi-messenger observations of CCSN events. In this study, we describe the monitoring concept and present the sensitivity of the system to pre-SN and SN neutrinos at the Jiangmen Underground Neutrino Observatory (JUNO), a 20 kton liquid scintillator detector currently under construction in South China. The real-time monitoring system is designed to ensure both prompt alert speed and comprehensive coverage of progenitor stars. It incorporates prompt monitors on the electronic board as well as online monitors at the data acquisition stage. Assuming a false alert rate of 1 per year, this monitoring system exhibits sensitivity to pre-SN neutrinos up to a distance of approximately 1.6 (0.9) kiloparsecs and SN neutrinos up to about 370 (360) kiloparsecs for a progenitor mass of 30 solar masses, considering both normal and inverted mass ordering scenarios. The pointing ability of the CCSN is evaluated by analyzing the accumulated event anisotropy of inverse beta decay interactions from pre-SN or SN neutrinos. This, along with the early alert, can play a crucial role in facilitating follow-up multi-messenger observations of the next galactic or nearby extragalactic CCSN. △ Less

Submitted 4 December, 2023; v1 submitted 13 September, 2023; originally announced September 2023.

Comments: 24 pages, 9 figures, accepted for the publication at JCAP

Observation of gamma rays up to 320 TeV from the middle-aged TeV pulsar wind nebula HESS J1849$-$000

Authors: M. Amenomori, S. Asano, Y. W. Bao, X. J. Bi, D. Chen, T. L. Chen, W. Y. Chen, Xu Chen, Y. Chen, Cirennima, S. W. Cui, Danzengluobu, L. K. Ding, J. H. Fang, K. Fang, C. F. Feng, Zhaoyang Feng, Z. Y. Feng, Qi Gao, A. Gomi, Q. B. Gou, Y. Q. Guo, Y. Y. Guo, Y. Hayashi, H. H. He , et al. (93 additional authors not shown)

Abstract: Gamma rays from HESS J1849$-$000, a middle-aged TeV pulsar wind nebula (PWN), are observed by the Tibet air shower array and the muon detector array. The detection significance of gamma rays reaches $4.0\, σ$ and $4.4\, σ$ levels above 25 TeV and 100 TeV, respectively, in units of Gaussian standard deviation $σ$. The energy spectrum measured between $40\, {\rm TeV} < E < 320\, {\rm TeV}$ for the f… ▽ More Gamma rays from HESS J1849$-$000, a middle-aged TeV pulsar wind nebula (PWN), are observed by the Tibet air shower array and the muon detector array. The detection significance of gamma rays reaches $4.0\, σ$ and $4.4\, σ$ levels above 25 TeV and 100 TeV, respectively, in units of Gaussian standard deviation $σ$. The energy spectrum measured between $40\, {\rm TeV} < E < 320\, {\rm TeV}$ for the first time is described with a simple power-law function of ${\rm d}N/{\rm d}E = (2.86 \pm 1.44) \times 10^{-16}(E/40\, {\rm TeV})^{-2.24 \pm 0.41}\, {\rm TeV}^{-1}\, {\rm cm}^{-2}\, {\rm s}^{-1}$. The gamma-ray energy spectrum from the sub-TeV ($E < 1\, {\rm TeV}$) to sub-PeV ($100\, {\rm TeV} < E < 1\, {\rm PeV}$) ranges including the results of previous studies can be modeled with the leptonic scenario, inverse Compton scattering by high-energy electrons accelerated by the PWN of PSR J1849$-$0001. On the other hand, the gamma-ray energy spectrum can also be modeled with the hadronic scenario in which gamma rays are generated from the decay of neutral pions produced by collisions between accelerated cosmic-ray protons and the ambient molecular cloud found in the gamma-ray emitting region. The cutoff energy of cosmic-ray protons $E_{\rm p\, cut}$, cut is estimated at ${\rm log}_{10}(E_{\rm p,\, cut}/{\rm TeV}) = 3.73^{+2.98}_{-0.66}$, suggesting that protons are accelerated up to the PeV energy range. Our study thus proposes that HESS J1849$-$000 should be further investigated as a new candidate for a Galactic PeV cosmic-ray accelerator, PeVatron. △ Less

Submitted 26 August, 2023; originally announced August 2023.

Comments: 10 pages, 2 figures, Accepted for publication from the Astrophysical Journal

Measurement of the Gamma-Ray Energy Spectrum beyond 100 TeV from the HESS J1843$-$033 Region

Authors: M. Amenomori, S. Asano, Y. W. Bao, X. J. Bi, D. Chen, T. L. Chen, W. Y. Chen, Xu Chen, Y. Chen, Cirennima, S. W. Cui, Danzengluobu, L. K. Ding, J. H. Fang, K. Fang, C. F. Feng, Zhaoyang Feng, Z. Y. Feng, Qi Gao, A. Gomi, Q. B. Gou, Y. Q. Guo, Y. Y. Guo, H. H. He, Z. T. He , et al. (91 additional authors not shown)

Abstract: HESS J1843$-$033 is a very-high-energy gamma-ray source whose origin remains unidentified. This work presents, for the first time, the energy spectrum of gamma rays beyond $100\, {\rm TeV}$ from the HESS J1843$-$033 region using the data recorded by the Tibet air shower array and its underground muon detector array. A gamma-ray source with an extension of $0.34^{\circ} \pm 0.12^{\circ}$ is success… ▽ More HESS J1843$-$033 is a very-high-energy gamma-ray source whose origin remains unidentified. This work presents, for the first time, the energy spectrum of gamma rays beyond $100\, {\rm TeV}$ from the HESS J1843$-$033 region using the data recorded by the Tibet air shower array and its underground muon detector array. A gamma-ray source with an extension of $0.34^{\circ} \pm 0.12^{\circ}$ is successfully detected above $25\, {\rm TeV}$ at $(α,\, δ) = (281.09^{\circ}\pm 0.10^{\circ},\, -3.76^{\circ}\pm 0.09^{\circ})$ near HESS J1843$-$033 with a statistical significance of $6.2\, σ$, and the source is named TASG J1844$-$038. The position of TASG J1844$-$038 is consistent with those of HESS J1843$-$033, eHWC J1842$-$035, and LHAASO J1843$-$0338. The measured gamma-ray energy spectrum in $25\, {\rm TeV} < E < 130\, {\rm TeV}$ is described with ${\rm d}N/{\rm d}E = (9.70\pm 1.89)\times 10^{-16} (E/40\, {\rm TeV})^{-3.26\pm 0.30}\, {\rm TeV}^{-1} {\rm cm}^{-2} {\rm s}^{-1}$, and the spectral fit to the combined spectra of HESS J1843$-$033, LHAASO J1843$-$0338, and TASG J1844$-$038 implies the existence of a cutoff at $49.5\pm 9.0\, {\rm TeV}$. Associations of TASG J1844-038 with SNR G28.6$-$0.1 and PSR J1844-0346 are also discussed in detail for the first time. △ Less

Submitted 26 August, 2023; originally announced August 2023.

Comments: 11 pages, 4 figures, 1 table

JUNO sensitivity to the annihilation of MeV dark matter in the galactic halo

Authors: JUNO Collaboration, Angel Abusleme, Thomas Adam, Shakeel Ahmad, Rizwan Ahmed, Sebastiano Aiello, Muhammad Akram, Abid Aleem, Tsagkarakis Alexandros, Fengpeng An, Qi An, Giuseppe Andronico, Nikolay Anfimov, Vito Antonelli, Tatiana Antoshkina, Burin Asavapibhop, João Pedro Athayde Marcondes de André, Didier Auguste, Weidong Bai, Nikita Balashov, Wander Baldini, Andrea Barresi, Davide Basilico, Eric Baussan, Marco Bellato , et al. (581 additional authors not shown)

Abstract: We discuss JUNO sensitivity to the annihilation of MeV dark matter in the galactic halo via detecting inverse beta decay reactions of electron anti-neutrinos resulting from the annihilation. We study possible backgrounds to the signature, including the reactor neutrinos, diffuse supernova neutrino background, charged- and neutral-current interactions of atmospheric neutrinos, backgrounds from muon… ▽ More We discuss JUNO sensitivity to the annihilation of MeV dark matter in the galactic halo via detecting inverse beta decay reactions of electron anti-neutrinos resulting from the annihilation. We study possible backgrounds to the signature, including the reactor neutrinos, diffuse supernova neutrino background, charged- and neutral-current interactions of atmospheric neutrinos, backgrounds from muon-induced fast neutrons and cosmogenic isotopes. A fiducial volume cut, as well as the pulse shape discrimination and the muon veto are applied to suppress the above backgrounds. It is shown that JUNO sensitivity to the thermally averaged dark matter annihilation rate in 10 years of exposure would be significantly better than the present-day best limit set by Super-Kamiokande and would be comparable to that expected by Hyper-Kamiokande. △ Less

Submitted 13 September, 2023; v1 submitted 15 June, 2023; originally announced June 2023.

Comments: 25 pages, 9 figures, matches the publised version

Journal ref: JCAP 09 (2023) 001

Dark matter search with a strongly-coupled hybrid spin system

Authors: Kai Wei, Zitong Xu, Yuxuan He, Xiaolin Ma, Xing Heng, Xiaofei Huang, Wei Quan, Wei Ji, Jia Liu, Xiaoping Wang, Jiancheng Fang, Dmitry Budker

Abstract: Observational evidence suggests the existence of dark matter (DM), which comprises approximately $84.4\%$ of matter in the universe. Recent advances in tabletop quantum sensor technology have enabled searches for nongravitational interactions of DM. Our experiment named ChangE utilizes Coupled Hot Atom eNsembles to search for liGht dark mattEr and new physics. We identify a strongly-coupled hybrid… ▽ More Observational evidence suggests the existence of dark matter (DM), which comprises approximately $84.4\%$ of matter in the universe. Recent advances in tabletop quantum sensor technology have enabled searches for nongravitational interactions of DM. Our experiment named ChangE utilizes Coupled Hot Atom eNsembles to search for liGht dark mattEr and new physics. We identify a strongly-coupled hybrid spin-resonance (HSR) regime that enhances the bandwidth of $^{21}$Ne nuclear spin by three orders of magnitude while maintaining high sensitivity. In combination with a self-compensating mode (SC) for low frequencies, we present a comprehensive broadband search for axion-like dark matter with Compton frequencies in the range of $[0.01, 1000]$ Hz. We set new constraints on the DM interactions with neutrons and protons, accounting for the stochastic effect. For the axion-neutron coupling, our results reach a low value of $|g_{ann}|\le 3\times 10^{-10}$ in the frequency range $[0.02 , 4]$ Hz surpassing astrophysical limits and provide the strongest laboratory constraints in the $[10, 100]$ Hz range. For the axion-proton coupling, we offer the best terrestrial constraints for the frequency below 100 Hz. △ Less

Submitted 13 June, 2023; originally announced June 2023.

Comments: 14 pages, 5 figures

The First LHAASO Catalog of Gamma-Ray Sources

Authors: Zhen Cao, F. Aharonian, Q. An, Axikegu, Y. X. Bai, Y. W. Bao, D. Bastieri, X. J. Bi, Y. J. Bi, J. T. Cai, Q. Cao, W. Y. Cao, Zhe Cao, J. Chang, J. F. Chang, A. M. Chen, E. S. Chen, Liang Chen, Lin Chen, Long Chen, M. J. Chen, M. L. Chen, Q. H. Chen, S. H. Chen, S. Z. Chen , et al. (255 additional authors not shown)

Abstract: We present the first catalog of very-high energy and ultra-high energy gamma-ray sources detected by the Large High Altitude Air Shower Observatory (LHAASO). The catalog was compiled using 508 days of data collected by the Water Cherenkov Detector Array (WCDA) from March 2021 to September 2022 and 933 days of data recorded by the Kilometer Squared Array (KM2A) from January 2020 to September 2022.… ▽ More We present the first catalog of very-high energy and ultra-high energy gamma-ray sources detected by the Large High Altitude Air Shower Observatory (LHAASO). The catalog was compiled using 508 days of data collected by the Water Cherenkov Detector Array (WCDA) from March 2021 to September 2022 and 933 days of data recorded by the Kilometer Squared Array (KM2A) from January 2020 to September 2022. This catalog represents the main result from the most sensitive large coverage gamma-ray survey of the sky above 1 TeV, covering declination from $-$20$^{\circ}$ to 80$^{\circ}$. In total, the catalog contains 90 sources with an extended size smaller than $2^\circ$ and a significance of detection at $> 5σ$. Based on our source association criteria, 32 new TeV sources are proposed in this study. Among the 90 sources, 43 sources are detected with ultra-high energy ($E > 100$ TeV) emission at $> 4σ$ significance level. We provide the position, extension, and spectral characteristics of all the sources in this catalog. △ Less

Submitted 27 November, 2023; v1 submitted 26 May, 2023; originally announced May 2023.

Comments: 40 pages, 13 figures, 4 tables

Journal ref: The Astrophysical Journal Supplement Series, 271 (2024) 25

Measurement of ultra-high-energy diffuse gamma-ray emission of the Galactic plane from 10 TeV to 1 PeV with LHAASO-KM2A

Authors: Zhen Cao, F. Aharonian, Q. An, Axikegu, Y. X. Bai, Y. W. Bao, D. Bastieri, X. J. Bi, Y. J. Bi, J. T. Cai, Q. Cao, W. Y. Cao, Zhe Cao, J. Chang, J. F. Chang, A. M. Chen, E. S. Chen, Liang Chen, Lin Chen, Long Chen, M. J. Chen, M. L. Chen, Q. H. Chen, S. H. Chen, S. Z. Chen , et al. (255 additional authors not shown)

Abstract: The diffuse Galactic $γ$-ray emission, mainly produced via interactions between cosmic rays and the interstellar medium and/or radiation field, is a very important probe of the distribution, propagation, and interaction of cosmic rays in the Milky Way. In this work we report the measurements of diffuse $γ$-rays from the Galactic plane between 10 TeV and 1 PeV energies, with the square kilometer ar… ▽ More The diffuse Galactic $γ$-ray emission, mainly produced via interactions between cosmic rays and the interstellar medium and/or radiation field, is a very important probe of the distribution, propagation, and interaction of cosmic rays in the Milky Way. In this work we report the measurements of diffuse $γ$-rays from the Galactic plane between 10 TeV and 1 PeV energies, with the square kilometer array of the Large High Altitude Air Shower Observatory (LHAASO). Diffuse emissions from the inner ($15^{\circ}<l<125^{\circ}$, $|b|<5^{\circ}$) and outer ($125^{\circ}<l<235^{\circ}$, $|b|<5^{\circ}$) Galactic plane are detected with $29.1σ$ and $12.7σ$ significance, respectively. The outer Galactic plane diffuse emission is detected for the first time in the very- to ultra-high-energy domain ($E>10$~TeV). The energy spectrum in the inner Galaxy regions can be described by a power-law function with an index of $-2.99\pm0.04$, which is different from the curved spectrum as expected from hadronic interactions between locally measured cosmic rays and the line-of-sight integrated gas content. Furthermore, the measured flux is higher by a factor of $\sim3$ than the prediction. A similar spectrum with an index of $-2.99\pm0.07$ is found in the outer Galaxy region, and the absolute flux for $10\lesssim E\lesssim60$ TeV is again higher than the prediction for hadronic cosmic ray interactions. The latitude distributions of the diffuse emission are consistent with the gas distribution, while the longitude distributions show clear deviation from the gas distribution. The LHAASO measurements imply that either additional emission sources exist or cosmic ray intensities have spatial variations. △ Less

Submitted 19 August, 2023; v1 submitted 9 May, 2023; originally announced May 2023.

Comments: 12 pages, 8 figures, 5 tables; accepted for publication in Physical Review Letters; source mask file provided as ancillary file

Journal ref: Phys. Rev. Lett. 131, 151001 (2023)

Two Transient Quasi-periodic Oscillations in $γ$-Ray Emission from the Blazar S4 0954+658

Authors: Yunlu Gong, Shiting Tian, Liancheng Zhou, Tingfeng Yi, Jun Fang

Abstract: In this work, we report periodicity search analyses in the gamma-ray light curve of the blazar S4 0954+658 monitoring undertaken by the Fermi Large Area Telescope (LAT). Four analytical methods and a tool are adopted to detect any periodic flux modulation and corresponding significance level, revealing that (i) a 66 d quasi-periodic oscillation (QPO) with the significance level of $> 5σ$ spanning… ▽ More In this work, we report periodicity search analyses in the gamma-ray light curve of the blazar S4 0954+658 monitoring undertaken by the Fermi Large Area Telescope (LAT). Four analytical methods and a tool are adopted to detect any periodic flux modulation and corresponding significance level, revealing that (i) a 66 d quasi-periodic oscillation (QPO) with the significance level of $> 5σ$ spanning over 600 d from 2015 to 2016 (MJD 57145--57745), resulting in continuous observation of nine cycles, which is one of the highest cycles discerned in blazar gamma-ray light curve; (ii) a possible QPO of 210 d at a moderate significance of $\sim3.5σ$ lasted for over 880 d from 2020 to 2022 (MJD 59035--59915), which lasted for four cycles. In addition, we discuss several physical models to explain the origin of the two transient QPOs and conclude that a geometrical scenario involving a plasma blob moving helically inside the jet can explain the time scale of the QPO. △ Less

Submitted 6 April, 2023; originally announced April 2023.

Comments: 10 pages, 11 figures; accepted for publication in ApJ

Investigating the radiative properties of LHAASO J1908+0621

Authors: Keyao Wu, Liancheng Zhou, Yunlu Gong, Jun Fang

Abstract: LHAASO J1908+0621 has been recently detected as a source emitting $γ$-rays with energies above 100 TeV, and the multiband observations show that a break around 1 TeV appears in the $γ$-ray spectrum. We have reanalyzed the GeV $γ$-ray properties for the 100 TeV source using 14 years of data recorded by the Fermi Large Area Telescope (Fermi-LAT). The spectrum in the energy range range 30-500 GeV has… ▽ More LHAASO J1908+0621 has been recently detected as a source emitting $γ$-rays with energies above 100 TeV, and the multiband observations show that a break around 1 TeV appears in the $γ$-ray spectrum. We have reanalyzed the GeV $γ$-ray properties for the 100 TeV source using 14 years of data recorded by the Fermi Large Area Telescope (Fermi-LAT). The spectrum in the energy range range 30-500 GeV has an index of 1.50 $\pm$ 0.26, which is much smaller than that detected in the TeV $γ$-rays. Additionally, the radiation properties of this source are investigated based on a one-zone time-dependent model. In the model, LHAASO J1908+0621 is associated with a pulsar wind nebula (PWN) powered by the pulsar PSR J1907$+$0602. High-energy particles composed of electrons and positrons are injected into the nebula. Multiband nonthermal emission is produced via synchrotron radiation and inverse Compton scattering (ICS). Taking the effect of radiative energy losses and adiabatic cooling into account, the spectral energy distribution from the model with a broken power-law for the distribution of the injected particles can explain the detected fluxes in the $γ$-ray bands. The results support that LHAASO J1908+0621 originates from the PWN powered by PSR J1907$+$0602, and the $γ$-rays with energy above 100 TeV are produced by the electrons/positrons in the nebula via ICS. △ Less

Submitted 7 December, 2022; originally announced December 2022.

Comments: 7 pages, 7 figures, MNRAS accepted

Model Independent Approach of the JUNO $^8$B Solar Neutrino Program

Authors: JUNO Collaboration, Jie Zhao, Baobiao Yue, Haoqi Lu, Yufeng Li, Jiajie Ling, Zeyuan Yu, Angel Abusleme, Thomas Adam, Shakeel Ahmad, Rizwan Ahmed, Sebastiano Aiello, Muhammad Akram, Abid Aleem, Tsagkarakis Alexandros, Fengpeng An, Qi An, Giuseppe Andronico, Nikolay Anfimov, Vito Antonelli, Tatiana Antoshkina, Burin Asavapibhop, João Pedro Athayde Marcondes de André, Didier Auguste, Weidong Bai , et al. (579 additional authors not shown)

Abstract: The physics potential of detecting $^8$B solar neutrinos will be exploited at the Jiangmen Underground Neutrino Observatory (JUNO), in a model independent manner by using three distinct channels of the charged-current (CC), neutral-current (NC) and elastic scattering (ES) interactions. Due to the largest-ever mass of $^{13}$C nuclei in the liquid-scintillator detectors and the {expected} low backg… ▽ More The physics potential of detecting $^8$B solar neutrinos will be exploited at the Jiangmen Underground Neutrino Observatory (JUNO), in a model independent manner by using three distinct channels of the charged-current (CC), neutral-current (NC) and elastic scattering (ES) interactions. Due to the largest-ever mass of $^{13}$C nuclei in the liquid-scintillator detectors and the {expected} low background level, $^8$B solar neutrinos would be observable in the CC and NC interactions on $^{13}$C for the first time. By virtue of optimized event selections and muon veto strategies, backgrounds from the accidental coincidence, muon-induced isotopes, and external backgrounds can be greatly suppressed. Excellent signal-to-background ratios can be achieved in the CC, NC and ES channels to guarantee the $^8$B solar neutrino observation. From the sensitivity studies performed in this work, we show that JUNO, with ten years of data, can reach the {1$σ$} precision levels of 5%, 8% and 20% for the $^8$B neutrino flux, $\sin^2θ_{12}$, and $Δm^2_{21}$, respectively. It would be unique and helpful to probe the details of both solar physics and neutrino physics. In addition, when combined with SNO, the world-best precision of 3% is expected for the $^8$B neutrino flux measurement. △ Less

Submitted 6 March, 2024; v1 submitted 15 October, 2022; originally announced October 2022.

Comments: 19 pages, 7 figures, accepted version to appear in The Astrophysical Journal. Yufeng Li and Jiajie Ling are corresponding authors

Journal ref: Astrophysical Journal 965 (2024) 122

Flux Variations of Cosmic Ray Air Showers Detected by LHAASO-KM2A During a Thunderstorm on 10 June 2021

Authors: LHAASO Collaboration, F. Aharonian, Q. An, Axikegu, L. X. Bai, Y. X. Bai, Y. W. Bao, D. Bastieri, X. J. Bi, Y. J. Bi, J. T. Cai, Zhe Cao, Zhen Cao, J. Chang, J. F. Chang, E. S. Chen, Liang Chen, Liang Chen, Long Chen, M. J. Chen, M. L. Chen, S. H. Chen, S. Z. Chen, T. L. Chen, X. J. Chen , et al. (248 additional authors not shown)

Abstract: The Large High Altitude Air Shower Observatory (LHAASO) has three sub-arrays, KM2A, WCDA and WFCTA. The flux variations of cosmic ray air showers were studied by analyzing the KM2A data during the thunderstorm on 10 June 2021. The number of shower events that meet the trigger conditions increases significantly in atmospheric electric fields, with maximum fractional increase of 20%. The variations… ▽ More The Large High Altitude Air Shower Observatory (LHAASO) has three sub-arrays, KM2A, WCDA and WFCTA. The flux variations of cosmic ray air showers were studied by analyzing the KM2A data during the thunderstorm on 10 June 2021. The number of shower events that meet the trigger conditions increases significantly in atmospheric electric fields, with maximum fractional increase of 20%. The variations of trigger rates (increases or decreases) are found to be strongly dependent on the primary zenith angle. The flux of secondary particles increases significantly, following a similar trend with that of the shower events. To better understand the observed behavior, Monte Carlo simulations are performed with CORSIKA and G4KM2A (a code based on GEANT4). We find that the experimental data (in saturated negative fields) are in good agreement with simulations, assuming the presence of a uniform upward electric field of 700 V/cm with a thickness of 1500 m in the atmosphere above the observation level. Due to the acceleration/deceleration and deflection by the atmospheric electric field, the number of secondary particles with energy above the detector threshold is modified, resulting in the changes in shower detection rate. △ Less

Submitted 6 December, 2022; v1 submitted 25 July, 2022; originally announced July 2022.

Comments: 18 pages, 11 figures

Journal ref: Chinese Phys. C 47 015001 (2023)

Prospects for Detecting the Diffuse Supernova Neutrino Background with JUNO

Authors: JUNO Collaboration, Angel Abusleme, Thomas Adam, Shakeel Ahmad, Rizwan Ahmed, Sebastiano Aiello, Muhammad Akram, Fengpeng An, Qi An, Giuseppe Andronico, Nikolay Anfimov, Vito Antonelli, Tatiana Antoshkina, Burin Asavapibhop, João Pedro Athayde Marcondes de André, Didier Auguste, Nikita Balashov, Wander Baldini, Andrea Barresi, Davide Basilico, Eric Baussan, Marco Bellato, Antonio Bergnoli, Thilo Birkenfeld, Sylvie Blin , et al. (577 additional authors not shown)

Abstract: We present the detection potential for the diffuse supernova neutrino background (DSNB) at the Jiangmen Underground Neutrino Observatory (JUNO), using the inverse-beta-decay (IBD) detection channel on free protons. We employ the latest information on the DSNB flux predictions, and investigate in detail the background and its reduction for the DSNB search at JUNO. The atmospheric neutrino induced n… ▽ More We present the detection potential for the diffuse supernova neutrino background (DSNB) at the Jiangmen Underground Neutrino Observatory (JUNO), using the inverse-beta-decay (IBD) detection channel on free protons. We employ the latest information on the DSNB flux predictions, and investigate in detail the background and its reduction for the DSNB search at JUNO. The atmospheric neutrino induced neutral current (NC) background turns out to be the most critical background, whose uncertainty is carefully evaluated from both the spread of model predictions and an envisaged \textit{in situ} measurement. We also make a careful study on the background suppression with the pulse shape discrimination (PSD) and triple coincidence (TC) cuts. With latest DSNB signal predictions, more realistic background evaluation and PSD efficiency optimization, and additional TC cut, JUNO can reach the significance of 3$σ$ for 3 years of data taking, and achieve better than 5$σ$ after 10 years for a reference DSNB model. In the pessimistic scenario of non-observation, JUNO would strongly improve the limits and exclude a significant region of the model parameter space. △ Less

Submitted 13 October, 2022; v1 submitted 18 May, 2022; originally announced May 2022.

Comments: 29 pages, 11 figures, final published version in JCAP

Journal ref: JCAP 10 (2022) 033

Quasi-periodic behaviour in the $γ$-ray light curve of the blazar PKS 0405-385

Authors: Yunlu Gong, Liancheng Zhou, Min Yuan, Haiyun Zhang, Tingfeng Yi, Jun Fang

Abstract: We analyze the quasi-periodic oscillation (QPO) of the historical light curve of FSRQs PKS 0405-385 detected by the Fermi LAT from August 2008 to November 2021. To identify and determine the QPO signal of PKS 0405-385 in the $γ$-ray light curve, we use four time series analysis techniques based on frequency and time domains, i.e., the Lomb-Scargle periodogram (LSP), the weighted wavelet z-transfor… ▽ More We analyze the quasi-periodic oscillation (QPO) of the historical light curve of FSRQs PKS 0405-385 detected by the Fermi LAT from August 2008 to November 2021. To identify and determine the QPO signal of PKS 0405-385 in the $γ$-ray light curve, we use four time series analysis techniques based on frequency and time domains, i.e., the Lomb-Scargle periodogram (LSP), the weighted wavelet z-transform (WWZ), the REDFIT and the epoch folding. The results show that PKS 0405-385 has a quasi-periodic behavior of $\sim$2.8 yr with the significance of $\sim$4.3$σ$ in Fermi long-term monitoring. Remarkably, we also performed QPO analysis in the G-band light curve observed from October 2014 to October 2021 using LSP and WWZ technology, and the results ($\sim$4$σ$ of significance) are consistent with the periodic detection in $γ$-ray. This may imply that the optical emission is radiated by an electron population same as the \gr\ emission. In discussing the possible mechanism of quasi-periodic behavior, either the helical motion within a jet or the supermassive black hole binary system provides a viable explanation for the QPO of 2.8 yr, and the relevant parameters have been estimated. △ Less

Submitted 4 May, 2022; originally announced May 2022.

Comments: 10 pages, 8 figures, Accepted for publication in ApJ

New Constraints on Exotic Spin-Velocity-Dependent Interactions

Authors: Kai Wei, Wei Ji, Changbo Fu, Arne Wickenbrock, Jiancheng Fang, Victor Flambaum, Dmitry Budker

Abstract: Experimental searches for new, "fifth" forces are attracting a lot of attention because they allow to test theoretical extensions to the standard model. Here, we report a new experimental search for possible fifth forces, specifically spin-and-velocity dependent forces, by using a K-Rb-$^{21}$Ne co-magnetometer and a tungsten ring featuring a high nucleon density. Taking advantage of the high sens… ▽ More Experimental searches for new, "fifth" forces are attracting a lot of attention because they allow to test theoretical extensions to the standard model. Here, we report a new experimental search for possible fifth forces, specifically spin-and-velocity dependent forces, by using a K-Rb-$^{21}$Ne co-magnetometer and a tungsten ring featuring a high nucleon density. Taking advantage of the high sensitivity of the co-magnetometer, the pseudomagnetic field from the fifth force is measured to be $<7$\,aT. This sets new limits on coupling constants for the neutron-nucleon and proton-nucleon interactions in the range of $\ge 0.1$ m. The coupling constant limits are established to be $|g_V^n|<6.6\times 10^{-11}$ and $|g_V^p|<3.0\times 10^{-10}$, which are more than one order of magnitude tighter than astronomical and cosmological limits on the coupling between the new gauge boson such as Z$'$ and standard model particles. △ Less

Submitted 14 March, 2022; originally announced March 2022.

Journal ref: Nature Communications | (2022)13:7387

Pasta phases in neutron stars under strong magnetic fields

Authors: Xiaopeng Wang, Jing Li, Jianjun Fang, Helena Pais, Constança Providência

Abstract: In the present work, we consider nuclear matter in the innermost crust of neutron stars under the presence of a strong magnetic field within the framework of a relativistic mean-field description. Two models with a different slope of the symmetry energy are considered in order to discuss the density-dependence of the equation of state on the crust structure. The non-homogeneous matter in $β$-equil… ▽ More In the present work, we consider nuclear matter in the innermost crust of neutron stars under the presence of a strong magnetic field within the framework of a relativistic mean-field description. Two models with a different slope of the symmetry energy are considered in order to discuss the density-dependence of the equation of state on the crust structure. The non-homogeneous matter in $β$-equilibrium is described within the coexisting phases method, and the effect of including the anomalous magnetic moment is discussed. Five different geometries for the pasta structures are considered. It is shown that strong magnetic fields cause an extension of the inner crust of the neutron stars, with the occurrence of a series of disconnected non-homogeneous matter regions above the one existing for a null magnetic field. Moreover, we observed that in these disconnected regions, for some values of the magnetic field, all five different cluster geometrical shapes occur, and the gas density is close to the cluster density. Also, the pressure at the neutron star crust-core transition much larger than the pressure obtained for a zero magnetic field. Another noticeable effect of the presence of strong magnetic fields is the increase of the proton fraction, favoring the appearance of protons in the gas background. △ Less

Submitted 11 February, 2022; originally announced February 2022.

Comments: 13 pages, 8 figures. Accepted for publication in PRD

Journal ref: Phys. Rev. D 105, 063004 (2022)

Kinetic Simulation on Electron, Proton and Helium Acceleration in a Nonrelativistic Quasiparallel shock

Authors: Jun Fang, Qi Xia, Shiting Tian, Liancheng Zhou, Huan Yu

Abstract: In addition to electrons and protons, nonrelativistic quasiparallel shocks are expected to possess the ability to accelerate heavy ions. The shocks in supernova remnants are generally supposed to be accelerators of the Galactic cosmic rays, which consist of many species of particles. We investigate diffusive shock acceleration (DSA) of electrons, protons and helium ions in a nonrelativistic quasip… ▽ More In addition to electrons and protons, nonrelativistic quasiparallel shocks are expected to possess the ability to accelerate heavy ions. The shocks in supernova remnants are generally supposed to be accelerators of the Galactic cosmic rays, which consist of many species of particles. We investigate diffusive shock acceleration (DSA) of electrons, protons and helium ions in a nonrelativistic quasiparallel shock through 1D particle-in-cell (PIC) simulation with a helium-to-proton number density ratio of $0.1$, which is relevant for the Galactic cosmic rays. The simulation indicates that waves can be excited by the flow of the energetic protons and helium ions upstream of the nonrelativistic quasiparallel shock with a sonic Mach number of 14 and an alfvén Mach number of 19.5 in the shock rest frame, and the charged particles are scattered by the self-generated waves and accelerated gradually. Moreover, the spectra of the charged particles downstream of the shock are thermal plus a nonthermal tail, and the acceleration is efficient with about $7\%$ and $5.4\%$ of the bulk kinetic energy transferred into the nonthermal protons and helium ions in the near downstream region at the end of the simulation, respectively. △ Less

Submitted 28 March, 2022; v1 submitted 5 February, 2022; originally announced February 2022.

Comments: 6 pages, 5 figures. Accepted for publication in MNRAS

Peta-electron volt gamma-ray emission from the Crab Nebula

Authors: The LHAASO Collaboration, Zhen Cao, F. Aharonian, Q. An, Axikegu, L. X. Bai, Y. X. Bai, Y. W. Bao, D. Bastieri, X. J. Bi, Y. J. Bi, H. Cai, J. T. Cai, Zhe Cao, J. Chang, J. F. Chang, B. M. Chen, E. S. Chen, J. Chen, Liang Chen, Liang Chen, Long Chen, M. J. Chen, M. L. Chen, Q. H. Chen , et al. (250 additional authors not shown)

Abstract: The Crab pulsar and the surrounding nebula powered by the pulsar's rotational energy through the formation and termination of a relativistic electron-positron wind is a bright source of gamma-rays carrying crucial information about this complex conglomerate. We report the detection of $γ$-rays with a spectrum showing gradual steepening over three energy decades, from $5\times 10^{-4}$ to $1.1$ pet… ▽ More The Crab pulsar and the surrounding nebula powered by the pulsar's rotational energy through the formation and termination of a relativistic electron-positron wind is a bright source of gamma-rays carrying crucial information about this complex conglomerate. We report the detection of $γ$-rays with a spectrum showing gradual steepening over three energy decades, from $5\times 10^{-4}$ to $1.1$ petaelectronvolt (PeV). The ultra-high-energy photons exhibit the presence of a PeV electron accelerator (a pevatron) with an acceleration rate exceeding 15% of the absolute theoretical limit. Assuming that unpulsed $γ$-rays are produced at the termination of the pulsar's wind, we constrain the pevatron's size, between $0.025$ and $0.1$ pc, and the magnetic field $\approx 110 μ$G. The production rate of PeV electrons, $2.5 \times 10^{36}$ erg $\rm s^{-1}$, constitutes 0.5% of the pulsar's spin-down luminosity, although we do not exclude a non-negligible contribution of PeV protons to the production of the highest energy $γ$-rays. △ Less

Submitted 11 November, 2021; originally announced November 2021.

Comments: 43 pages, 13 figures, 2 tables; Published in Science

Journal ref: Science, 2021, Vol 373, Issue 6553, pp. 425-430

Potential PeVatron supernova remnant G106.3+2.7 seen in the highest-energy gamma rays

Authors: M. Amenomori, Y. W. Bao, X. J. Bi, D. Chen, T. L. Chen, W. Y. Chen, Xu Chen, Y. Chen, Cirennima, S. W. Cui, Danzengluobu, L. K. Ding, J. H. Fang, K. Fang, C. F. Feng, Zhaoyang Feng, Z. Y. Feng, Qi Gao, Q. B. Gou, Y. Q. Guo, Y. Y. Guo, H. H. He, Z. T. He, K. Hibino, N. Hotta , et al. (70 additional authors not shown)

Abstract: Cosmic rays (protons and other atomic nuclei) are believed to gain energies of petaelectronvolts (PeV) and beyond at astrophysical particle accelerators called 'PeVatrons' inside our Galaxy. Although a characteristic feature of a PeVatron is expected to be a hard gamma-ray energy spectrum that extends beyond 100 teraelectronvolts (TeV) without a cutoff, none of the currently known sources exhibits… ▽ More Cosmic rays (protons and other atomic nuclei) are believed to gain energies of petaelectronvolts (PeV) and beyond at astrophysical particle accelerators called 'PeVatrons' inside our Galaxy. Although a characteristic feature of a PeVatron is expected to be a hard gamma-ray energy spectrum that extends beyond 100 teraelectronvolts (TeV) without a cutoff, none of the currently known sources exhibits such a spectrum due to the low maximum energy of accelerated cosmic rays or insufficient detector sensitivity around 100 TeV. Here we report the observation of gamma-ray emission from the supernova remnant G106.3+2.7 above 10 TeV. This work provides flux data points up to and above 100 TeV and indicates that the very-high-energy gamma-ray emission above 10 TeV is well correlated with a molecular cloud rather than the pulsar PSR J2229+6114. Regarding the gamma-ray emission mechanism of G106.3+2.7, this morphological feature appears to favor a hadronic origin via the π0 decay caused by accelerated relativistic protons over a leptonic one via the inverse-Compton scattering by relativistic electrons. Furthermore, we point out that an X-ray flux upper limit on the synchrotron spectrum would provide important information to firmly establish the hadronic scenario as the mechanism of particle acceleration at the source. △ Less

Submitted 7 September, 2021; originally announced September 2021.

Comments: Published in Nature Astronomy

Journal ref: Nature Astronomy, 5, 460-464 (2021)

Investigating the energy distribution of the high-energy particles in the Crab nebula

Authors: Lu Wen, Ke-Yao Wu, Huan Yu, Jun Fang

Abstract: The Crab nebula is a prominent pulsar wind nebula (PWN) detected in multiband observations ranging from radio to very high-energy (VHE) $γ$-rays. Recently, $γ$-rays with energies above $1 \mathrm{PeV}$ had been detected by the Large High Altitude Air Shower Observatory (LHAASO), and the energy of the most energetic particles in the nebula can be constrained. In this paper, we investigate the broad… ▽ More The Crab nebula is a prominent pulsar wind nebula (PWN) detected in multiband observations ranging from radio to very high-energy (VHE) $γ$-rays. Recently, $γ$-rays with energies above $1 \mathrm{PeV}$ had been detected by the Large High Altitude Air Shower Observatory (LHAASO), and the energy of the most energetic particles in the nebula can be constrained. In this paper, we investigate the broadest spectral energy distribution of the Crab nebula and the energy distribution of the electrons emitting the multiwavelength nonthermal emission based on a one-zone time-dependent model. The nebula is powered by the pulsar, and high-energy electrons/positrons with a broken power-law spectrum are continually injected in the nebula as the pulsar spins down. Multiwavelength nonthermal emission is generated by the leptons through synchrotron radiation and inverse Compton scattering. Using appropriate parameters, the detected fluxes for the nebula can be well reproduced, especially for the $γ$-rays from $10^2\,\mathrm{MeV}$ to $1\,\mathrm{PeV}$. The results show that the detected $γ$-rays can be produced by the leptons via the inverse Compton scattering, and the lower limit of the Lorentz factor of the most energetic leptons is $\sim 8.5\times10^{9}$. It can be concluded that there are electrons/positrons with energies higher than $4.3$\,PeV in the Crab nebula. △ Less

Submitted 14 September, 2021; v1 submitted 2 September, 2021; originally announced September 2021.

Comments: 8 pages, 5 figures, RAA Accepted

Gamma-ray Observation of the Cygnus Region in the 100 TeV Energy Region

Authors: M. Amenomori, Y. W. Bao, X. J. Bi, D. Chen, T. L. Chen, W. Y. Chen, Xu Chen, Y. Chen, Cirennima, S. W. Cui, Danzengluobu, L. K. Ding, J. H. Fang, K. Fang, C. F. Feng, Zhaoyang Feng, Z. Y. Feng, Qi Gao, A. Gomi, Q. B. Gou, Y. Q. Guo, Y. Y. Guo, H. H. He, Z. T. He, K. Hibino , et al. (88 additional authors not shown)

Abstract: We report observations of gamma-ray emissions with energies in the 100 TeV energy region from the Cygnus region in our Galaxy. Two sources are significantly detected in the directions of the Cygnus OB1 and OB2 associations. Based on their positional coincidences, we associate one with a pulsar PSR J2032+4127 and the other mainly with a pulsar wind nebula PWN G75.2+0.1 with the pulsar moving away f… ▽ More We report observations of gamma-ray emissions with energies in the 100 TeV energy region from the Cygnus region in our Galaxy. Two sources are significantly detected in the directions of the Cygnus OB1 and OB2 associations. Based on their positional coincidences, we associate one with a pulsar PSR J2032+4127 and the other mainly with a pulsar wind nebula PWN G75.2+0.1 with the pulsar moving away from its original birthplace situated around the centroid of the observed gamma-ray emission. This work would stimulate further studies of particle acceleration mechanisms at these gamma-ray sources. △ Less

Submitted 2 July, 2021; originally announced July 2021.

Comments: Accepted for publication in the Physical Review Letters

Strong magnetic fields: neutron stars with an extended inner crust

Authors: Helena Pais, Bruno Bertolino, Jianjun Fang, Xiaopeng Wang, Constança Providência

Abstract: Using relativistic mean-field models, the formation of clusterized matter, as the one expected to exist in the inner crust of neutron stars, is determined under the effect of strong magnetic fields. As already predicted from a calculation of the unstable modes resulting from density fluctuations at subsaturation densities, we confirm in the present work that for magnetic field intensities of the o… ▽ More Using relativistic mean-field models, the formation of clusterized matter, as the one expected to exist in the inner crust of neutron stars, is determined under the effect of strong magnetic fields. As already predicted from a calculation of the unstable modes resulting from density fluctuations at subsaturation densities, we confirm in the present work that for magnetic field intensities of the order of $\approx 5 \times 10^{16}$ G to $5 \times 10^{17}$ G, pasta phases may occur for densities well above the zero-field crust-core transition density. This confirms that the extension of the crust may be larger than expected. It is also verified that the equilibrium structure of the clusterized matter is very sensitive to the intensity of the magnetic fields. As a result, the decay of the magnetic field may give rise to internal stresses which may result on the yield and fracture of the inner crust lattice. △ Less

Submitted 1 June, 2021; originally announced June 2021.

Comments: 9 pages, 4 figures, accepted for publication in the EPJ A Topical Issue "The QCD Phase Diagram in Strong Magnetic Fields"

First Detection of sub-PeV Diffuse Gamma Rays from the Galactic Disk: Evidence for Ubiquitous Galactic Cosmic Rays beyond PeV Energies

Authors: M. Amenomori, Y. W. Bao, X. J. Bi, D. Chen, T. L. Chen, W. Y. Chen, Xu Chen, Y. Chen, Cirennima, S. W. Cui, Danzengluobu, L. K. Ding, J. H. Fang, K. Fang, C. F. Feng, Zhaoyang Feng, Z. Y. Feng, Qi Gao, Q. B. Gou, Y. Q. Guo, Y. Y. Guo, H. H. He, Z. T. He, K. Hibino, N. Hotta , et al. (70 additional authors not shown)

Abstract: We report, for the first time, the long-awaited detection of diffuse gamma rays with energies between 100 TeV and 1 PeV in the Galactic disk. Particularly, all gamma rays above 398 TeV are observed apart from known TeV gamma-ray sources and compatible with expectations from the hadronic emission scenario in which gamma rays originate from the decay of $π^0$'s produced through the interaction of pr… ▽ More We report, for the first time, the long-awaited detection of diffuse gamma rays with energies between 100 TeV and 1 PeV in the Galactic disk. Particularly, all gamma rays above 398 TeV are observed apart from known TeV gamma-ray sources and compatible with expectations from the hadronic emission scenario in which gamma rays originate from the decay of $π^0$'s produced through the interaction of protons with the interstellar medium in the Galaxy. This is strong evidence that cosmic rays are accelerated beyond PeV energies in our Galaxy and spread over the Galactic disk. △ Less

Submitted 17 May, 2021; v1 submitted 11 April, 2021; originally announced April 2021.

Comments: Accepted for publication in the Physical Review Letters

Journal ref: Phys. Rev. Lett. 126, 141101 (2021)

Calibration of the Air Shower Energy Scale of the Water and Air Cherenkov Techniques in the LHAASO experiment

Authors: F. Aharonian, Q. An, Axikegu, L. X. Bai, Y. X. Bai, Y. W. Bao, D. Bastieri, X. J. Bi, Y. J. Bi, H. Cai, J. T. Cai, Z. Cao Z. Cao, J. Chang, J. F. Chang, X. C. Chang, B. M. Chen, J. Chen, L. Chen, L. Chen, L. Chen, M. J. Chen, M. L. Chen, Q. H. Chen, S. H. Chen, S. Z. Chen , et al. (233 additional authors not shown)

Abstract: The Wide Field-of-View Cherenkov Telescope Array (WFCTA) and the Water Cherenkov Detector Arrays (WCDA) of LHAASO are designed to work in combination for measuring the energy spectra of various cosmic ray species over a very wide energy range from a few TeV to 10 PeV. The energy calibration of WCDA can be achieved with a proven technique of measuring the westward shift of the Moon shadow of galact… ▽ More The Wide Field-of-View Cherenkov Telescope Array (WFCTA) and the Water Cherenkov Detector Arrays (WCDA) of LHAASO are designed to work in combination for measuring the energy spectra of various cosmic ray species over a very wide energy range from a few TeV to 10 PeV. The energy calibration of WCDA can be achieved with a proven technique of measuring the westward shift of the Moon shadow of galactic cosmic rays due to the geomagnetic field. This deflection angle $Δ$ is inversely proportional to the energy of the cosmic rays. The precise measurements of the shifts by WCDA allows us to calibrate its energy scale for energies as high as 35 TeV. The energy scale measured by WCDA can be used to cross calibrate the energy reconstructed by WFCTA, which spans the whole energy range up to 10 PeV. In this work, we will demonstrate the feasibility of the method using the data collected from April 2019 to January 2020 by the WFCTA array and WCDA-1 detector, the first of the three water Cherenkov ponds, already commissioned at LHAASO site. △ Less

Submitted 13 April, 2021; v1 submitted 11 April, 2021; originally announced April 2021.

Journal ref: Phys. Rev. D 104, 062007 (2021)

Construction and On-site Performance of the LHAASO WFCTA Camera

Authors: F. Aharonian, Q. An, Axikegu, L. X. Bai, Y. X. Bai, Y. W. Bao, D. Bastieri, X. J. Bi, Y. J. Bi, H. Cai, J. T. Cai, Z. Cao, Z. Cao, J. Chang, J. F. Chang, X. C. Chang, B. M. Chen, J. Chen, L. Chen, L. Chen, L. Chen, M. J. Chen, M. L. Chen, Q. H. Chen, S. H. Chen , et al. (234 additional authors not shown)

Abstract: The focal plane camera is the core component of the Wide Field-of-view Cherenkov/fluorescence Telescope Array (WFCTA) of the Large High-Altitude Air Shower Observatory (LHAASO). Because of the capability of working under moonlight without aging, silicon photomultipliers (SiPM) have been proven to be not only an alternative but also an improvement to conventional photomultiplier tubes (PMT) in this… ▽ More The focal plane camera is the core component of the Wide Field-of-view Cherenkov/fluorescence Telescope Array (WFCTA) of the Large High-Altitude Air Shower Observatory (LHAASO). Because of the capability of working under moonlight without aging, silicon photomultipliers (SiPM) have been proven to be not only an alternative but also an improvement to conventional photomultiplier tubes (PMT) in this application. Eighteen SiPM-based cameras with square light funnels have been built for WFCTA. The telescopes have collected more than 100 million cosmic ray events and preliminary results indicate that these cameras are capable of working under moonlight. The characteristics of the light funnels and SiPMs pose challenges (e.g. dynamic range, dark count rate, assembly techniques). In this paper, we present the design features, manufacturing techniques and performances of these cameras. Finally, the test facilities, the test methods and results of SiPMs in the cameras are reported here. △ Less

Submitted 4 July, 2021; v1 submitted 29 December, 2020; originally announced December 2020.

Comments: 45 pages, 21 figures, article

Journal ref: Eur. Phys. J. C 81, 657 (2021)

The observation of the Crab Nebula with LHAASO-KM2A for the performance study

Authors: F. Aharonian, Q. An, Axikegu, L. X. Bai, Y. X. Bai, Y. W. Bao, D. Bastieri, X. J. Bi, Y. J. Bi, H. Cai, J. T. Cai, Z. Cao, Z. Cao, J. Chang, J. F. Chang, X. C. Chang, B. M. Chen, J. Chen, L. Chen, L. Chen, L. Chen, M. J. Chen, M. L. Chen, Q. H. Chen, S. H. Chen , et al. (234 additional authors not shown)

Abstract: As a sub-array of the Large High Altitude Air Shower Observatory (LHAASO), KM2A is mainly designed to cover a large fraction of the northern sky to hunt for gamma-ray sources at energies above 10 TeV. Even though the detector construction is still underway, a half of the KM2A array has been operating stably since the end of 2019. In this paper, we present the pipeline of KM2A data analysis and the… ▽ More As a sub-array of the Large High Altitude Air Shower Observatory (LHAASO), KM2A is mainly designed to cover a large fraction of the northern sky to hunt for gamma-ray sources at energies above 10 TeV. Even though the detector construction is still underway, a half of the KM2A array has been operating stably since the end of 2019. In this paper, we present the pipeline of KM2A data analysis and the first observation on the Crab Nebula, a standard candle in very high energy gamma-ray astronomy. We detect gamma-ray signals from the Crab Nebula in both energy ranges of 10$-$100 TeV and $>$100 TeV with high significance, by analyzing the KM2A data of 136 live days between December 2019 and May 2020. With the observations, we test the detector performance including angular resolution, pointing accuracy and cosmic ray background rejection power. The energy spectrum of the Crab Nebula in the energy range 10-250 TeV fits well with a single power-law function dN/dE =(1.13$\pm$0.05$_{stat}$$\pm$0.08$_{sys}$)$\times$10$^{-14}$$\cdot$(E/20TeV)$^{-3.09\pm0.06_{stat}\pm0.02_{sys}}$ cm$^{-2}$ s$^{-1}$ TeV$^{-1}$. It is consistent with previous measurements by other experiments. This opens a new window of gamma-ray astronomy above 0.1 PeV through which ultrahigh-energy gamma-ray new phenomena, such as cosmic PeVatrons, might be discovered. △ Less

Submitted 13 October, 2020; originally announced October 2020.

Comments: 13 pages, 15 figures,submitted to CPC

Journal ref: Chinese Phys. C 45 025002 (2021)

Numerically investigating the peculiar periphery of a supernova remnant in the medium with a density gradient: the case of RCW 103

Authors: Chun-Yan Lu, Jing-Wen Yan, Lu Wen, Jun Fang

Abstract: The young shell-type supernova remnant RCW 103 has peculiar properties in the X-ray morphology obtained with Chandra. The southeastern shell is more brighter in the X-rays, and the curved border of the shell in this region is more flatten than the other part. We investigate the formation of the peculiar periphery of the supernova remnant RCW 103 using 3D hydrodynamical simulation. Assuming that th… ▽ More The young shell-type supernova remnant RCW 103 has peculiar properties in the X-ray morphology obtained with Chandra. The southeastern shell is more brighter in the X-rays, and the curved border of the shell in this region is more flatten than the other part. We investigate the formation of the peculiar periphery of the supernova remnant RCW 103 using 3D hydrodynamical simulation. Assuming that the supernova ejecta has been evolved in the medium with a density gradient, the detected shape of the periphery can be generally reproduced. For RCW 103, with the ejecta mass of $3.0~M_{\odot}$, the density of the background material of $2.0~\mathrm{cm}^{-3}$, and a gradient of $3.3 - 4.0 ~\mathrm{cm}^{-3} \mathrm{pc}^{-1}$, the X-ray periphery can be generally reproduced. The simulation turned out that the asymmetry of the SNR RCW 103 is mainly due to the inhomogeneous medium with a density gradient. △ Less

Submitted 6 August, 2020; originally announced August 2020.

Comments: Accepted for publication in RAA,7 pages, 3 figures

Investigating the multiband nonthermal emission of the 100 TeV source eHWC J2019$+$368 with a pulsar wind nebula scenario

Authors: Jun Fang, Lu Wen, Huan Yu, Chen Songzhan

Abstract: eHWC J2019+368 is one of the sources emitting $γ$-rays with energies higher than 100 TeV based on the recent measurement with the High Altitude Water Cherenkov Observatory (HAWC), and the origin is still in debate. The pulsar PSR J2021$+$3651 is spatially coincident with the TeV source. We investigate theoretically whether the multiband nonthermal emission of eHWC J2019+368 can originate from the… ▽ More eHWC J2019+368 is one of the sources emitting $γ$-rays with energies higher than 100 TeV based on the recent measurement with the High Altitude Water Cherenkov Observatory (HAWC), and the origin is still in debate. The pulsar PSR J2021$+$3651 is spatially coincident with the TeV source. We investigate theoretically whether the multiband nonthermal emission of eHWC J2019+368 can originate from the pulsar wind nebula (PWN) G75.2$+$0.1 powered by PSR J2021$+$3651. In the model, the spin-down power of the pulsar is transferred to high-energy particles and magnetic field in the nebula. As the particles with an energy distribution of either a broken power-law or a power-law continually injected into the nebula, the multiband nonthermal emission is produced via synchrotron radiation and inverse Compton scattering. The spectral energy distribution of the nebula from the model with the reasonable parameters is generally consistent with the detected radio, X-ray and TeV $γ$-ray fluxes. Our study supports that the PWN has the ability to produce the TeV $γ$-rays of eHWC J2019+368, and the most energetic particles in the nebula have energies up to about $0.4$ PeV. △ Less

Submitted 27 July, 2020; originally announced July 2020.

Comments: 6 pages, 9 figures

Investigating the morphology of the supernova remnant G349.7+00.2 in the medium with a density gradient

Authors: Jing-Wen Yan, Chun-Yan Lu, Lu Wen, Huan Yu, Jun Fang

Abstract: G349.7 + 00.2 is a young Galactic supernova remnant (SNR) with a mushroom morphology in radio and X-rays, and it has been detected across the entire electromagnetic spectrum from radio to high energy $γ$-rays. Moreover, the remnant is interacting with a molecular cloud based on the observations in the radio and infrared band. The reason for the formation of the periphery and the dynamical evolutio… ▽ More G349.7 + 00.2 is a young Galactic supernova remnant (SNR) with a mushroom morphology in radio and X-rays, and it has been detected across the entire electromagnetic spectrum from radio to high energy $γ$-rays. Moreover, the remnant is interacting with a molecular cloud based on the observations in the radio and infrared band. The reason for the formation of the periphery and the dynamical evolution of the remnant are investigated using 3D hydrodynamical (HD) simulations. Under the assumption that the supernova ejecta is evolved in the medium with a density gradient, the shell is composed of two hemispheres with different radiuses, and the smaller hemisphere is in relatively dense media. The resulting periphery of remnant is consistent with detected ones, and it can be concluded that the peculiar periphery of G349.7+00.2 can be reproduced as the remnants interacting with the medium with a density gradient. △ Less

Submitted 15 April, 2020; originally announced April 2020.

Comments: 8 pages, 5 figures

Numerically investigating the morphology of the supernova remnant SN 1006 in the ambient medium with a density discontinuity

Authors: Jun Fang, Jingwen Yan, Lu Wen, Chunyan Lu, Huan Yu

Abstract: Multiband observations on the type Ia supernova remnant SN 1006 indicate peculiar properties in its morphologies of emission in the radio, optical and X-ray bands. In the hard X-rays, the remnant is bilateral with two opposite bright limbs with prominent protrusions. Moreover, a filament has been detected at the radio, optical and soft X-ray wavelengthes. The reason for these peculiar features in… ▽ More Multiband observations on the type Ia supernova remnant SN 1006 indicate peculiar properties in its morphologies of emission in the radio, optical and X-ray bands. In the hard X-rays, the remnant is bilateral with two opposite bright limbs with prominent protrusions. Moreover, a filament has been detected at the radio, optical and soft X-ray wavelengthes. The reason for these peculiar features in the morphologies of the remnant is investigated using 3D HD simulations. With the assumption that the supernova ejecta is evolved in the ambient medium with a density discontinuity, the radius of the remnant's boundary is smaller in the tenuous medium, and the shell consists of two hemispheres with different radiuses. Along particular line of sights, protrusions appear on the periphery of the remnants since the emission from the edge of the hemisphere with a larger radius is located outside that from the shell of the small hemisphere. Furthermore, the northwest filament of SN 1006 arises as a result of the intersection of the line of sight and the shocked material near the edges of the two hemispheres. It can be concluded that the features that the protrusions on the northeast and southwest limbs and the northwest filament in the morphologies of SN 1006 can be reproduced as the remnants interacting with the medium with a density discontinuity. △ Less

Submitted 7 November, 2019; originally announced November 2019.

Comments: 6 pages, 5 figures, submitted to MNRAS

The SFR-radius connection: data and implications for wind strength and halo concentration

Authors: Lin Lin, S. M. Faber, David C. Koo, Samir Salim, Aaron A. Dutton, Jerome J. Fang, Fangzhou Jiang, Cristoph T. Lee, Aldo Rodríguez-Puebla, A. van der Wel, Yicheng Guo, Guillermo Barro, Joel R. Primack, Avishai Dekel, Zhu Chen, Yifei Luo, Viraj Pandya, Rachel S. Somerville, Henry C. Ferguson, Susan Kassin, Anton M. Koekemoer, Norman A. Grogin, Audrey Galametz, P. Santini, Hooshang Nayyeri , et al. (4 additional authors not shown)

Abstract: This paper is one in a series that explores the importance of radius as a second parameter in galaxy evolution. The topic investigated here is the relationship between star formation rate (SFR) and galaxy radius ($R_{\rm e}$) for main-sequence star-forming galaxies. The key observational result is that, over a wide range of stellar mass and redshift in both CANDELS and SDSS, there is little trend… ▽ More This paper is one in a series that explores the importance of radius as a second parameter in galaxy evolution. The topic investigated here is the relationship between star formation rate (SFR) and galaxy radius ($R_{\rm e}$) for main-sequence star-forming galaxies. The key observational result is that, over a wide range of stellar mass and redshift in both CANDELS and SDSS, there is little trend between SFR and $R_{\rm e}$ at fixed stellar mass. The Kennicutt-Schmidt law, or any similar density-related star formation law, then implies that smaller galaxies must have lower gas fractions than larger galaxies (at fixed $M_{\ast}$), and this is supported by observations of local star-forming galaxies. We investigate the implication by adopting the equilibrium "bathtub" model: the ISM gas mass is assumed to be constant over time and the net star formation rate is the difference between the accretion rate of gas onto the galaxy from the halo and the outflow rate due to winds. To match the observed null correlation between SFR and radius, the bathtub model requires that smaller galaxies at fixed mass have weaker galactic winds. Our hypothesis is that galaxies are a 2-dimensional family whose properties are set mainly by halo mass and concentration. Galaxy radius and accretion rate plausibly both depend on halo concentration, which predicts how wind strength should vary with $R_{\rm e}$ and SFR. △ Less

Submitted 20 July, 2020; v1 submitted 24 October, 2019; originally announced October 2019.

Comments: 21 pages, 9 figures, ApJ accepted

Quenching as a Contest between Galaxy Halos and their Central Black Holes

Authors: Zhu Chen, S. M. Faber, David. C. Koo, Rachel S. Somerville, Joel R. Primack, Avishai Dekel, Aldo Rodríguez-Puebla, Yicheng Guo, Guillermo Barro, Dale D. Kocevski, A. van der Wel, Joanna Woo, Eric F. Bell, Jerome J. Fang, Henry C. Ferguson, Mauro Giavalisco, Marc Huertas-Company, Fangzhou Jiang, Susan Kassin, Lin Lin, Fengshan Liu, Yifei Luo, Zhijian Luo, Camilla Pacifici, Viraj Pandya , et al. (5 additional authors not shown)

Abstract: Existing models of galaxy formation have not yet explained striking correlations between structure and star-formation activity in galaxies, notably the sloped and moving boundaries that divide star-forming from quenched galaxies in key structural diagrams. This paper uses these and other relations to ``reverse-engineer'' the quenching process for central galaxies. The basic idea is that star-formi… ▽ More Existing models of galaxy formation have not yet explained striking correlations between structure and star-formation activity in galaxies, notably the sloped and moving boundaries that divide star-forming from quenched galaxies in key structural diagrams. This paper uses these and other relations to ``reverse-engineer'' the quenching process for central galaxies. The basic idea is that star-forming galaxies with larger radii (at a given stellar mass) have lower black-hole masses due to lower central densities. Galaxies cross into the green valley when the cumulative effective energy radiated by their black hole equals $\sim4\times$ their halo-gas binding energy. Since larger-radii galaxies have smaller black holes, one finds they must evolve to higher stellar masses in order to meet this halo-energy criterion, which explains the sloping boundaries. A possible cause of radii differences among star-forming galaxies is halo concentration. The evolutionary tracks of star-forming galaxies are nearly parallel to the green-valley boundaries, and it is mainly the sideways motions of these boundaries with cosmic time that cause galaxies to quench. BH-scaling laws for star-forming, quenched, and green-valley galaxies are different, and most BH mass growth takes place in the green valley. Implications include: the radii of star-forming galaxies are an important second parameter in shaping their black holes; black holes are connected to their halos but in different ways for star-forming, quenched, and green-valley galaxies; and the same BH-halo quenching mechanism has been in place since $z \sim 3$. We conclude with a discussion of black hole-galaxy co-evolution, the origin and interpretation of BH scaling laws. △ Less

Submitted 26 May, 2020; v1 submitted 24 September, 2019; originally announced September 2019.

Comments: 45 pages, 20 figures, 2 tables, accepted for publication in ApJ

Early acceleration of electrons and protons at the nonrelativistic quasiparallel shocks with different obliquity angles

Authors: Jun Fang, Chun-Yan Lu, Jing-Wen Yan, Huan Yu

Abstract: The early acceleration of protons and electrons in the nonrelativistic collisionless shocks with three obliquities are investigated through 1D particle-in-cell simulations. In the simulations, the charged particles possessing a velocity of $0.2\, c$ flow towards a reflecting boundary, and the shocks with a sonic Mach number of $13.4$ and a Alfven Mach number of $16.5$ in the downstream shock frame… ▽ More The early acceleration of protons and electrons in the nonrelativistic collisionless shocks with three obliquities are investigated through 1D particle-in-cell simulations. In the simulations, the charged particles possessing a velocity of $0.2\, c$ flow towards a reflecting boundary, and the shocks with a sonic Mach number of $13.4$ and a Alfven Mach number of $16.5$ in the downstream shock frame are generated. In these quasi-parallel shocks with the obliquity angles $θ= 15^\circ$, $30^\circ$, and $45^\circ$, some of the protons and the electrons can be injected into the acceleration processes, and their downstream spectra in the momentum space show a power law tail at a time of $1.89\times10^5 ω_{\rm pe}^{-1}$, where $ω_{\rm pe}$ is the electron plasma frequency. Moreover, the charged particles reflected at the shock excite magnetic waves upstream of the shock. The shock drift acceleration is more prominent with a larger obliquity angle for the shocks, but the accelerated particles diffuse parallel to the shock propagation direction more easily to participate in the diffusive shock acceleration. At the time still in the early acceleration stage, more energetic protons and electrons appear in the downstream of the shock for $θ= 15^\circ$ compared with the other two obliquities; moreover, in the upstream region, the spectrum of the accelerated electrons is the hardest for $θ_{\rm nB} = 45^\circ$ among the three obliquities, whereas the proton spectra for $θ_{\rm nB} = 15^\circ$ and $45^\circ$ are similar as a result of the competition of the effectiveness of the shock drift acceleration and the diffusive shock acceleration. △ Less

Submitted 21 August, 2019; originally announced August 2019.

Comments: 10 pages, 10 figures, accepted for publication in RAA

Structural and Stellar Population Properties vs. Bulge Types in Sloan Digital Sky Survey Central Galaxies

Authors: Yifei Luo, S. M. Faber, Aldo Rodriguez-Puebla, Joanna Woo, Yicheng Guo, David C. Koo, Joel R. Primack, Zhu Chen, Hassen M. Yesuf, Lin Lin, Guillermo Barro, Jerome J. Fang, Viraj Pandya, M. Huertas-Company, Shude Mao

Abstract: This paper studies pseudo-bulges (P-bulges) and classical bulges (C-bulges) in Sloan Digital Sky Survey central galaxies using the new bulge indicator $ΔΣ_1$, which measures relative central stellar-mass surface density within 1 kpc. We compare $ΔΣ_1$ to the established bulge-type indicator $Δ\langleμ_e\rangle$ from Gadotti (2009) and show that classifying by $ΔΣ_1$ agrees well with… ▽ More This paper studies pseudo-bulges (P-bulges) and classical bulges (C-bulges) in Sloan Digital Sky Survey central galaxies using the new bulge indicator $ΔΣ_1$, which measures relative central stellar-mass surface density within 1 kpc. We compare $ΔΣ_1$ to the established bulge-type indicator $Δ\langleμ_e\rangle$ from Gadotti (2009) and show that classifying by $ΔΣ_1$ agrees well with $Δ\langleμ_e\rangle$. $ΔΣ_1$ requires no bulge-disk decomposition and can be measured on SDSS images out to $z = 0.07$. Bulge types using it are mapped onto twenty different structural and stellar-population properties for 12,000 SDSS central galaxies with masses 10.0 < log $M_*$/$M_{\odot}$ < 10.4. New trends emerge from this large sample. Structural parameters show fairly linear log-log relations vs. $ΔΣ_1$ and $Δ\langleμ_e\rangle$ with only moderate scatter, while stellar-population parameters show a highly non-linear "elbow" in which specific star-formation rate remains roughly flat with increasing central density and then falls rapidly at the elbow, where galaxies begin to quench. P-bulges occupy the low-density end of the horizontal arm of the elbow and are universally star-forming, while C-bulges occupy the elbow and the vertical branch and exhibit a wide range of star-formation rates at fixed density. The non-linear relation between central density and star-formation rate has been seen before, but this mapping onto bulge class is new. The wide range of star-formation rates in C-bulges helps to explain why bulge classifications using different parameters have sometimes disagreed in the past. The elbow-shaped relation between density and stellar indices suggests that central structure and stellar-populations evolve at different rates as galaxies begin to quench. △ Less

Submitted 9 March, 2020; v1 submitted 21 August, 2019; originally announced August 2019.

Comments: 23 pages, 17 figures, published in MNRAS

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 493, Issue 2, 2020, p.1686-1707

First Detection of Photons with Energy Beyond 100 TeV from an Astrophysical Source

Authors: M. Amenomori, Y. W. Bao, X. J. Bi, D. Chen, T. L. Chen, W. Y. Chen, Xu Chen, Y. Chen, Cirennima, S. W. Cui, Danzengluobu, L. K. Ding, J. H. Fang, K. Fang, C. F. Feng, Zhaoyang Feng, Z. Y. Feng, Qi Gao, Q. B. Gou, Y. Q. Guo, H. H. He, Z. T. He, K. Hibino, N. Hotta, Haibing Hu , et al. (66 additional authors not shown)

Abstract: We report on the highest energy photons from the Crab Nebula observed by the Tibet air shower array with the underground water-Cherenkov-type muon detector array. Based on the criterion of muon number measured in an air shower, we successfully suppress 99.92% of the cosmic-ray background events with energies $E>100$ TeV. As a result, we observed 24 photon-like events with $E>100$ TeV against 5.5 b… ▽ More We report on the highest energy photons from the Crab Nebula observed by the Tibet air shower array with the underground water-Cherenkov-type muon detector array. Based on the criterion of muon number measured in an air shower, we successfully suppress 99.92% of the cosmic-ray background events with energies $E>100$ TeV. As a result, we observed 24 photon-like events with $E>100$ TeV against 5.5 background events, which corresponds to 5.6$σ$ statistical significance. This is the first detection of photons with $E>100$ TeV from an astrophysical source. △ Less

Submitted 13 June, 2019; originally announced June 2019.

Comments: April 4, 2019; Submitted to the Physical Review Letters

Journal ref: Phys. Rev. Lett. 123, 051101 (2019)

Deep ugrizY Imaging and DEEP2/3 Spectroscopy: A Photometric Redshift Testbed for LSST and Public Release of Data from the DEEP3 Galaxy Redshift Survey

Authors: Rongpu Zhou, Michael C. Cooper, Jeffrey A. Newman, Matthew L. N. Ashby, James Aird, Christopher J. Conselice, Marc Davis, Aaron A. Dutton, S. M. Faber, Jerome J. Fang, G. G. Fazio, Puragra Guhathakurta, Dale Kocevski, David C. Koo, Kirpal Nandra, Andrew C. Phillips, David J. Rosario, Edward F. Schlafly, Jonathan R. Trump, Benjamin Weiner, Christopher N. A. Willmer, Renbin Yan

Abstract: We present catalogs of calibrated photometry and spectroscopic redshifts in the Extended Groth Strip, intended for studies of photometric redshifts (photo-z's). The data includes ugriz photometry from CFHTLS and Y-band photometry from the Subaru Suprime camera, as well as spectroscopic redshifts from the DEEP2, DEEP3 and 3D-HST surveys. These catalogs incorporate corrections to produce effectively… ▽ More We present catalogs of calibrated photometry and spectroscopic redshifts in the Extended Groth Strip, intended for studies of photometric redshifts (photo-z's). The data includes ugriz photometry from CFHTLS and Y-band photometry from the Subaru Suprime camera, as well as spectroscopic redshifts from the DEEP2, DEEP3 and 3D-HST surveys. These catalogs incorporate corrections to produce effectively matched-aperture photometry across all bands, based upon object size information available in the catalog and Moffat profile point spread function fits. We test this catalog with a simple machine learning-based photometric redshift algorithm based upon Random Forest regression, and find that the corrected aperture photometry leads to significant improvement in photo-z accuracy compared to the original SExtractor catalogs from CFHTLS and Subaru. The deep ugrizY photometry and spectroscopic redshifts are well-suited for empirical tests of photometric redshift algorithms for LSST. The resulting catalogs are publicly available. We include a basic summary of the strategy of the DEEP3 Galaxy Redshift Survey to accompany the recent public release of DEEP3 data. △ Less

Submitted 12 August, 2019; v1 submitted 19 March, 2019; originally announced March 2019.

Comments: 20 pages, 9 figures; published in MNRAS; for associated catalogs, see http://d-scholarship.pitt.edu/36064/

Journal ref: MNRAS 488, 4565-4584 (2019)

A gamma-ray periodic modulation in Globular Cluster 47 Tucanae

Authors: Peng-Fei Zhang, Jia-Neng Zhou, Da-Hai Yan, Jing-Zhi Yan, Yi-Zhong Fan, Jun Fang, Li Zhang

Abstract: The Globular Cluster 47 Tucanae was firstly detected in gamma-rays by the Large Area Telescope (LAT) onboard the \emph{Fermi} Gamma-ray Space Telescope, and the gamma-ray emission has been widely attributed to the millisecond pulsars. In this work, we analyze the Fermi-LAT pass 8 data ranging from 2008 August to 2017 May and report the detection of a modulation with a period of $18.416\pm0.008$ ho… ▽ More The Globular Cluster 47 Tucanae was firstly detected in gamma-rays by the Large Area Telescope (LAT) onboard the \emph{Fermi} Gamma-ray Space Telescope, and the gamma-ray emission has been widely attributed to the millisecond pulsars. In this work, we analyze the Fermi-LAT pass 8 data ranging from 2008 August to 2017 May and report the detection of a modulation with a period of $18.416\pm0.008$ hours at a significance level of $\sim 4.8σ$. This is the first time to detect a significant modulation with a period much longer than that of millisecond pulsars in gamma-rays from Globular Clusters. The periodic modulation signal appears in the {\it Swift}-BAT data as well. The phase-folded Chandra X-ray light curve of a point source may have provided an additional clue. △ Less

Submitted 2 December, 2020; v1 submitted 24 October, 2018; originally announced October 2018.

Comments: 7 pages, 6 figures. Comments welcome

Journal ref: ApJL, 904, L29 (2020)

Stability of the neutron-proton-electron matter under strong magnetic fields: the covariant Vlasov approach

Authors: Sidney Avancini, B. P. Bertolino, Aziz Rabhi, Jianjun Fang, Helena Pais, Constança Providência

Abstract: The neutron-proton-electron(npe) matter under a strong magnetic field is studied in the context of the covariant Vlasov approach. We use a Walecka-type hadronic model and the dispersion relations for the longitudinal and transverse modes are obtained. The instability regions for longitudinal and transverse modes are also studied. The crust-core transition of a magnetized neutron star is discussed. The neutron-proton-electron(npe) matter under a strong magnetic field is studied in the context of the covariant Vlasov approach. We use a Walecka-type hadronic model and the dispersion relations for the longitudinal and transverse modes are obtained. The instability regions for longitudinal and transverse modes are also studied. The crust-core transition of a magnetized neutron star is discussed. △ Less

Submitted 31 August, 2018; originally announced September 2018.

Comments: 14pages, 4 figures, accepted for publication in PRC

Journal ref: Phys. Rev. C 98, 025805, 2018

The evolution of galaxy shapes in CANDELS: from prolate to oblate

Authors: Haowen Zhang, Joel R. Primack, S. M. Faber, David C. Koo, Avishai Dekel, Zhu Chen, Daniel Ceverino, Yu-Yen Chang, Jerome J. Fang, Yicheng Guo, Lin Lin, Arjen van der Wel

Abstract: We model the projected b/a-log a distributions of CANDELS main sequence star-forming galaxies, where a (b) is the semi-major (semi-minor) axis of the galaxy images. We find that smaller-a galaxies are rounder at all stellar masses M and redshifts, so we include a when analyzing b/a distributions. Approximating intrinsic shapes of the galaxies as triaxial ellipsoids and assuming a multivariate norm… ▽ More We model the projected b/a-log a distributions of CANDELS main sequence star-forming galaxies, where a (b) is the semi-major (semi-minor) axis of the galaxy images. We find that smaller-a galaxies are rounder at all stellar masses M and redshifts, so we include a when analyzing b/a distributions. Approximating intrinsic shapes of the galaxies as triaxial ellipsoids and assuming a multivariate normal distribution of galaxy size and two shape parameters, we construct their intrinsic shape and size distributions to obtain the fractions of prolate, oblate and spheroidal galaxies in each redshift and mass bin. We find that galaxies tend to be prolate at low m and high redshifts, and oblate at high M and low redshifts, qualitatively consistent with van der Wel et al. (2014), implying that galaxies tend to evolve from prolate to oblate. These results are consistent with the predictions from simulations (Ceverino et al. 2015, Tomassetti et al. 2016) that the transition from prolate to oblate is caused by a compaction event at a characteristic mass range, making the galaxy center baryon dominated. We give probabilities of a galaxy's being prolate, oblate or spheroidal as a function of its M, redshift, projected b/a and a, which can facilitate target selections of galaxies with specific shapes at hight redshifts. We also give predicted optical depths of galaxies, which are qualitatively consistent with the expected correlation that AV should be higher for edge-on disk galaxies in each log a slice at low redshift and high mass bins. △ Less

Submitted 31 May, 2018; originally announced May 2018.

Comments: 24 pages, 25 figures, submitted to MNRAS

On the Transition of the Galaxy Quenching Mode at 0.5<z<1 in CANDELS

Authors: F. S. Liu, Meng Jia, Hassen M. Yesuf, S. M. Faber, David C. Koo, Yicheng Guo, Eric F. Bell, Dongfei Jiang, Weichen Wang, Anton M. Koekemoer, Xianzhong Zheng, Jerome J. Fang, Guillermo Barro, Pablo G. Pérez-González, Avishai Dekel, Dale Kocevski, Nimish P. Hathi, Darren Croton, M. Huertas-Company, Xianmin Meng, Wei Tong, Lu Liu

Abstract: We investigate the galaxy quenching process at intermediate redshift using a sample of $\sim4400$ galaxies with $M_{\ast} > 10^{9}M_{\odot}$ between redshift 0.5 and 1.0 in all five CANDELS fields. We divide this sample, using the integrated specific star formation rate (sSFR), into four sub-groups: star-forming galaxies (SFGs) above and below the ridge of the star-forming main sequence (SFMS), tr… ▽ More We investigate the galaxy quenching process at intermediate redshift using a sample of $\sim4400$ galaxies with $M_{\ast} > 10^{9}M_{\odot}$ between redshift 0.5 and 1.0 in all five CANDELS fields. We divide this sample, using the integrated specific star formation rate (sSFR), into four sub-groups: star-forming galaxies (SFGs) above and below the ridge of the star-forming main sequence (SFMS), transition galaxies and quiescent galaxies. We study their $UVI$ ($U-V$ versus $V-I$) color gradients to infer their sSFR gradients out to twice effective radii. We show that on average both star-forming and transition galaxies at all masses are not fully quenched at any radii, whereas quiescent galaxies are fully quenched at all radii. We find that at low masses ($M_{\ast} = 10^{9}-10^{10}M_{\odot}$) SFGs both above and below the SFMS ridge generally have flat sSFR profiles, whereas the transition galaxies at the same masses generally have sSFRs that are more suppressed in their outskirts. In contrast, at high masses ($M_{\ast} > 10^{10.5}M_{\odot}$), SFGs above and below the SFMS ridge and transition galaxies generally have varying degrees of more centrally-suppressed sSFRs relative to their outskirts. These findings indicate that at $z\sim~0.5-1.0$ the main galaxy quenching mode depends on its already formed stellar mass, exhibiting a transition from "the outside-in" at $M_{\ast} \leq 10^{10}M_{\odot}$ to "the inside-out" at $M_{\ast} > 10^{10.5}M_{\odot}$. In other words, our findings support that internal processes dominate the quenching of massive galaxies, whereas external processes dominate the quenching of low-mass galaxies. △ Less

Submitted 1 May, 2018; originally announced May 2018.

Comments: 19 pages, 12 figures, 2 tables, accepted for publication in ApJ