Showing 1–50 of 127 results for author: Ma, H

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

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

Density-based clustering algorithm for galaxy group/cluster identification

Authors: Hai-Xia Ma, Tsutomu T. Takeuchi, Suchetha Cooray, Yongda Zhu

Abstract: A direct approach to studying the galaxy-halo connection is the analysis of observed groups and clusters of galaxies that trace the underlying dark matter halos, making identifying galaxy clusters and their associated brightest cluster galaxies (BCGs) crucial. We test and propose a robust density-based clustering algorithm that outperforms the traditional Friends-of-Friends (FoF) algorithm in the… ▽ More A direct approach to studying the galaxy-halo connection is the analysis of observed groups and clusters of galaxies that trace the underlying dark matter halos, making identifying galaxy clusters and their associated brightest cluster galaxies (BCGs) crucial. We test and propose a robust density-based clustering algorithm that outperforms the traditional Friends-of-Friends (FoF) algorithm in the currently available galaxy group/cluster catalogs. Our new approach is a modified version of the Ordering Points To Identify the Clustering Structure (OPTICS) algorithm, which accounts for line-of-sight positional uncertainties due to redshift space distortions by incorporating a scaling factor, and is thereby referred to as sOPTICS. When tested on both a galaxy group catalog based on semi-analytic galaxy formation simulations and observational data, our algorithm demonstrated robustness to outliers and relative insensitivity to hyperparameter choices. In total, we compared the results of eight clustering algorithms. The proposed density-based clustering method, sOPTICS, outperforms FoF in accurately identifying giant galaxy clusters and their associated BCGs in various environments with higher purity and recovery rate, also successfully recovering 115 BCGs out of 118 reliable BCGs from a large galaxy sample. △ Less

Submitted 16 May, 2024; originally announced May 2024.

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

The evolutionary pathways of disk galaxies with different sizes

Authors: Hong-Chuan Ma, Min Du, Luis C. Ho, Ming-jie Sheng, Shihong Liao

Abstract: From the IllustrisTNG-50 simulation, a sample of 836 central disk galaxies with tiny stellar halos is chosen to study the inherent evolution of galaxies driven by nature. These galaxies are classified as compact, normal, or extended by referencing their locations on the mass-size ($M_\star-R_{\rm 1/2}$) diagram. This research demonstrates the distinctive evolutionary pathways of galaxies with diff… ▽ More From the IllustrisTNG-50 simulation, a sample of 836 central disk galaxies with tiny stellar halos is chosen to study the inherent evolution of galaxies driven by nature. These galaxies are classified as compact, normal, or extended by referencing their locations on the mass-size ($M_\star-R_{\rm 1/2}$) diagram. This research demonstrates the distinctive evolutionary pathways of galaxies with different sizes in IllustrisTNG simulations, primarily driven by nature. It is confirmed that disk galaxies inherit the angular momentum of their parent dark matter halos. More compact galaxies form earlier within halos possessing lower specific angular momentum through heightened star formation during the early phase at redshifts above 2. During the later phase, the size of extended galaxies experiences more pronounced growth by accreting gas with high angular momentum. Additionally, we reveal that many key characteristics of galaxies are linked to their mass and size: (1) compact galaxies tend to exhibit higher metal content, proportional to the potential well $\frac{M_\star}{R_{\rm 1/2}}$, (2) compact galaxies host more massive bulges and black holes, and higher central concentration. Furthermore, our analysis indicates that galaxies of all types continue to actively engage in star formation, with no evident signs of quenching attributed to their varying sizes and angular momenta. △ Less

Submitted 2 July, 2024; v1 submitted 16 April, 2024; originally announced April 2024.

Comments: 11 pages, 11 figures. Accepted for publication in A&A

Disentanglement of mixed interference fringes in optical interferometers: theory and applications

Authors: Kaiyuan Yang, Weilong Wei, Xiafei Ma, Botao Chen, Junqiu Chu, Xinling Liu, Yuhua Cheng, Hu Yang, Haotong Ma, Bo Qi, Zongliang Xie

Abstract: Optical interferometric imaging enables astronomical observation at extremely high angular resolution. The necessary optical information for imaging, such as the optical path differences and visibilities, is easy to extract from fringes generated by the combination of two beams. With more than two apertures, the image-plane interference pattern becomes an increasingly indistinguishable mixture of… ▽ More Optical interferometric imaging enables astronomical observation at extremely high angular resolution. The necessary optical information for imaging, such as the optical path differences and visibilities, is easy to extract from fringes generated by the combination of two beams. With more than two apertures, the image-plane interference pattern becomes an increasingly indistinguishable mixture of fringe spacings and directions. For decades, the state-of-the-art approaches for obtaining two-aperture fringes from an interferometer array composed of many apertures are limited to pairwise combinations using bulk optics. Here, we derive and demonstrate a fringe disentanglement theory that can digitally transform the interference pattern of N apertures to N(N-1)/2 pairwise fringes without any optics, thus providing straightforward methods of information acquisition for interferometers. We demonstrate applications of our technique by both simulation and experiment, showing that this theory can be used for simultaneously sensing pistons and determining the individual visibilities of all combining apertures. Furthermore, we use the proposed theory to phase a 1.5-meter segmented flat telescope, demonstrating its validity for engineering implementation. This theory may not only benefit optical imaging but also interferometry-based measurements, by providing an exceptional capability to simplify the interferometric output generated by a system of many apertures. △ Less

Submitted 9 April, 2024; originally announced April 2024.

Comments: 21 pages, 9 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.

The Physical Origin of the Mass-Size Relation and Its Scatter of Disk Galaxies

Authors: Min Du, Hong-Chuan Ma, Wen-Yu Zhong, Luis C. Ho, Shihong Liao, Yingjie Peng

Abstract: Utilizing a kinematic decomposition of simulated galaxies, we focus on galaxies with tiny kinematically inferred stellar halos, indicative of weak external influences. We investigate the intricate interplay between internal (natural) and external (nurture) processes in shaping the scaling relationships of specific angular momentum ($j_\star$), stellar mass ($M_\star$), and size of disk galaxies wi… ▽ More Utilizing a kinematic decomposition of simulated galaxies, we focus on galaxies with tiny kinematically inferred stellar halos, indicative of weak external influences. We investigate the intricate interplay between internal (natural) and external (nurture) processes in shaping the scaling relationships of specific angular momentum ($j_\star$), stellar mass ($M_\star$), and size of disk galaxies within the IllustrisTNG simulation. The correlation among mass, size, and angular momentum of galaxies is examined by comparing simulations with observations and the theoretical predictions of the exponential hypothesis. Galaxies with tiny stellar halos exhibit a large scatter in the $j_\star$-$M_\star$ relation, which suggests that it is inherently present in their initial conditions. The analysis reveals that the disks of these galaxies adhere to the exponential hypothesis, resulting in a tight fiducial $j_\star$-$M_\star$-scale length (size) relation that is qualitatively consistent with observations. The inherent scatter in $j_\star$ provides a robust explanation for the mass-size relation and its substantial variability. Notably, galaxies that are moderately influenced by external processes closely adhere to a scaling relation akin to that of galaxies with tiny stellar halos. This result underscores the dominant role of internal processes in shaping the overall $j_\star$-$M_\star$ and mass-size relation, with external effects playing a relatively minor role in disk galaxies. Furthermore, the correlation between galaxy size and the virial radius of the dark matter halo exists but fails to provide strong evidence of the connection between galaxies and their parent dark matter halos. △ Less

Submitted 25 March, 2024; originally announced March 2024.

Comments: 11 pages, 9 figures, accepted for publication in A&A

Journal ref: A&A 686, A168 (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)

Spin speed correlations and the evolution of galaxy-halo systems

Authors: Ming-Jie Sheng, Lin Zhu, Hao-Ran Yu, Hongchuan Ma, Haikun Li, Peng Wang, Xi Kang

Abstract: Galaxy angular momenta (spins) contain valuable cosmological information, complementing their positions and velocities. The baryonic spin direction of galaxies has been probed as a reliable tracer of their host halos and the primordial spin modes. Here we use the TNG100 simulation of the IllustrisTNG project to study the spin magnitude correlations between dark matter, gas, and stellar components… ▽ More Galaxy angular momenta (spins) contain valuable cosmological information, complementing their positions and velocities. The baryonic spin direction of galaxies has been probed as a reliable tracer of their host halos and the primordial spin modes. Here we use the TNG100 simulation of the IllustrisTNG project to study the spin magnitude correlations between dark matter, gas, and stellar components of galaxy-halo systems and their evolutions across cosmic history. We find that these components generate similar initial spin magnitudes from the same tidal torque in Lagrangian space. At low redshifts, the gas component still traces the spin magnitude of the dark matter halo and the primordial spin magnitude. However, the traceability of the stellar component depends on the $ex$ $situ$ stellar mass fraction, $f_{\rm acc}$. Our results suggest that the galaxy baryonic spin magnitude can also serve as a tracer of their host halo and the initial perturbations, and the galaxy-halo correlations are affected by the similarity of their evolution histories. △ Less

Submitted 2 July, 2024; v1 submitted 14 November, 2023; originally announced November 2023.

Comments: 10 pages, 10 figures. Matches the accepted version in Physical Review D

Journal ref: Physical Review D, 109, 123548, 2024

Cloud-cloud collision in S235: triggered the formation of high-mass stars and young star clusters

Authors: En Chen, Yu Gao, Shiyu Zhang, Xuepeng Chen, Min Fang, Qianru He, Xuejian Jiang, Yan Sun, Xiaolong Wang, Hongjun Ma

Abstract: We present the analysis of cloud-cloud collision (CCC) process in the Galactic molecular complex S235. Our new CO observations performed with the PMO-13.7m telescope reveal two molecular clouds, namely the S235-Main and the S235-ABC, with $\sim$ 4 km s$^{-1}$ velocity separation. The bridge feature, the possible colliding interface and the complementary distribution of the two clouds are significa… ▽ More We present the analysis of cloud-cloud collision (CCC) process in the Galactic molecular complex S235. Our new CO observations performed with the PMO-13.7m telescope reveal two molecular clouds, namely the S235-Main and the S235-ABC, with $\sim$ 4 km s$^{-1}$ velocity separation. The bridge feature, the possible colliding interface and the complementary distribution of the two clouds are significant observational signatures of cloud-cloud collision in S235. The most direct evidence of cloud-cloud collision process in S235 is that the S235-Main (in a distance of 1547$^{+44}_{-43}$ pc) and S235-ABC (1567$^{+33}_{-39}$ pc) meet at almost the same position (within 1$σ$ error range) at a supersonic relative speed. We identified ten $^{13}$CO clumps from PMO-13.7m observations, 22 dust cores from the archival SCUBA-2 data, and 550 YSOs from NIR-MIR data. 63$\%$ of total YSOs are clustering in seven MST groups (M1$-$M7). The tight association between the YSO groups (M1 $\&$ M7) and the bridge feature suggests that the CCC process triggers star formation there. The collisional impact subregion (the South) shows $3\sim5$ times higher CFE and SFE (average value of 12.3$\%$ and 10.6$\%$, respectively) than the non-collisional impact subregion (2.4$\%$ and 2.6$\%$, respectively), suggesting that the CCC process may have enhanced the CFE and SFE of the clouds compared to those without collision influence. △ Less

Submitted 7 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

3D hydrodynamics simulations of a 3 $M_{\odot}$ core-helium burning star

Authors: Simon Blouin, Falk Herwig, Huaqing Mao, Pavel Denissenkov, Paul R. Woodward

Abstract: The inner structure of core-helium burning (CHeB) stars remains uncertain due to the yet unknown nature of mixing at the boundary of their cores. Large convective cores beyond a bare Schwarzschild model are favoured both from theoretical arguments and from asteroseismological constraints. However, the exact nature of this extra mixing, and in particular the possible presence of semiconvective laye… ▽ More The inner structure of core-helium burning (CHeB) stars remains uncertain due to the yet unknown nature of mixing at the boundary of their cores. Large convective cores beyond a bare Schwarzschild model are favoured both from theoretical arguments and from asteroseismological constraints. However, the exact nature of this extra mixing, and in particular the possible presence of semiconvective layers, is still debated. In this work, we approach this problem through a new avenue by performing the first full-sphere 3D hydrodynamics simulations of the interiors of CHeB stars. We use the PPMstar explicit gas dynamics code to simulate the inner 0.45 $M_{\odot}$ of a 3 $M_{\odot}$ CHeB star. Simulations are performed using different Cartesian grid resolutions (768$^3$, 1152$^3$ and 1728$^3$) and heating rates. We use two different initial states, one based on MESA's predictive mixing scheme (which significantly extends the core beyond the Schwarzschild boundary) and one based on the convective premixing approach (which exhibits a semiconvective interface). The general behaviour of the flow in the convective core and in the stable envelope (where internal gravity waves are observed) is consistent with our recent simulations of core convection in massive main-sequence stars, and so are the various luminosity scaling relations. The semiconvective layers are dominated by strong internal gravity waves that do not produce measurable species mixing, but overshooting motions from the convective core gradually homogenize the semiconvective interface. This process can possibly completely erase the semiconvective layers, which would imply that CHeB stars do not harbour a semiconvection zone. △ Less

Submitted 17 November, 2023; v1 submitted 21 August, 2023; originally announced August 2023.

Comments: Accepted for publication in MNRAS

Probing Ultra-late Reionization: Direct Measurements of the Mean Free Path over $5<z<6$

Authors: Yongda Zhu, George D. Becker, Holly M. Christenson, Anson D'Aloisio, Sarah E. I. Bosman, Tom Bakx, Valentina D'Odorico, Manuela Bischetti, Christopher Cain, Frederick B. Davies, Rebecca L. Davies, Anna-Christina Eilers, Xiaohui Fan, Prakash Gaikwad, Martin G. Haehnelt, Laura C. Keating, Girish Kulkarni, Samuel Lai, Hai-Xia Ma, Andrei Mesinger, Yuxiang Qin, Sindhu Satyavolu, Tsutomu T. Takeuchi, Hideki Umehata, Jinyi Yang

Abstract: The mean free path of ionizing photons, $λ_{\rm mfp}$, is a critical parameter for modeling the intergalactic medium (IGM) both during and after reionization. We present direct measurements of $λ_{\rm mfp}$ from QSO spectra over the redshift range $5<z<6$, including the first measurements at $z\simeq5.3$ and 5.6. Our sample includes data from the XQR-30 VLT large program, as well as new Keck/ESI o… ▽ More The mean free path of ionizing photons, $λ_{\rm mfp}$, is a critical parameter for modeling the intergalactic medium (IGM) both during and after reionization. We present direct measurements of $λ_{\rm mfp}$ from QSO spectra over the redshift range $5<z<6$, including the first measurements at $z\simeq5.3$ and 5.6. Our sample includes data from the XQR-30 VLT large program, as well as new Keck/ESI observations of QSOs near $z \sim 5.5$, for which we also acquire new [C II] 158$μ$m redshifts with ALMA. By measuring the Lyman continuum transmission profile in stacked QSO spectra, we find $λ_{\rm mfp} = 9.33_{-1.80}^{+2.06}$, $5.40_{-1.40}^{+1.47}$, $3.31_{-1.34}^{+2.74}$, and $0.81_{-0.48}^{+0.73}$ pMpc at $z=5.08$, 5.31, 5.65, and 5.93, respectively. Our results demonstrate that $λ_{\rm mfp}$ increases steadily and rapidly with time over $5<z<6$. Notably, we find that $λ_{\rm mfp}$ deviates significantly from predictions based on a fully ionized and relaxed IGM as late as $z=5.3$. By comparing our results to model predictions and indirect $λ_{\rm mfp}$ constraints based on IGM Ly$α$ opacity, we find that the $λ_{\rm mfp}$ evolution is consistent with scenarios wherein the IGM is still undergoing reionization and/or retains large fluctuations in the ionizing UV background well below redshift six. △ Less

Submitted 8 August, 2023; originally announced August 2023.

Comments: 17 pages, 6 figures, 2 tables; accepted for publication in ApJ

Wave generation and energetic electron scattering in solar flares

Authors: Hanqing Ma, James F. Drake, Marc Swisdak

Abstract: We conduct two-dimensional particle-in-cell simulations to investigate the scattering of electron heat flux by self-generated oblique electromagnetic waves. The heat flux is modeled as a bi-kappa distribution with a T_parallel > T_perp temperature anisotropy maintained by continuous injection at the boundaries. The anisotropic distribution excites oblique whistler waves and filamentary-like Weibel… ▽ More We conduct two-dimensional particle-in-cell simulations to investigate the scattering of electron heat flux by self-generated oblique electromagnetic waves. The heat flux is modeled as a bi-kappa distribution with a T_parallel > T_perp temperature anisotropy maintained by continuous injection at the boundaries. The anisotropic distribution excites oblique whistler waves and filamentary-like Weibel instabilities. Electron velocity distributions taken after the system has reached a steady state show that these in stabilities inhibit the heat flux and drive the total distributions towards isotropy. Electron trajectories in velocity space show a circular-like diffusion along constant energy surfaces in the wave frame. The key parameter controlling the scattering rate is the average speed, or drift speed vd, of the heat flux compared with the electron Alfven speed vAe, with higher drift speeds producing stronger fluctua tions and a more significant reduction of the heat flux. Reducing the density of the electrons carrying the heat flux by 50% does not significantly affect the scattering rate. A scaling law for the electron scattering rate versus vd/vAe is deduced from the simulations. The implications of these results for understanding energetic electron transport during solar flare energy release are discussed. △ Less

Submitted 6 July, 2023; originally announced July 2023.

On the detection of the electromagnetic counterparts from lensed gravitational wave events by binary neutron star mergers

Authors: Hao Ma, Youjun Lu, Xiao Guo, Siqi Zhang, Qingbo Chu

Abstract: Future ground-based gravitational wave (GW) detectors, i.e., Einstein telescope (ET) and Cosmic Explorer (CE), are expected to detect a significant number of lensed binary neutron star (BNS) mergers, which may provide a unique tool to probe cosmology. In this paper, we investigate the detectability of the optical/infrared electromagnetic (EM) counterparts (kilonovae/afterglows) from these lensed B… ▽ More Future ground-based gravitational wave (GW) detectors, i.e., Einstein telescope (ET) and Cosmic Explorer (CE), are expected to detect a significant number of lensed binary neutron star (BNS) mergers, which may provide a unique tool to probe cosmology. In this paper, we investigate the detectability of the optical/infrared electromagnetic (EM) counterparts (kilonovae/afterglows) from these lensed BNS mergers by future GW detectors and EM telescopes using simple kilonova, afterglow, and lens models. ET and CE are expected to detect $\sim5.32^{+26.1}_{-5.10}$ and $67.3^{+332}_{-64.7}$ lensed BNS mergers per year. We find that the EM counterparts associated with all these mergers will be detectable by an all sky-survey in the H-band with the limiting magnitude $m_{\textrm{lim}}\gtrsim27$, while the detectable fraction is $\lesssim0.4\%$ in the g-/z-band if with $m_{\textrm{lim}}\lesssim24$. Generally it is more efficient to search the lensed EM counterparts by adopting the infrared bands than the optical/UV bands with the same $m_{\textrm{lim}}$. Future telescopes like Vera C. Rubin Observatory, China Space Station Telescope, and Euclid can hardly detect the EM counterparts of even one lensed BNS merger. Roman Space Telescope (RST) and James Webb Space Telescope (JWST) have the capability to detect about a few or more such events per year. Moreover, the time delays and separations between the lensed image pairs are typically in the ranges from minutes to months and from $0.1$ to $1$\,arcsec, suggesting that both the GW and EM images of most lensed BNS mergers can be well resolved by not only CE/ET in the time domain but also RST/JWST spatially. △ Less

Submitted 5 July, 2023; originally announced July 2023.

Comments: 16 pages, 10 figures, MNRAS in press

Diffractive lensing of nano-Hertz gravitational waves emitted from supermassive binary black holes by intervening galaxies

Authors: Hao Ma, Youjun Lu, Zhiwei Chen, Yunfeng Chen

Abstract: Pulsar timing array (PTA) experiments are expected to detect nano-Hertz gravitational waves (GWs) emitted from individual inspiralling supermassive binary black holes (SMBBHs). The GW signals from a small fraction of these SMBBHs may be diffractively lensed by intervening galaxies. In this paper, we investigate the diffractive lensing effects on the continuous GW signals from the lensed SMBBHs and… ▽ More Pulsar timing array (PTA) experiments are expected to detect nano-Hertz gravitational waves (GWs) emitted from individual inspiralling supermassive binary black holes (SMBBHs). The GW signals from a small fraction of these SMBBHs may be diffractively lensed by intervening galaxies. In this paper, we investigate the diffractive lensing effects on the continuous GW signals from the lensed SMBBHs and estimate the detectable number of such signals by PTAs, such as the Chinese PTA (CPTA) and the Square Kilometer Array (SKA) PTA. We find that the amplitude of the lensed GW signals may be only amplified by a factor of $\sim 1.01-1.14$ ($16\%-84\%$ range) and the phase of the signals may shift somewhat due to the lensing, significantly different from those strongly lensed high frequency GW signals from compact binary mergers in the geometric optics. We estimate that $\sim 0.01\%$ of all detected nano-Hertz GW signals from individual SMBBHs by future PTA experiments are lensed by foreground galaxies (i.e., up to $\sim 106$ for CPTA and up to $\sim 289$ for SKA-PTA). However, the lensed nano-Hertz GW signals are difficult to be distinguished from those without lensing by the PTA observations only. We further discuss the possibility about the identification of the lensed nano-Hertz GW signals from SMBBHs via the electromagnetic detection of their host galaxies or active galactic nuclei. △ Less

Submitted 5 July, 2023; originally announced July 2023.

Comments: 12 pages, 8 figures, Accepted for publication in MNRAS

Origin of superionic state in Earth's inner core

Authors: Ina Park, Yu He, Ho-kwang Mao, Ji Hoon Shim, Duck Young Kim

Abstract: Earth's inner core (IC) serves as a reservoir for volatile elements, which significantly affects its behavior and properties. Recent studies suggest that superionicity can be observed in ice and iron hydrides under high-pressure and temperature conditions, providing an alternative understanding of the planet's interior. In this study, we demonstrated that electride formation drives the superionic… ▽ More Earth's inner core (IC) serves as a reservoir for volatile elements, which significantly affects its behavior and properties. Recent studies suggest that superionicity can be observed in ice and iron hydrides under high-pressure and temperature conditions, providing an alternative understanding of the planet's interior. In this study, we demonstrated that electride formation drives the superionic state in iron hydride under IC pressure conditions. The electride stabilizes the iron lattice and provides a pathway for volatile diffusion. The coupling between lattice stability and superionicity is triggered near 100 GPa and enhanced at higher pressures. The electride-driven superionicity can also be generalized for volatiles in other rocky planetary cores. These findings provide new insights into the mechanisms of core formation and evolution of rocky planets. △ Less

Submitted 13 June, 2023; originally announced June 2023.

Comments: 26 pages, 4 main figures and 6 supplementary 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

The Lobster Eye Imager for Astronomy Onboard the SATech-01 Satellite

Authors: Z. X. Ling, X. J. Sun, C. Zhang, S. L. Sun, G. Jin, S. N. Zhang, X. F. Zhang, J. B. Chang, F. S. Chen, Y. F. Chen, Z. W. Cheng, W. Fu, Y. X. Han, H. Li, J. F. Li, Y. Li, Z. D. Li, P. R. Liu, Y. H. Lv, X. H. Ma, Y. J. Tang, C. B. Wang, R. J. Xie, Y. L. Xue, A. L. Yan , et al. (101 additional authors not shown)

Abstract: The Lobster Eye Imager for Astronomy (LEIA), a pathfinder of the Wide-field X-ray Telescope of the Einstein Probe (EP) mission, was successfully launched onboard the SATech-01 satellite of the Chinese Academy of Sciences on 27 July 2022. In this paper, we introduce the design and on-ground test results of the LEIA instrument. Using state-of-the-art Micro-Pore Optics (MPO), a wide field-of-view (Fo… ▽ More The Lobster Eye Imager for Astronomy (LEIA), a pathfinder of the Wide-field X-ray Telescope of the Einstein Probe (EP) mission, was successfully launched onboard the SATech-01 satellite of the Chinese Academy of Sciences on 27 July 2022. In this paper, we introduce the design and on-ground test results of the LEIA instrument. Using state-of-the-art Micro-Pore Optics (MPO), a wide field-of-view (FoV) of 346 square degrees (18.6 degrees * 18.6 degrees) of the X-ray imager is realized. An optical assembly composed of 36 MPO chips is used to focus incident X-ray photons, and four large-format complementary metal-oxide semiconductor (CMOS) sensors, each of 6 cm * 6 cm, are used as the focal plane detectors. The instrument has an angular resolution of 4 - 8 arcmin (in FWHM) for the central focal spot of the point spread function, and an effective area of 2 - 3 cm2 at 1 keV in essentially all the directions within the field of view. The detection passband is 0.5 - 4 keV in the soft X-rays and the sensitivity is 2 - 3 * 10-11 erg s-1 cm-2 (about 1 mini-Crab) at 1,000 second observation. The total weight of LEIA is 56 kg and the power is 85 W. The satellite, with a design lifetime of 2 years, operates in a Sun-synchronous orbit of 500 km with an orbital period of 95 minutes. LEIA is paving the way for future missions by verifying in flight the technologies of both novel focusing imaging optics and CMOS sensors for X-ray observation, and by optimizing the working setups of the instrumental parameters. In addition, LEIA is able to carry out scientific observations to find new transients and to monitor known sources in the soft X-ray band, albeit limited useful observing time available. △ Less

Submitted 24 May, 2023; originally announced May 2023.

Comments: Accepted by RAA

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)

3D hydrodynamic simulations of massive main-sequence stars. III. The effect of radiation pressure and diffusion leading to a 1D equilibrium model

Authors: Huaqing Mao, Paul Woodward, Falk Herwig, Pavel A. Denissenkov, Simon Blouin, William Thompson

Abstract: We present 3D hydrodynamical simulations of core convection with a stably stratified envelope of a 25 $\mathrm{M}_\odot$ star in the early phase of the main-sequence. We use the explicit gas-dynamics code $\texttt{PPMstar}$ which tracks two fluids and includes radiation pressure and radiative diffusion. Multiple series of simulations with different luminosities and radiative thermal conductivities… ▽ More We present 3D hydrodynamical simulations of core convection with a stably stratified envelope of a 25 $\mathrm{M}_\odot$ star in the early phase of the main-sequence. We use the explicit gas-dynamics code $\texttt{PPMstar}$ which tracks two fluids and includes radiation pressure and radiative diffusion. Multiple series of simulations with different luminosities and radiative thermal conductivities are presented. The entrainment rate at the convective boundary, internal gravity waves in and above the boundary region, and the approach to dynamical equilibrium shortly after a few convective turnovers are investigated. From the results of these simulations we extrapolate to find the entrainment rate at the nominal heating rate and thermal diffusion given by the $\texttt{MESA}$ stellar evolution model on which the 3D stratification is based. Further, to study the effect of radiative diffusion on the thermal timescale, we perform very long simulations accelerated by 10000 times their nominal luminosities. In these simulations the growing penetrative convection reduces the initially unrealistically large entrainment. This reduction is enabled by a spatial separation that develops between the entropy gradient and the composition gradient. The convective boundary moves outward much more slowly at the end of these simulations. Finally, we present a method to predict the extent and character of penetrative convection beyond the Schwarzschild boundary. This method is intended to be ultimately deployed in 1D stellar evolution calculations and is based on the properties of penetrative convection in our simulations carried forward through the local thermal timescale. △ Less

Submitted 20 April, 2023; originally announced April 2023.

Comments: 33 pages, 31 figures; submitted to ApJ

3D hydrodynamics simulations of core convection in supermassive main-sequence stars

Authors: Simon Blouin, Huaqing Mao, Tyrone E. Woods, Pavel Denissenkov, Paul R. Woodward, Falk Herwig

Abstract: Supermassive stars are Population III stars with masses exceeding $10^4\,M_{\odot}$ that could be the progenitors of the first supermassive black holes. Their interiors are in a regime where radiation pressure dominates the equation of state. In this work, we use the explicit gas dynamics code PPMstar to simulate the hydrogen-burning core of a $10^4\,M_{\odot}$ supermassive main-sequence star. The… ▽ More Supermassive stars are Population III stars with masses exceeding $10^4\,M_{\odot}$ that could be the progenitors of the first supermassive black holes. Their interiors are in a regime where radiation pressure dominates the equation of state. In this work, we use the explicit gas dynamics code PPMstar to simulate the hydrogen-burning core of a $10^4\,M_{\odot}$ supermassive main-sequence star. These are the first 3D hydrodynamics simulations of core convection in supermassive stars. We perform a series of ten simulations at different heating rates and on Cartesian grids with resolutions of $768^3$, $1152^3$ and $1728^3$. We examine different properties of the convective flow, including its large-scale morphology, its velocity spectrum and its mixing properties. We conclude that the radiation pressure-dominated nature of the interior does not noticeably affect the behaviour of convection compared to the case of core convection in a massive main-sequence star where gas pressure dominates. Our simulations also offer support for the use of mixing-length theory in 1D models of supermassive stars. △ Less

Submitted 17 March, 2023; originally announced March 2023.

Comments: 9 pages, 9 figures, movies at https://www.ppmstar.org/, accepted for publication in MNRAS

Measurement of the Compton scattering in germanium with a p-type point-contact germanium detector for dark matter detection

Authors: J. W. Hu, L. T. Yang, Q. Yue, X. P. Geng, H. B. Li, Y. F. Liang, S. T. Lin, S. K. Liu, H. Ma, L. Wang, H. T. Wong, B. T. Zhang

Abstract: Low-energy background through Compton scattering from the ambient $γ$ rays can be contaminated in direct dark matter search experiments. In this paper, we report comparable measurements on low-energy spectra via Compton scattering from several $γ$-ray sources with a p-type point-contact germanium detector. The spectra between 500 eV and 18 keV have been measured and analyzed. Moreover, the feature… ▽ More Low-energy background through Compton scattering from the ambient $γ$ rays can be contaminated in direct dark matter search experiments. In this paper, we report comparable measurements on low-energy spectra via Compton scattering from several $γ$-ray sources with a p-type point-contact germanium detector. The spectra between 500 eV and 18 keV have been measured and analyzed. Moreover, the features of the electron binding effect, particularly at the edges of the K- and L-shells in the germanium atom, were observed with different gamma sources and were consistent with the models in the Geant4 simulation. An empirical background model is proposed that provides insights into understanding the low-energy background in germanium for direct dark matter experiments. △ Less

Submitted 14 March, 2023; originally announced March 2023.

Comments: 7 pages, 5 figures and 2 tables; version accepted by JINST

3D hydrodynamics simulations of internal gravity waves in red giant branch stars

Authors: Simon Blouin, Huaqing Mao, Falk Herwig, Pavel Denissenkov, Paul R. Woodward, William R. Thompson

Abstract: We present the first 3D hydrodynamics simulations of the excitation and propagation of internal gravity waves (IGWs) in the radiative interiors of low-mass stars on the red giant branch (RGB). We use the PPMstar explicit gas dynamics code to simulate a portion of the convective envelope and all the radiative zone down to the hydrogen-burning shell of a 1.2$M_{\odot}$ upper RGB star. We perform sim… ▽ More We present the first 3D hydrodynamics simulations of the excitation and propagation of internal gravity waves (IGWs) in the radiative interiors of low-mass stars on the red giant branch (RGB). We use the PPMstar explicit gas dynamics code to simulate a portion of the convective envelope and all the radiative zone down to the hydrogen-burning shell of a 1.2$M_{\odot}$ upper RGB star. We perform simulations for different grid resolutions (768$^3$, 1536$^3$ and 2880$^3$), a range of driving luminosities, and two different stratifications (corresponding to the bump luminosity and the tip of the RGB). Our RGB tip simulations can be directly performed at the nominal luminosity, circumventing the need for extrapolations to lower luminosities. A rich, continuous spectrum of IGWs is observed, with a significant amount of total power contained at high wavenumbers. By following the time evolution of a passive dye in the stable layers, we find that IGW mixing in our simulations is weaker than predicted by a simple analytical prescription based on shear mixing and not efficient enough to explain the missing RGB extra mixing. However, we may be underestimating the efficiency of IGW mixing given that our simulations include a limited portion of the convective envelope. Quadrupling its radial extent compared to our fiducial setup increases convective velocities by up to a factor 2 and IGW velocities by up to a factor 4. We also report the formation of a $\sim 0.2\,H_P$ penetration zone and evidence that IGWs are excited by plumes that overshoot into the stable layers. △ Less

Submitted 13 March, 2023; originally announced March 2023.

Comments: 20 pages, 29 figures, movies at https://www.ppmstar.org, resubmitted to MNRAS

3D hydrodynamic simulations of massive main-sequence stars II. Convective excitation and spectra of internal gravity waves

Authors: William Thompson, Falk Herwig, Paul R. Woodward, Huaqing Mao, Pavel Denissenkov, Dominic M. Bowman, Simon Blouin

Abstract: Recent photometric observations of massive stars have identified a low-frequency power excess which appears as stochastic low-frequency variability in light curve observations. We present the oscillation properties of high resolution hydrodynamic simulations of a 25 $\mathrm{M}_\odot$ star performed with the PPMStar code. The model star has a convective core mass of… ▽ More Recent photometric observations of massive stars have identified a low-frequency power excess which appears as stochastic low-frequency variability in light curve observations. We present the oscillation properties of high resolution hydrodynamic simulations of a 25 $\mathrm{M}_\odot$ star performed with the PPMStar code. The model star has a convective core mass of $\approx\, 12\, \mathrm{M}_\odot$ and approximately half of the envelope simulated. From this simulation, we extract light curves from several directions, average them over each hemisphere, and process them as if they were real photometric observations. We show how core convection excites waves with a similar frequency as the convective time scale in addition to significant power across a forest of low and high angular degree $l$ modes. We find that the coherence of these modes is relatively low as a result of their stochastic excitation by core convection, with lifetimes on the order of 10s of days. Thanks to the still significant power at higher $l$ and this relatively low coherence, we find that integrating over a hemisphere produces a power spectrum that still contains measurable power up to the Brunt--Väisälä frequency. These power spectra extracted from the stable envelope are qualitatively similar to observations, with same order of magnitude yet lower characteristic frequency. This work further shows the potential of long-duration, high-resolution hydrodynamic simulations for connecting asteroseismic observations to the structure and dynamics of core convection and the convective boundary. △ Less

Submitted 29 April, 2024; v1 submitted 10 March, 2023; originally announced March 2023.

Comments: Accepted version (2024 April 26). Originally submitted 2022 November 04

3D hydrodynamic simulations of massive main-sequence stars. I. Dynamics and mixing of convection and internal gravity waves

Authors: Falk Herwig, Paul R. Woodward, Huaqing Mao, William R. Thompson, Pavel Denissenkov, Josh Lau, Simon Blouin, Robert Andrassy, Adam Paul

Abstract: We performed 3D hydrodynamic simulations of the inner $\approx 50\%$ radial extent of a $25\ \mathrm{M_\odot}$ star in the early phase of the main sequence and investigate core convection and internal gravity waves in the core-envelope boundary region. Simulations for different grid resolutions and driving luminosities establish scaling relations to constrain models of mixing for 1D applications.… ▽ More We performed 3D hydrodynamic simulations of the inner $\approx 50\%$ radial extent of a $25\ \mathrm{M_\odot}$ star in the early phase of the main sequence and investigate core convection and internal gravity waves in the core-envelope boundary region. Simulations for different grid resolutions and driving luminosities establish scaling relations to constrain models of mixing for 1D applications. As in previous works, the turbulent mass entrainment rate extrapolated to nominal heating is unrealistically high ($1.58\times 10^{-4}\ \mathrm{M_\odot/yr}$), which is discussed in terms of the non-equilibrium response of the simulations to the initial stratification. We measure quantitatively the effect of mixing due to internal gravity waves excited by core convection interacting with the boundary in our simulations. The wave power spectral density as a function of frequency and wavelength agrees well with the GYRE eigenmode predictions based on the 1D spherically averaged radial profile. A diffusion coefficient profile that reproduces the spherically averaged abundance distribution evolution is determined for each simulation. Through a combination of eigenmode analysis and scaling relations it is shown that in the $N^2$-peak region, mixing is due to internal gravity waves and follows the scaling relation $D_\mathrm{IGW-hydro} \propto L^{4/3}$ over a $\gtrapprox 2\ \mathrm{dex}$ range of heating factors. Different extrapolations of the mixing efficiency down to nominal heating are discussed. If internal gravity wave mixing is due to thermally-enhanced shear mixing, an upper limit is $D_\mathrm{IGW} \lessapprox 2$ to $3\times 10^4\ \mathrm{cm^2/s}$ at nominal heating in the $N^2$-peak region above the convective core. △ Less

Submitted 19 July, 2023; v1 submitted 9 March, 2023; originally announced March 2023.

Comments: MNRAS. Clarifications and improvements following referee recommendations. Movies and data access https://www.ppmstar.org. 38 Figs. In v2 the bottom panel of Fig 28 is showing the wrong case (M116 instead of M114). Corrected in v3

How Close Dark Matter Halos and MOND Are to Each Other: Three-Dimensional Tests Based on Gaia DR2

Authors: Yongda Zhu, Hai-Xia Ma, Xiao-Bo Dong, Yang Huang, Tobias Mistele, Bo Peng, Qian Long, Tianqi Wang, Liang Chang, Xi Jin

Abstract: Aiming at discriminating different gravitational potential models of the Milky Way, we perform tests based on the kinematic data powered by the Gaia DR2 astrometry, over a large range of $(R,z)$ locations. Invoking the complete form of Jeans equations that admit three integrals of motion, we use the independent $R$- and $z$-directional equations as two discriminators ($T_R$ and $T_z$). We apply th… ▽ More Aiming at discriminating different gravitational potential models of the Milky Way, we perform tests based on the kinematic data powered by the Gaia DR2 astrometry, over a large range of $(R,z)$ locations. Invoking the complete form of Jeans equations that admit three integrals of motion, we use the independent $R$- and $z$-directional equations as two discriminators ($T_R$ and $T_z$). We apply the formula for spatial distributions of radial and vertical velocity dispersions proposed by Binney et al., and successfully extend it to azimuthal components, $σ_θ(R,z)$ and $V_θ(R,z)$; the analytic form avoids the numerical artifacts caused by numerical differentiation in Jeans-equations calculation given the limited spatial resolutions of observations, and more importantly reduces the impact of kinematic substructures in the Galactic disk. It turns out that whereas the current kinematic data are able to reject Moffat's Modified Gravity (let alone the Newtonian baryon-only model), Milgrom's MOND is still not rejected. In fact, both the carefully calibrated fiducial model invoking a spherical dark matter (DM) halo and MOND are equally consistent with the data at almost all spatial locations (except that probably both have respective problems at low-$|z|$ locations), no matter which a tracer population or which a meaningful density profile is used. Because there is no free parameter at all in the quasi-linear MOND model we use, and the baryonic parameters are actually fine-tuned in the DM context, such an effective equivalence is surprising, and might be calling forth a transcending synthesis of the two paradigms. △ Less

Submitted 23 November, 2022; originally announced November 2022.

Comments: 20 pages, 13 figures; accepted for publication in MNRAS

First wide field-of-view X-ray observations by a lobster eye focusing telescope in orbit

Authors: C. Zhang, Z. X. Ling, X. J. Sun, S. L. Sun, Y. Liu, Z. D. Li, Y. L. Xue, Y. F. Chen, Y. F. Dai, Z. Q. Jia, H. Y. Liu, X. F. Zhang, Y. H. Zhang, S. N. Zhang, F. S. Chen, Z. W. Cheng, W. Fu, Y. X. Han, H. Li, J. F. Li, Y. Li, P. R. Liu, X. H. Ma, Y. J. Tang, C. B. Wang , et al. (53 additional authors not shown)

Abstract: As a novel X-ray focusing technology, lobster eye micro-pore optics (MPO) feature both a wide observing field of view and true imaging capability, promising sky monitoring with significantly improved sensitivity and spatial resolution in soft X-rays. Since first proposed by Angel (1979), the optics have been extensively studied, developed and trialed over the past decades. In this Letter, we repor… ▽ More As a novel X-ray focusing technology, lobster eye micro-pore optics (MPO) feature both a wide observing field of view and true imaging capability, promising sky monitoring with significantly improved sensitivity and spatial resolution in soft X-rays. Since first proposed by Angel (1979), the optics have been extensively studied, developed and trialed over the past decades. In this Letter, we report on the first-light results from a flight experiment of the Lobster Eye Imager for Astronomy ($LEIA$), a pathfinder of the wide-field X-ray telescope of the Einstein Probe mission. The piggyback imager, launched in July 2022, has a mostly un-vignetted field of view of $18.6^\circ \times 18.6^\circ $. Its spatial resolution is in the range of 4$-$7 arcmin in FWHM and the focal spot effective area is 2$-$3 cm$^2$, both showing only mild fluctuations across the field of view. We present images of the Galactic center region, Sco X-1 and the diffuse Cygnus Loop nebular taken in snapshot observations over 0.5$-$4 keV. These are truly wide-field X-ray images of celestial bodies observed, for the first time, by a focusing imaging telescope. Initial analyses of the in-flight data show excellent agreement between the observed images and the on-ground calibration and simulations. The instrument and its characterization are briefly described, as well as the flight experiment. The results provide a solid basis for the development of the present and proposed wide-field X-ray missions using lobster eye MPO. △ Less

Submitted 17 November, 2022; originally announced November 2022.

Comments: 11 pages, 4 figures. Accepted for publication in Astrophysical Journal Letter

Galaxy Manifold: Characterizing and understanding galaxies with two parameters

Authors: Suchetha Cooray, Tsutomu T. Takeuchi, Daichi Kashino, Shuntaro A. Yoshida, Hai-Xia Ma, Kai T. Kono

Abstract: We report the discovery of a two-dimensional Galaxy Manifold within the multi-dimensional luminosity space of local galaxies. The multi-dimensional luminosity space is constructed using 11 bands that span from far ultraviolet to near-infrared for redshift < 0.1 galaxies observed with GALEX, SDSS, and UKIDSS. The two latent parameters are sufficient to express 93.2% of the variance in the galaxy sa… ▽ More We report the discovery of a two-dimensional Galaxy Manifold within the multi-dimensional luminosity space of local galaxies. The multi-dimensional luminosity space is constructed using 11 bands that span from far ultraviolet to near-infrared for redshift < 0.1 galaxies observed with GALEX, SDSS, and UKIDSS. The two latent parameters are sufficient to express 93.2% of the variance in the galaxy sample, suggesting that this Galaxy Manifold is one of the most efficient representations of galaxies. The transformation between the observed luminosities and the manifold parameters as an analytic mapping is provided. The manifold representation provides accurate (85%) morphological classifications with a simple linear boundary, and galaxy properties can be estimated with minimal scatter (0.12 dex and 0.04 dex for star formation rate and stellar mass, respectively) by calibrating with the two-dimensional manifold location. Under the assumption that the manifold expresses the possible parameter space of galaxies, the evolution on the manifold is considered. We find that constant and exponentially decreasing star formation histories form almost orthogonal modes of evolution on the manifold. Through these simple models, we understand that the two modes are closely related to gas content, which suggests the close relationship of the manifold to gas accretion. Without assuming a star formation history, a gas-regulated model reproduces an exponentially declining star formation history with a timescale of $\sim$1.2 Gyrs on the manifold. Lastly, the found manifold suggests a paradigm where galaxies are characterized by their mass/scale and specific SFR, which agrees with previous studies of dimensionality reduction. △ Less

Submitted 11 October, 2022; originally announced October 2022.

Comments: Submitted to MNRAS. 15 pages and 17 figures. All comments are welcome

Dark Matter: DAMA/LIBRA and its perspectives

Authors: R. Bernabei, P. Belli, F. Cappella, V. Caracciolo, R. Cerulli, C. J. Dai, A. d'Angelo, A. Incicchitti, A. Leoncini, X. H. Ma, V. Merlo, F. Montecchia, X. D. Sheng, Z. P. Ye

Abstract: The long-standing model-independent annual modulation effect measured by DAMA deep underground at Gran Sasso Laboratory with different experimental configurations is summarized and perspectives will be highlighted. DAMA/LIBRA-phase2 set-up, $\simeq$ 250 kg highly radio-pure NaI(Tl) confirms the evidence of a signal that meets all the requirements of the model independent Dark Matter annual modulat… ▽ More The long-standing model-independent annual modulation effect measured by DAMA deep underground at Gran Sasso Laboratory with different experimental configurations is summarized and perspectives will be highlighted. DAMA/LIBRA-phase2 set-up, $\simeq$ 250 kg highly radio-pure NaI(Tl) confirms the evidence of a signal that meets all the requirements of the model independent Dark Matter annual modulation signature at high C.L.; the full exposure is 2.86 ton $\times$ yr over 22 annual cycles. The experiment is currently collecting data in the DAMA/LIBRA-phase2 empowered configuration with an even lower software energy threshold. Other recent claims are shortly commented. △ Less

Submitted 2 September, 2022; originally announced September 2022.

Comments: Proceedings of IDM 2022, July 18-22, 2022. A Section is dedicated to rebut the arguments of arXiv:2208.05158

Distributed YSOs in the Perseus Molecular Cloud from the Gaia and LAMOST Surveys

Authors: Xiao-Long Wang, Min Fang, Yu Gao, Hong-Xin Zhang, Gregory J. Herczeg, Hong-Jun Ma, En Chen, Xing-Yu Zhou

Abstract: Identifying the young optically visible population in a star-forming region is essential for fully understanding the star formation event. In this paper, We identify 211 candidate members of the Perseus molecular cloud based on Gaia astronomy. We use LAMOST spectra to confirm that 51 of these candidates are new members, bringing the total census of known members to 856. The newly confirmed members… ▽ More Identifying the young optically visible population in a star-forming region is essential for fully understanding the star formation event. In this paper, We identify 211 candidate members of the Perseus molecular cloud based on Gaia astronomy. We use LAMOST spectra to confirm that 51 of these candidates are new members, bringing the total census of known members to 856. The newly confirmed members are less extincted than previously known members. Two new stellar aggregates are identified in our updated census. With the updated member list, we obtain a statistically significant distance gradient of $\rm 4.84\;pc\;deg^{-1}$ from west to east. Distances and extinction corrected color-magnitude diagrams indicate that NGC 1333 is significantly younger than IC 348 and the remaining cloud regions. The disk fraction in NGC 1333 is higher than elsewhere, consistent with its youngest age. The star formation scenario in the Perseus molecular cloud is investigated and the bulk motion of the distributed population is consistent with the cloud being swept away by the Per-Tau Shell. △ Less

Submitted 27 July, 2022; originally announced July 2022.

Comments: 27 pages, 14 figures, 1 table; Accepted for publication in ApJ

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)

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

Dynamics in a stellar convective layer and at its boundary: Comparison of five 3D hydrodynamics codes

Authors: R. Andrassy, J. Higl, H. Mao, M. Mocák, D. G. Vlaykov, W. D. Arnett, I. Baraffe, S. W. Campbell, T. Constantino, P. V. F. Edelmann, T. Goffrey, T. Guillet, F. Herwig, R. Hirschi, L. Horst, G. Leidi, C. Meakin, J. Pratt, F. Rizzuti, F. K. Roepke, P. Woodward

Abstract: Our ability to predict the structure and evolution of stars is in part limited by complex, 3D hydrodynamic processes such as convective boundary mixing. Hydrodynamic simulations help us understand the dynamics of stellar convection and convective boundaries. However, the codes used to compute such simulations are usually tested on extremely simple problems and the reliability and reproducibility o… ▽ More Our ability to predict the structure and evolution of stars is in part limited by complex, 3D hydrodynamic processes such as convective boundary mixing. Hydrodynamic simulations help us understand the dynamics of stellar convection and convective boundaries. However, the codes used to compute such simulations are usually tested on extremely simple problems and the reliability and reproducibility of their predictions for turbulent flows is unclear. We define a test problem involving turbulent convection in a plane-parallel box, which leads to mass entrainment from, and internal-wave generation in, a stably stratified layer. We compare the outputs from the codes FLASH, MUSIC, PPMSTAR, PROMPI, and SLH, which have been widely employed to study hydrodynamic problems in stellar interiors. The convection is dominated by the largest scales that fit into the simulation box. All time-averaged profiles of velocity components, fluctuation amplitudes, and fluxes of enthalpy and kinetic energy are within $\lesssim 3σ$ of the mean of all simulations on a given grid ($128^3$ and $256^3$ grid cells), where $σ$ describes the statistical variation due to the flow's time dependence. They also agree well with a $512^3$ reference run. The $128^3$ and $256^3$ simulations agree within $9\%$ and $4\%$, respectively, on the total mass entrained into the convective layer. The entrainment rate appears to be set by the amount of energy that can be converted to work in our setup and details of the small-scale flows in the boundary layer seem to be largely irrelevant. Our results lend credence to hydrodynamic simulations of flows in stellar interiors. We provide in electronic form all outputs of our simulations as well as all information needed to reproduce or extend our study. △ Less

Submitted 26 January, 2022; v1 submitted 1 November, 2021; originally announced November 2021.

Comments: 20 pages, 19 figures, new Fig. B.2, new clarifications, electronic materials updated, accepted for publication in A&A

Journal ref: A&A 659, A193 (2022)

The dark matter: DAMA/LIBRA and its perspectives

Authors: R. Bernabei, P. Belli, V. Caracciolo, R. Cerulli, V. Merlo, F. Cappella, A. d'Angelo, A. Incicchitti, C. J. Dai, X. H. Ma, X. D. Sheng, F. Montecchia, Z. P. Ye

Abstract: Experimental observations and theoretical arguments point out that Dark Matter (DM) particles are one of the most prominent component of the Universe. This motivated the pioneer DAMA experiment to investigate the presence of these particles in the galactic halo, by exploiting the model independent signature of the DM annual modulation of the rate and very highly radio-pure apparatus in underground… ▽ More Experimental observations and theoretical arguments point out that Dark Matter (DM) particles are one of the most prominent component of the Universe. This motivated the pioneer DAMA experiment to investigate the presence of these particles in the galactic halo, by exploiting the model independent signature of the DM annual modulation of the rate and very highly radio-pure apparatus in underground site. In this paper the results obtained by other two annual cycles of DAMA/LIBRA-phase2 are presented and the long-standing model-independent annual modulation effect measured by DAMA deep underground at the Gran Sasso National Laboratory (LNGS) of the I.N.F.N. with different experimental configurations is summarized. The improved experimental configuration of DAMA/LIBRA-phase2, $\simeq$ 250 kg highly radio-pure NaI(Tl), allowed to lower the software energy threshold. The total exposure of DAMA/LIBRA-phase2 over 8 annual cycles is 1.53 ton $\times$ yr. DAMA/LIBRA-phase2 confirms the evidence of a signal that meets all the requirements of the model independent Dark Matter annual modulation signature, at 11.8 $σ$ C.L. in the energy region (1-6) keV. In the energy region between 2 and 6 keV, where data are also available from DAMA/NaI and DAMA/LIBRA-phase1 (2.86 ton $\times$ yr), the achieved C.L. is 13.7 $σ$; the modulation amplitude of the single-hit scintillation events is: $(0.01014 \pm 0.00074)$ cpd/kg/keV, the measured phase is $(142.4 \pm 4.2)$ days and the measured period is $(0.99834 \pm 0.00067)$ yr, all these values are well in agreement with those expected for DM particles. No systematics or side reaction able to mimic the exploited DM signature (i.e. to account for the whole measured modulation amplitude and to simultaneously satisfy all the requirements of the signature), has been found or suggested by anyone throughout some decades thus far. △ Less

Submitted 10 October, 2021; originally announced October 2021.

Comments: Proceedings of the Sixteenth Marcel Grossmann Meeting - MG16, July 5-10, 2021. arXiv admin note: substantial text overlap with arXiv:1805.10486, arXiv:1308.5109, arXiv:1002.1028, arXiv:1306.1411

Alternative Sample Mass Measurement Technique for OSIRIS-REX Sample Collection Phase

Authors: Huikang Ma, Michael Skeen, Ryan Olds, Brennen Miller, Dante S. Lauretta

Abstract: The Origins, Spectral Interpretation, Resource Identification, and Security-Regolith Explorer (OSIRIS-REx) spacecraft is the third NASA New Frontiers Program mission and arrived at the near-Earth asteroid (101955) Bennu in December 2018. Following completion of sample collection in October 2020, otherwise known as Touch-And-Go (TAG), the OSIRIS-REx spacecraft was set to verify its collected sample… ▽ More The Origins, Spectral Interpretation, Resource Identification, and Security-Regolith Explorer (OSIRIS-REx) spacecraft is the third NASA New Frontiers Program mission and arrived at the near-Earth asteroid (101955) Bennu in December 2018. Following completion of sample collection in October 2020, otherwise known as Touch-And-Go (TAG), the OSIRIS-REx spacecraft was set to verify its collected sample mass requirement (> 60g of material). The thoroughly tested Sample Mass Measurement (SMM) method was to be used for this verification. Imaging of the Touch-And-Go Sample Acquisition Mechanism (TAGSAM) was received shortly following the TAG event, intended to ensure mechanism health prior to moving forward with the SMM activity. These images displayed sample leakage, prompting discussion for alternative paths forward. Risk of continued sample loss and a desire to retain as much material as possible lead the team to pursue an accelerated sample stow schedule and forgo the planned SMM activity. Once the sample was safely stowed in the return capsule an alternative SMM method was proposed. The alternative SMM technique utilized reaction wheel momentum data from identical TAGSAM movements prior to and following the TAG event to estimate changes in spacecraft moment of inertia. Conservation of momentum was used to isolate the sample mass from this inertia change. Using this new method, the spacecraft team was able to successfully estimate collected sample mass to be 250.37 +/- 101 g. △ Less

Submitted 12 September, 2021; originally announced September 2021.

Comments: Paper to be submitted to the 44th Annual AAS Guidance, Navigation and Control (GN&C) Conference

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)

CHIMPS2: Survey description and $^{12}$CO emission in the Galactic Centre

Authors: D. J. Eden, T. J. T. Moore, M. J. Currie, A. J. Rigby, E. Rosolowsky, Y. Su, Kee-Tae Kim, H. Parsons, O. Morata, H. -R. Chen, T. Minamidani, Geumsook Park, S. E. Ragan, J. S. Urquhart, R. Rani, K. Tahani, S. J. Billington, S. Deb, C. Figura, T. Fujiyoshi, G. Joncas, L. W. Liao, T. Liu, H. Ma, P. Tuan-Anh , et al. (81 additional authors not shown)

Abstract: The latest generation of Galactic-plane surveys is enhancing our ability to study the effects of galactic environment upon the process of star formation. We present the first data from CO Heterodyne Inner Milky Way Plane Survey 2 (CHIMPS2). CHIMPS2 is a survey that will observe the Inner Galaxy, the Central Molecular Zone (CMZ), and a section of the Outer Galaxy in $^{12}$CO, $^{13}$CO, and C… ▽ More The latest generation of Galactic-plane surveys is enhancing our ability to study the effects of galactic environment upon the process of star formation. We present the first data from CO Heterodyne Inner Milky Way Plane Survey 2 (CHIMPS2). CHIMPS2 is a survey that will observe the Inner Galaxy, the Central Molecular Zone (CMZ), and a section of the Outer Galaxy in $^{12}$CO, $^{13}$CO, and C$^{18}$O $(J = 3\rightarrow2)$ emission with the Heterodyne Array Receiver Program on the James Clerk Maxwell Telescope (JCMT). The first CHIMPS2 data presented here are a first look towards the CMZ in $^{12}$CO J = 3$\rightarrow$2 and cover $-3^{\circ}\leq\,\ell\,\leq\,5^{\circ}$ and $\mid$b$\mid \leq 0.5^{\circ}$ with angular resolution of 15 arcsec, velocity resolution of 1 km s$^{-1}$, and rms $ΔT_A ^\ast =$ 0.58 K at these resolutions. Such high-resolution observations of the CMZ will be a valuable data set for future studies, whilst complementing the existing Galactic Plane surveys, such as SEDIGISM, the Herschel infrared Galactic Plane Survey, and ATLASGAL. In this paper, we discuss the survey plan, the current observations and data, as well as presenting position-position maps of the region. The position-velocity maps detect foreground spiral arms in both absorption and emission. △ Less

Submitted 10 September, 2020; originally announced September 2020.

Comments: 18 pages, 14 figures, accepted for publication in MNRAS

3D1D hydro-nucleosynthesis simulations. I. Advective-reactive post-processing method and its application to H ingestion into He-shell flash convection in rapidly accreting white dwarfs

Authors: David Stephens, Falk Herwig, Paul Woodward, Pavel Denissenkov, Robert Andrassy, Huaqing Mao

Abstract: We present two mixing models for post-processing of 3D hydrodynamic simulations applied to convective-reactive i-process nucleosynthesis in a rapidly accreting white dwarf (RAWD) with $\mathrm{[Fe/H]} = -2.6$, in which H is ingested into a convective He shell. A 1D advective two-stream model adopts physically motivated radial and horizontal mixing coefficients constrained by 3D hydrodynamic simula… ▽ More We present two mixing models for post-processing of 3D hydrodynamic simulations applied to convective-reactive i-process nucleosynthesis in a rapidly accreting white dwarf (RAWD) with $\mathrm{[Fe/H]} = -2.6$, in which H is ingested into a convective He shell. A 1D advective two-stream model adopts physically motivated radial and horizontal mixing coefficients constrained by 3D hydrodynamic simulations. A simpler approach uses diffusion coefficients calculated from the same simulations. All 3D simulations include the energy feedback of the $^{12}$C(p,$γ$)$^{13}$N reaction from the H entrainment. Global oscillations of shell H ingestion in two of the RAWD simulations cause bursts of entrainment of H and non-radial hydrodynamic feedback. With the same nuclear network as in the 3D simulations, the 1D advective two-stream model reproduces the rate and location of the H burning within the He shell closely matching the 3D simulation predictions, as well as qualitatively displaying the asymmetry of the $X_{\mathrm{H}}$ profiles between the up- and downstream. With a full i-process network the advective mixing model captures the difference in the n-capture nucleosynthesis in the up- and downstream. For example, $^{89}$Kr and $^{90}$Kr with half-lives of 3.18 min and 32.3 s differ by a factor 2-10 in the two streams. In this particular application the diffusion approach provides globally the same abundance distribution as the advective two-stream mixing model. The resulting i-process yields are in excellent agreement with observations of the exemplary CEMP-r/s star CS31062-050. △ Less

Submitted 17 February, 2021; v1 submitted 29 January, 2020; originally announced January 2020.

Comments: 18 pages, 17 figures, 2 tables, accepted for publication in MNRAS

Non-Gaussianity in a general kind noncanonical warm inflation

Authors: Xiao-Min Zhang, Kai Li, Hong-Yang Ma, Qian Liu, Jian-Yang Zhu

Abstract: The general kind noncanonical warm inflation is introduced and its total non-Gaussianity of perturbation on uniform energy density hypersurfaces is studied. The total non-Gaussianity contains two complementary parts: the three- and the four-point correlations. The three-point correlation part denoted by $f^{int}_{NL}$, comes from the three-point correlation of inflaton field, while the four-point… ▽ More The general kind noncanonical warm inflation is introduced and its total non-Gaussianity of perturbation on uniform energy density hypersurfaces is studied. The total non-Gaussianity contains two complementary parts: the three- and the four-point correlations. The three-point correlation part denoted by $f^{int}_{NL}$, comes from the three-point correlation of inflaton field, while the four-point correlation part, denoted by $f_{NL}^{δN}$, is the contribution due to the four-point correlation function of the inflaton field. The aforementioned two parts of non-Gaussianity in general noncanonical warm inflation are calculated and analysed respectively, and the comparisons and discussions of the two parts are finally carried out in this study. △ Less

Submitted 6 October, 2019; v1 submitted 3 October, 2019; originally announced October 2019.

Comments: 8pages,0figures

Journal ref: Phys. Rev. D 100, 123521 (2019)

Improved model-dependent corollary analyses after the first six annual cycles of DAMA/LIBRA-phase2

Authors: R. Bernabei, P. Belli, F. Cappella, V. Caracciolo, R. Cerulli, C. J. Dai, A. d'Angelo, A. Di Marco, H. L. He, A. Incicchitti, X. H. Ma, V. Merlo, F. Montecchia, X. D. Sheng, Z. P. Ye

Abstract: Several of the many proposed Dark Matter candidate particles, already investigated with lower exposure and a higher software energy threshold, are further analyzed including the first DAMA/LIBRA--phase2 data release, with an exposure of 1.13 ton $\times$ yr and a lower software energy threshold (1 keV). The cumulative exposure above 2 keV considering also DAMA/NaI and DAMA/LIBRA--phase1 results is… ▽ More Several of the many proposed Dark Matter candidate particles, already investigated with lower exposure and a higher software energy threshold, are further analyzed including the first DAMA/LIBRA--phase2 data release, with an exposure of 1.13 ton $\times$ yr and a lower software energy threshold (1 keV). The cumulative exposure above 2 keV considering also DAMA/NaI and DAMA/LIBRA--phase1 results is now 2.46 ton $\times$ yr. The analysis permits to constraint the parameters' space of the considered candidates restricting their values -- with respect to previous analyses -- thanks to the increase of the exposure and to the lower energy threshold. △ Less

Submitted 23 December, 2019; v1 submitted 15 July, 2019; originally announced July 2019.

Comments: 48 pages, 33 figures; in press on Nucl. Phys. At. Energy (2019)

Journal ref: Nucl. Phys. At. Energy 20 (2019) 317-348

Simulating 3-D Stellar Hydrodynamics using PPM and PPB Multifluid Gas Dynamics on CPU and CPU+GPU Nodes

Authors: Paul R. Woodward, Pei-Hung Lin, Huaqing Mao, Robert Andrassy, Falk Herwig

Abstract: The special computational challenges of simulating 3-D hydrodynamics in deep stellar interiors are discussed, and numerical algorithmic responses described. Results of recent simulations carried out at scale on the NSF's Blue Waters machine at the University of Illinois are presented, with a special focus on the computational challenges they address. Prospects for future work using GPU-accelerated… ▽ More The special computational challenges of simulating 3-D hydrodynamics in deep stellar interiors are discussed, and numerical algorithmic responses described. Results of recent simulations carried out at scale on the NSF's Blue Waters machine at the University of Illinois are presented, with a special focus on the computational challenges they address. Prospects for future work using GPU-accelerated nodes such as those on the DoE's new Summit machine at Oak Ridge National Laboratory are described, with a focus on numerical algorithmic accommodations that we believe will be necessary. △ Less

Submitted 31 October, 2018; originally announced October 2018.

Comments: 14 pages, 7 figures

H-He Shell Interactions and Nucleosynthesis in Massive Population III Stars

Authors: Ondrea Clarkson, Falk Herwig, Robert Andrassy, Marco Pignatari, Paul Woodward, Huaqing Mao

Abstract: We report on our ongoing investigation into the nucleosynthetic and hydrodynamic nature of mixing at the interface between the H- and He-convection zones in massive Pop III stars. Studying a grid of 26 1D stellar evolution simulations with different mixing assumptions, we find that H-He interactions occur in 23/26 cases. We demonstrate the nucleosynthesis expected in a H-He interaction in an 80M… ▽ More We report on our ongoing investigation into the nucleosynthetic and hydrodynamic nature of mixing at the interface between the H- and He-convection zones in massive Pop III stars. Studying a grid of 26 1D stellar evolution simulations with different mixing assumptions, we find that H-He interactions occur in 23/26 cases. We demonstrate the nucleosynthesis expected in a H-He interaction in an 80M$_\odot$. Finally, we describe our progress in simulating a Pop III double convection zone in the PPMStar hydrodynamics code. △ Less

Submitted 29 October, 2018; originally announced October 2018.

Comments: NIC 2018 Proceeding

The Large Area Detector onboard the eXTP mission

Authors: Marco Feroci, Mahdi Ahangarianabhari, Giovanni Ambrosi, Filippo Ambrosino, Andrea Argan, Marco Barbera, Joerg Bayer, Pierluigi Bellutti, Bruna Bertucci, Giuseppe Bertuccio, Giacomo Borghi, Enrico Bozzo, Franck Cadoux, Riccardo Campana, Francesco Ceraudo, Tianxiang Chen, Daniela Cirrincione, Alessandra De Rosa, Ettore Del Monte, Sergio Di Cosimo, Sebastian Diebold, Yuri Evangelista, Qingmei Fan, Yannick Favre, Francesco Ficorella , et al. (46 additional authors not shown)

Abstract: The eXTP (enhanced X-ray Timing and Polarimetry) mission is a major project of the Chinese Academy of Sciences (CAS) and China National Space Administration (CNSA) currently performing an extended phase A study and proposed for a launch by 2025 in a low-earth orbit. The eXTP scientific payload envisages a suite of instruments (Spectroscopy Focusing Array, Polarimetry Focusing Array, Large Area Det… ▽ More The eXTP (enhanced X-ray Timing and Polarimetry) mission is a major project of the Chinese Academy of Sciences (CAS) and China National Space Administration (CNSA) currently performing an extended phase A study and proposed for a launch by 2025 in a low-earth orbit. The eXTP scientific payload envisages a suite of instruments (Spectroscopy Focusing Array, Polarimetry Focusing Array, Large Area Detector and Wide Field Monitor) offering unprecedented simultaneous wide-band X-ray spectral, timing and polarimetry sensitivity. A large European consortium is contributing to the eXTP study and it is expected to provide key hardware elements, including a Large Area Detector (LAD). The LAD instrument for eXTP is based on the design originally proposed for the LOFT mission within the ESA context. The eXTP/LAD envisages a deployed 3.4 m2 effective area in the 2-30 keV energy range, achieved through the technology of the large-area Silicon Drift Detectors - offering a spectral resolution of up to 200 eV FWHM at 6 keV - and of capillary plate collimators - limiting the field of view to about 1 degree. In this paper we provide an overview of the LAD instrument design, including new elements with respect to the earlier LOFT configuration. △ Less

Submitted 30 July, 2018; originally announced July 2018.

Comments: Paper presented at SPIE Space Telescopes and Instrumentation 2018

Journal ref: Proc. of SPIE Vol. 10699, 106991C - Space Telescopes and Instrumentation 2018: Ultraviolet to Gamma Ray, edited by Jan-Willem A. den Herder, Shouleh Nikzad, Kazuhiro Nakazawa

First model independent results from DAMA/LIBRA-phase2

Authors: R. Bernabei, P. Belli, A. Bussolotti, F. Cappella, V. Caracciolo, R. Cerulli, C. J. Dai, A. d'Angelo, A. Di Marco, H. L. He, A. Incicchitti, X. H. Ma, A. Mattei, V. Merlo, F. Montecchia, X. D. Sheng, Z. P. Ye

Abstract: The first model independent results obtained by the DAMA/LIBRA-phase2 experiment are presented. The data have been collected over 6 annual cycles corresponding to a total exposure of 1.13 ton $\times$ yr, deep underground at the Gran Sasso National Laboratory (LNGS) of the I.N.F.N. The DAMA/LIBRA-phase2 apparatus, $\simeq$ 250 kg highly radio-pure NaI(Tl), profits from a second generation high qua… ▽ More The first model independent results obtained by the DAMA/LIBRA-phase2 experiment are presented. The data have been collected over 6 annual cycles corresponding to a total exposure of 1.13 ton $\times$ yr, deep underground at the Gran Sasso National Laboratory (LNGS) of the I.N.F.N. The DAMA/LIBRA-phase2 apparatus, $\simeq$ 250 kg highly radio-pure NaI(Tl), profits from a second generation high quantum efficiency photomultipliers and of new electronics with respect to DAMA/LIBRA-phase1. The improved experimental configuration has also allowed to lower the software energy threshold. New data analysis strategies are presented. The DAMA/LIBRA-phase2 data confirm the evidence of a signal that meets all the requirements of the model independent Dark Matter (DM) annual modulation signature, at 9.5 $σ$ C.L. in the energy region (1-6) keV. In the energy region between 2 and 6 keV, where data are also available from DAMA/NaI and DAMA/LIBRA-phase1 (exposure $1.33$ ton $\times$ yr, collected over 14 annual cycles), the achieved C.L. for the full exposure (2.46 ton $\times$ yr) is 12.9 $σ$; the modulation amplitude of the single-hit scintillation events is: $(0.0103 \pm 0.0008)$ cpd/kg/keV, the measured phase is $(145 \pm 5)$ days and the measured period is $(0.999 \pm 0.001)$ yr, all these values are well in agreement with those expected for DM particles. No systematics or side reaction able to mimic the exploited DM signature (i.e. to account for the whole measured modulation amplitude and to simultaneously satisfy all the requirements of the signature), has been found or suggested by anyone throughout some decades thus far. △ Less

Submitted 11 January, 2019; v1 submitted 26 May, 2018; originally announced May 2018.

Comments: 28 pages, 17 figures, 5 tables. In press on Journal Nucl. Phys. and Atomic Energy. arXiv admin note: substantial text overlap with arXiv:1308.5109, arXiv:1002.1028, arXiv:1306.1411

Journal ref: Nucl. Phys. At. Energy 19 (2018) 307-325

Galactic Cosmic-Ray Anisotropy in the Northern hemisphere from the ARGO-YBJ Experiment during 2008-2012

Authors: B. Bartoli, P. Bernardini, X. J. Bi, Z. Cao, S. Catalanotti, S. Z. Chen, T. L. Chen, S. W. Cui, B. Z. Dai, A. D'Amone, Danzengluobu, I. De Mitri, B. D'Ettorre Piazzoli, T. Di Girolamo, G. Di Sciascio, C. F. Feng, Z. Y. Feng, W. Gao, Q. B. Gou, Y. Q. Guo, H. H. He, Haibing Hu, Hongbo Hu, M. Iacovacci, R. Iuppa , et al. (47 additional authors not shown)

Abstract: This paper reports on the observation of the sidereal large-scale anisotropy of cosmic rays using data collected by the ARGO-YBJ experiment over 5 years (2008$-$2012). This analysis extends previous work limited to the period from 2008 January to 2009 December,near the minimum of solar activity between cycles 23 and 24.With the new data sample the period of solar cycle 24 from near minimum to maxi… ▽ More This paper reports on the observation of the sidereal large-scale anisotropy of cosmic rays using data collected by the ARGO-YBJ experiment over 5 years (2008$-$2012). This analysis extends previous work limited to the period from 2008 January to 2009 December,near the minimum of solar activity between cycles 23 and 24.With the new data sample the period of solar cycle 24 from near minimum to maximum is investigated. A new method is used to improve the energy reconstruction, allowing us to cover a much wider energy range, from 4 to 520 TeV. Below 100 TeV, the anisotropy is dominated by two wide regions, the so-called "tail-in" and "loss-cone" features. At higher energies, a dramatic change of the morphology is confirmed. The yearly time dependence of the anisotropy is investigated. Finally, no noticeable variation of cosmic-ray anisotropy with solar activity is observed for a median energy of 7 TeV. △ Less

Submitted 23 May, 2018; originally announced May 2018.