Showing 1–50 of 234 results for author: Sun, X

Spectropolarimetric Inversion in Four Dimensions with Deep Learning (SPIn4D): I. Overview, Magnetohydrodynamic Modeling, and Stokes Profile Synthesis

Authors: Kai E. Yang, Lucas A. Tarr, Matthias Rempel, S. Curt Dodds, Sarah A. Jaeggli, Peter Sadowski, Thomas A. Schad, Ian Cunnyngham, Jiayi Liu, Yannik Glaser, Xudong Sun

Abstract: The National Science Foundation's Daniel K. Inouye Solar Telescope (DKIST) will provide high-resolution, multi-line spectropolarimetric observations that are poised to revolutionize our understanding of the Sun. Given the massive data volume, novel inference techniques are required to unlock its full potential. Here, we provide an overview of our "SPIn4D" project, which aims to develop deep convol… ▽ More The National Science Foundation's Daniel K. Inouye Solar Telescope (DKIST) will provide high-resolution, multi-line spectropolarimetric observations that are poised to revolutionize our understanding of the Sun. Given the massive data volume, novel inference techniques are required to unlock its full potential. Here, we provide an overview of our "SPIn4D" project, which aims to develop deep convolutional neural networks (CNNs) for estimating the physical properties of the solar photosphere from DKIST spectropolarimetric observations. We describe the magnetohydrodynamic (MHD) modeling and the Stokes profile synthesis pipeline that produce the simulated output and input data, respectively. These data will be used to train a set of CNNs that can rapidly infer the four-dimensional MHD state vectors by exploiting the spatiotemporally coherent patterns in the Stokes profile time series. Specifically, our radiative MHD model simulates the small-scale dynamo actions that are prevalent in quiet-Sun and plage regions. Six cases with different mean magnetic fields have been conducted; each case covers six solar-hours, totaling 109 TB in data volume. The simulation domain covers at least $25\times25\times8$ Mm with $16\times16\times12$ km spatial resolution, extending from the upper convection zone up to the temperature minimum region. The outputs are stored at a 40 s cadence. We forward model the Stokes profile of two sets of Fe I lines at 630 and 1565 nm, which will be simultaneously observed by DKIST and can better constrain the parameter variations along the line of sight. The MHD model output and the synthetic Stokes profiles are publicly available. △ Less

Submitted 29 July, 2024; originally announced July 2024.

Comments: Submitted to AAS Journals. Comments welcome

ANDES, the high resolution spectrograph for the ELT: science goals, project overview and future developments

Authors: A. Marconi, M. Abreu, V. Adibekyan, V. Alberti, S. Albrecht, J. Alcaniz, M. Aliverti, C. Allende Prieto, J. D. Alvarado Gómez, C. S. Alves, P. J. Amado, M. Amate, M. I. Andersen, S. Antoniucci, E. Artigau, C. Bailet, C. Baker, V. Baldini, A. Balestra, S. A. Barnes, F. Baron, S. C. C. Barros, S. M. Bauer, M. Beaulieu, O. Bellido-Tirado , et al. (264 additional authors not shown)

Abstract: The first generation of ELT instruments includes an optical-infrared high-resolution spectrograph, indicated as ELT-HIRES and recently christened ANDES (ArmazoNes high Dispersion Echelle Spectrograph). ANDES consists of three fibre-fed spectrographs ([U]BV, RIZ, YJH) providing a spectral resolution of $\sim$100,000 with a minimum simultaneous wavelength coverage of 0.4-1.8 $μ$m with the goal of ex… ▽ More The first generation of ELT instruments includes an optical-infrared high-resolution spectrograph, indicated as ELT-HIRES and recently christened ANDES (ArmazoNes high Dispersion Echelle Spectrograph). ANDES consists of three fibre-fed spectrographs ([U]BV, RIZ, YJH) providing a spectral resolution of $\sim$100,000 with a minimum simultaneous wavelength coverage of 0.4-1.8 $μ$m with the goal of extending it to 0.35-2.4 $μ$m with the addition of a U arm to the BV spectrograph and a separate K band spectrograph. It operates both in seeing- and diffraction-limited conditions and the fibre feeding allows several, interchangeable observing modes including a single conjugated adaptive optics module and a small diffraction-limited integral field unit in the NIR. Modularity and fibre-feeding allow ANDES to be placed partly on the ELT Nasmyth platform and partly in the Coudé room. ANDES has a wide range of groundbreaking science cases spanning nearly all areas of research in astrophysics and even fundamental physics. Among the top science cases, there are the detection of biosignatures from exoplanet atmospheres, finding the fingerprints of the first generation of stars, tests on the stability of Nature's fundamental couplings, and the direct detection of the cosmic acceleration. The ANDES project is carried forward by a large international consortium, composed of 35 Institutes from 13 countries, forming a team of almost 300 scientists and engineers which include the majority of the scientific and technical expertise in the field that can be found in ESO member states. △ Less

Submitted 19 July, 2024; originally announced July 2024.

Comments: SPIE astronomical telescope and instrumentation 2024, in press

Extended GeV $γ$-ray emission around the star forming region of the W3 complex

Authors: Qihang Wu, Xiaona Sun, Ruizhi Yang, Tingting Ge, Yunfeng Liang, Enwei Liang

Abstract: We analyze the GeV $γ$-ray emission from the W3 complex using about 14 years of Pass 8 data recorded by the $\it Fermi$ Large Area Telescope (\textit{Fermi}-LAT). We resolve the $γ$-ray emissions around W3 into two components: an elliptical Gaussian overlapping with the molecular gas and a point-like source near the cluster W3 Main. The pion-bump feature of SED for the elliptical Gaussian together… ▽ More We analyze the GeV $γ$-ray emission from the W3 complex using about 14 years of Pass 8 data recorded by the $\it Fermi$ Large Area Telescope (\textit{Fermi}-LAT). We resolve the $γ$-ray emissions around W3 into two components: an elliptical Gaussian overlapping with the molecular gas and a point-like source near the cluster W3 Main. The pion-bump feature of SED for the elliptical Gaussian together with the better fitting result of pion decay model favor the hadronic origin. We further argue that the cosmic rays (CRs) could originate from the interactions between cluster winds and the shock produced by the SNR HB3. The point-like source positionally coincident with the star cluster W3 Main indicates it may be directly powered by near clusters, while its fainter $γ$-ray emissions below 10 GeV is possibly due to the shelter from dense gas making the low-energy CRs incapable of penetrating the dense materials. Meanwhile, we cannot rule out that the $γ$-ray emissions originate from the interaction of accelerated protons in SNR with the ambient gas. △ Less

Submitted 27 June, 2024; originally announced June 2024.

Observational characteristics of circum-planetary-mass-object disks in the era of James Webb Space Telescope

Authors: Xilei Sun, Pinghui Huang, Ruobing Dong, Shang-Fei Liu

Abstract: Recent observations have confirmed circumplanetary disks (CPDs) embedded in parental protoplanetary disks (PPDs). On the other hand, planetary-mass companions (PMCs) and planetary-mass objects (PMOs) are likely to harbor their own accretion disks. Unlike PPDs, CPDs and other disks around planet analogues are generally too compact to be spatially resolved by current instrumentation. In this study,… ▽ More Recent observations have confirmed circumplanetary disks (CPDs) embedded in parental protoplanetary disks (PPDs). On the other hand, planetary-mass companions (PMCs) and planetary-mass objects (PMOs) are likely to harbor their own accretion disks. Unlike PPDs, CPDs and other disks around planet analogues are generally too compact to be spatially resolved by current instrumentation. In this study, we generate over 4,000 spectral energy distributions (SEDs) of circum-PMO-disks (CPMODs) with various host temperature and disk properties, which can be categorized into four prototypes, i.e., full, pre-transitional, transitional and evolved CPMODs. We propose a classification scheme based on their near-to-mid-infrared colors. Using those CPMOD models, we synthesize JWST (NIRCam and MIRI) photometry for F444W, F1000W and F2550W wide filters. We show F444W - F1000W and F444 - F2550W colors can be applied to distinguish different types of CPMODs, especially for those around hot hosts. Our results indicate that the ongoing and future JWST observations are promising to unveil structures and properties of CPMODs. △ Less

Submitted 13 June, 2024; originally announced June 2024.

Comments: 18 pages, 7 figures, accepted for publication in ApJ

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

Faraday moments of the Southern Twenty-centimeter All-sky Polarization Survey (STAPS)

Authors: N. Raycheva, M. Haverkorn, S. Ideguchi, J. M. Stil, X. Sun, J. L. Han, E. Carretti, X. Y. Gao, A. Bracco, S. E. Clark, J. M. Dickey, B. M. Gaensler, A. Hill, T. Landecker, A. Ordog, A. Seta, M. Tahani, M. Wolleben

Abstract: Faraday tomography of broadband radio polarization surveys enables us to study magnetic fields and their interaction with the interstellar medium (ISM). Such surveys include the Global Magneto-Ionic Medium Survey (GMIMS), which covers the northern and southern hemispheres at $\sim$ 300-1800 MHz. In this work, we used the GMIMS High Band South (1328-1768 MHz), also named the Southern Twenty-centi… ▽ More Faraday tomography of broadband radio polarization surveys enables us to study magnetic fields and their interaction with the interstellar medium (ISM). Such surveys include the Global Magneto-Ionic Medium Survey (GMIMS), which covers the northern and southern hemispheres at $\sim$ 300-1800 MHz. In this work, we used the GMIMS High Band South (1328-1768 MHz), also named the Southern Twenty-centimeter All-sky Polarization Survey (STAPS), which observes the southern sky at a resolution of 18$\arcmin$. To extract the key parameters of the magnetized ISM from STAPS, we computed the Faraday moments of the tomographic data cubes. These moments include the total polarized intensity, the mean Faraday depth weighted by the polarized intensity, the weighted dispersion of the Faraday spectrum, and its skewness. We compared the Faraday moments to those calculated over the same frequency range in the northern sky (using the Dominion Radio Astrophysical Observatory, DRAO), in a strip of $360\degr \times 30\degr$ that overlaps with STAPS coverage. We find that the total polarized intensity is generally dominated by diffuse emission that decreases at longitudes of $l \leq 300\degr$. The Faraday moments reveal a variety of polarization structures. Low-intensity regions at high latitudes usually have a single Faraday depth component. Due to its insufficiently large frequency coverage, STAPS cannot detect Faraday thick structures. Comparing the Faraday depths from STAPS to total rotation measures from extragalactic sources suggests that STAPS frequencies are high enough that the intervening ISM causes depolarization to background emission at intermediate and high Galactic latitudes. Where they overlap, the STAPS and DRAO surveys exhibit broad correspondence but differ in polarized intensity by a factor of $\sim$1.8. △ Less

Submitted 10 June, 2024; originally announced June 2024.

Comments: Revision in progress in Astronomy and Astrophysics

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

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

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

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

Comments: 15 pages, 9 figures

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

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

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

Submitted 13 May, 2024; originally announced May 2024.

Comments: 11 pages, 5 figures

Soft X-ray prompt emission from a high-redshift gamma-ray burst EP240315a

Authors: Y. Liu, H. Sun, D. Xu, D. S. Svinkin, J. Delaunay, N. R. Tanvir, H. Gao, C. Zhang, Y. Chen, X. -F. Wu, B. Zhang, W. Yuan, J. An, G. Bruni, D. D. Frederiks, G. Ghirlanda, J. -W. Hu, A. Li, C. -K. Li, J. -D. Li, D. B. Malesani, L. Piro, G. Raman, R. Ricci, E. Troja , et al. (170 additional authors not shown)

Abstract: Long gamma-ray bursts (GRBs) are believed to originate from core collapse of massive stars. High-redshift GRBs can probe the star formation and reionization history of the early universe, but their detection remains rare. Here we report the detection of a GRB triggered in the 0.5--4 keV band by the Wide-field X-ray Telescope (WXT) on board the Einstein Probe (EP) mission, designated as EP240315a,… ▽ More Long gamma-ray bursts (GRBs) are believed to originate from core collapse of massive stars. High-redshift GRBs can probe the star formation and reionization history of the early universe, but their detection remains rare. Here we report the detection of a GRB triggered in the 0.5--4 keV band by the Wide-field X-ray Telescope (WXT) on board the Einstein Probe (EP) mission, designated as EP240315a, whose bright peak was also detected by the Swift Burst Alert Telescope and Konus-Wind through off-line analyses. At a redshift of $z=4.859$, EP240315a showed a much longer and more complicated light curve in the soft X-ray band than in gamma-rays. Benefiting from a large field-of-view ($\sim$3600 deg$^2$) and a high sensitivity, EP-WXT captured the earlier engine activation and extended late engine activity through a continuous detection. With a peak X-ray flux at the faint end of previously known high-$z$ GRBs, the detection of EP240315a demonstrates the great potential for EP to study the early universe via GRBs. △ Less

Submitted 25 April, 2024; originally announced April 2024.

Comments: 41 pages, 8 figures, 7 tables

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.

GeV gamma-ray emission in the field of young massive star cluster RCW 38

Authors: Ting-Ting Ge, Xiao-Na Sun, Rui-Zhi Yang, Pak-Hin Thomas Tam, Ming-Xuan Lu, En-Wei Liang

Abstract: We report the detection of gamma-ray emission by the Fermi Large Area Telescope (Fermi-LAT) towards the young massive star cluster RCW 38 in the 1-500 GeV photon energy range. We found spatially extended GeV emission towards the direction of RCW 38, which is best modelled by a Gaussian disc of 0.23$°$ radius with a significance of the extension is $\sim 11.4 σ$. Furthermore, the spatial correlatio… ▽ More We report the detection of gamma-ray emission by the Fermi Large Area Telescope (Fermi-LAT) towards the young massive star cluster RCW 38 in the 1-500 GeV photon energy range. We found spatially extended GeV emission towards the direction of RCW 38, which is best modelled by a Gaussian disc of 0.23$°$ radius with a significance of the extension is $\sim 11.4 σ$. Furthermore, the spatial correlation with the ionized and molecular gas content favors the hadronic origin of the gamma-ray emission. The gamma-ray spectrum of RCW 38 has a relatively hard photon index of $2.44 \pm 0.03$, which is similar to other young massive star clusters. We argue that the diffuse GeV gamma-ray emission in this region likely originates from the interaction of accelerated protons in the stellar cluster with the ambient gas. △ Less

Submitted 28 March, 2024; originally announced March 2024.

Comments: 8 pages, 5 figures, accepted for publication in MNRAS. arXiv admin note: text overlap with arXiv:2210.01352

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)

Simulation Studies for the First Pathfinder of the CATCH Space Mission

Authors: Yiming Huang, Juan Zhang, Lian Tao, Zhengwei Li, Donghua Zhao, Qian-Qing Yin, Xiangyang Wen, Jingyu Xiao, Chen Zhang, Shuang-Nan Zhang, Shaolin Xiong, Qingcui Bu, Jirong Cang, Dezhi Cao, Wen Chen, Siran Ding, Min Gao, Yang Gao, Shujin Hou, Liping Jia, Ge Jin, Dalin Li, Jinsong Li, Panping Li, Yajun Li , et al. (20 additional authors not shown)

Abstract: The Chasing All Transients Constellation Hunters (CATCH) space mission is an intelligent constellation consisting of 126 micro-satellites in three types (A, B, and C), designed for X-ray observation with the objective of studying the dynamic universe. Currently, we are actively developing the first Pathfinder (CATCH-1) for the CATCH mission, specifically for type-A satellites. CATCH-1 is equipped… ▽ More The Chasing All Transients Constellation Hunters (CATCH) space mission is an intelligent constellation consisting of 126 micro-satellites in three types (A, B, and C), designed for X-ray observation with the objective of studying the dynamic universe. Currently, we are actively developing the first Pathfinder (CATCH-1) for the CATCH mission, specifically for type-A satellites. CATCH-1 is equipped with Micro Pore Optics (MPO) and a 4-pixel Silicon Drift Detector (SDD) array. To assess its scientific performance, including the effective area of the optical system, on-orbit background, and telescope sensitivity, we employ the Monte Carlo software Geant4 for simulation in this study. The MPO optics exhibit an effective area of $41$ cm$^2$ at the focal spot for 1 keV X-rays, while the entire telescope system achieves an effective area of $29$ cm$^2$ at 1 keV when taking into account the SDD detector's detection efficiency. The primary contribution to the background is found to be from the Cosmic X-ray Background. Assuming a 625 km orbit with an inclination of $29^\circ$, the total background for CATCH-1 is estimated to be $8.13\times10^{-2}$ counts s$^{-1}$ in the energy range of 0.5--4 keV. Based on the background within the central detector and assuming a Crab-like source spectrum, the estimated ideal sensitivity could achieve $1.9\times10^{-12}$ erg cm$^{-2}$ s$^{-1}$ for an exposure of 10$^4$ s in the energy band of 0.5--4 keV. Furthermore, after simulating the background caused by low-energy charged particles near the geomagnetic equator, we have determined that there is no need to install a magnetic deflector. △ Less

Submitted 23 February, 2024; originally announced February 2024.

Nuclear mass table in deformed relativistic Hartree-Bogoliubov theory in continuum, II: Even-$Z$ nuclei

Authors: DRHBc Mass Table Collaboration, Peng Guo, Xiaojie Cao, Kangmin Chen, Zhihui Chen, Myung-Ki Cheoun, Yong-Beom Choi, Pak Chung Lam, Wenmin Deng, Jianmin Dong, Pengxiang Du, Xiaokai Du, Kangda Duan, Xiaohua Fan, Wei Gao, Lisheng Geng, Eunja Ha, Xiao-Tao He, Jinniu Hu, Jingke Huang, Kun Huang, Yanan Huang, Zidan Huang, Kim Da Hyung, Hoi Yat Chan , et al. (58 additional authors not shown)

Abstract: The mass table in the deformed relativistic Hartree-Bogoliubov theory in continuum (DRHBc) with the PC-PK1 density functional has been established for even-$Z$ nuclei with $8\le Z\le120$, extended from the previous work for even-even nuclei [Zhang $\it{et.~al.}$ (DRHBc Mass Table Collaboration), At. Data Nucl. Data Tables 144, 101488 (2022)]. The calculated binding energies, two-nucleon and one-ne… ▽ More The mass table in the deformed relativistic Hartree-Bogoliubov theory in continuum (DRHBc) with the PC-PK1 density functional has been established for even-$Z$ nuclei with $8\le Z\le120$, extended from the previous work for even-even nuclei [Zhang $\it{et.~al.}$ (DRHBc Mass Table Collaboration), At. Data Nucl. Data Tables 144, 101488 (2022)]. The calculated binding energies, two-nucleon and one-neutron separation energies, root-mean-square (rms) radii of neutron, proton, matter, and charge distributions, quadrupole deformations, and neutron and proton Fermi surfaces are tabulated and compared with available experimental data. A total of 4829 even-$Z$ nuclei are predicted to be bound, with an rms deviation of 1.477 MeV from the 1244 mass data. Good agreement with the available experimental odd-even mass differences, $α$ decay energies, and charge radii is also achieved. The description accuracy for nuclear masses and nucleon separation energies as well as the prediction for drip lines is compared with the results obtained from other relativistic and nonrelativistic density functional. The comparison shows that the DRHBc theory with PC-PK1 provides an excellent microscopic description for the masses of even-$Z$ nuclei. The systematics of the nucleon separation energies, odd-even mass differences, pairing energies, two-nucleon gaps, $α$ decay energies, rms radii, quadrupole deformations, potential energy curves, neutron density distributions, and neutron mean-field potentials are discussed. △ Less

Submitted 10 June, 2024; v1 submitted 5 February, 2024; originally announced February 2024.

Comments: 394 pages, 17 figures, 2 tables, published in Atomic Data and Nuclear Data Tables, data file in the TXT form is available for download under "Ancillary files"

Journal ref: Peng Guo, et. al. (DRHBc Mass Table Collaboration), Atomic Data and Nuclear Data Tables 158 (2024) 101661

Spatial resolution effects on the solar open flux estimates

Authors: Ivan Milic, Rebecca Centeno, Xudong Sun, Matthias Rempel, Jaime de la Cruz Rodriguez

Abstract: Spectropolarimetric observations used to infer the solar magnetic fields are obtained with a limited spatial resolution. The effects of this limited resolution on the inference of the open flux over the observed region have not been extensively studied. We aim to characterize the biases that arise in the inference of the mean flux density by performing an end-to-end study that involves the generat… ▽ More Spectropolarimetric observations used to infer the solar magnetic fields are obtained with a limited spatial resolution. The effects of this limited resolution on the inference of the open flux over the observed region have not been extensively studied. We aim to characterize the biases that arise in the inference of the mean flux density by performing an end-to-end study that involves the generation of synthetic data, its interpretation (inversion), and a comparison of the results with the original model. We synthesized polarized spectra of the two magnetically sensitive lines of neutral iron around 630\,nm from a state-of-the-art numerical simulation of the solar photosphere. We then performed data degradation to simulate the effect of the telescope with a limited angular resolution and interpreted (inverted) the data using a Milne-Eddington spectropolarimetric inversion code. We then studied the dependence of the inferred parameters on the telescope resolution. The results show a significant decrease in the mean magnetic flux density -- related to the open flux observed at the disk center -- with decreasing telescope resolution. The original net magnetic field flux is fully resolved by a 1m telescope, but a 20\,cm aperture telescope yields a 30\% smaller value. Even in the fully resolved case, the result is still biased due to the corrugation of the photospheric surface. Even the spatially averaged quantities, such as the open magnetic flux in the observed region, are underestimated when the magnetic structures are unresolved. The reason for this is the presence of nonlinearities in the magnetic field inference process. This effect might have implications for the modeling of large-scale solar magnetic fields; for example, those corresponding to the coronal holes, or the polar magnetic fields, which are relevant to our understanding of the solar cycle. △ Less

Submitted 4 February, 2024; originally announced February 2024.

Comments: To be published in A&A

A Machine Learning made Catalog of FR-II Radio Galaxies from the FIRST Survey

Authors: Bao-Qiang Lao, Xiao-Long Yang, Sumit Jaiswal, Prashanth Mohan, Xiao-Hui Sun, Sheng-Li Qin, Ru-Shuang Zhao

Abstract: We present an independent catalog (FRIIRGcat) of 45,241 Fanaroff-Riley Type II (FR-II) radio galaxies compiled from the Very Large Array Faint Images of the Radio Sky at Twenty-centimeters (FIRST) survey and employed the deep learning method. Among them, optical and/or infrared counterparts are identified for 41,425 FR-IIs. This catalog spans luminosities… ▽ More We present an independent catalog (FRIIRGcat) of 45,241 Fanaroff-Riley Type II (FR-II) radio galaxies compiled from the Very Large Array Faint Images of the Radio Sky at Twenty-centimeters (FIRST) survey and employed the deep learning method. Among them, optical and/or infrared counterparts are identified for 41,425 FR-IIs. This catalog spans luminosities $2.63\times10^{22}\leq L_{\rm rad}\leq6.76\times10^{29}\,{\rm W}\,{\rm Hz}^{-1}$ and redshifts up to $z=5.01$. The spectroscopic classification indicates that there are 1431 low-excitation radio galaxies and 260 high-excitation radio galaxies. Among the spectroscopically identified sources, black hole masses are estimated for 4837 FR-IIs, which are in $10^{7.5}\lesssim M_{\rm BH}\lesssim 10^{9.5}$ $M_{\odot}$. Interestingly, this catalog reveals a couple of giant radio galaxies (GRGs), which are already in the existing GRG catalog, confirming the efficiency of this FR-II catalog. Furthermore, 284 new GRGs are unveiled in this new FR-II sample; they have the largest projected sizes ranging from 701 to 1209 kpc and are located at redshifts $0.31<z<2.42$. Finally, we explore the distribution of the jet position angle and it shows that the faint Images of the FIRST images are significantly affected by the systematic effect (the observing beams). The method presented in this work is expected to be applicable to the radio sky surveys that are currently being conducted because they have finely refined telescope arrays. On the other hand, we are expecting that further new methods will be dedicated to solving this problem. △ Less

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

Comments: 24 pages, 17 figures, Accepted by Research in Astronomy and Astrophysics

Journal ref: 2024, Research in Astronomy and Astrophysics Volume 24, Pages 035021

Detector performance of the Gamma-ray Transient Monitor onboard DRO-A Satellite

Authors: Pei-Yi Feng, Zheng-Hua An, Da-Li Zhang, Chen-Wei Wang, Chao Zheng, Sheng Yang, Shao-Lin Xiong, Jia-Cong Liu, Xin-Qiao Li, Ke Gong, Xiao-Jing Liu, Min Gao, Xiang-Yang Wen, Ya-Qing liu, Xiao-Yun Zhao, Fan Zhang, Xi-Lei Sun, Hong Lu

Abstract: Gamma-ray Transient Monitor (GTM) is an all-sky monitor onboard the Distant Retrograde Orbit-A (DRO-A) satellite with the scientific objective of detecting gamma-ray transients ranging from 20 keV to 1 MeV. GTM is equipped with 5 Gamma-ray Transient Probe (GTP) detector modules, utilizing the NaI(Tl) scintillator coupled with a SiPM array. To reduce the SiPM noise, GTP makes use of a dedicated dua… ▽ More Gamma-ray Transient Monitor (GTM) is an all-sky monitor onboard the Distant Retrograde Orbit-A (DRO-A) satellite with the scientific objective of detecting gamma-ray transients ranging from 20 keV to 1 MeV. GTM is equipped with 5 Gamma-ray Transient Probe (GTP) detector modules, utilizing the NaI(Tl) scintillator coupled with a SiPM array. To reduce the SiPM noise, GTP makes use of a dedicated dual-channel coincident readout design. In this work, we firstly studied the impact of different coincidence times on detection efficiency and ultimately selected the 500 ns time coincidence window for offline data processing. To test the performance of GTPs and validate the Monte Carlo simulated energy response, we conducted comprehensive ground calibration tests using Hard X-ray Calibration Facility (HXCF) and radioactive sources, including energy response, detection efficiency, spatial response, bias-voltage response, and temperature dependence. We extensively presented the ground calibration results, and validated the design and mass model of GTP detector. These work paved the road for the in-flight observation and science data analysis. △ Less

Submitted 15 January, 2024; originally announced January 2024.

Comments: 13 pages, 25 figures

The Decay Process of an α-configuration Sunspot

Authors: Yang Peng, Zhi-Ke Xue, Xiao-Li Yan, Aimee A. Norton, Zhong-Quan Qu, Jin-Cheng Wang, Zhe Xu, Li-Heng Yang, Qiao-Ling Li, Li-Ping Yang, Xia Sun

Abstract: The decay of sunspot plays a key role in magnetic flux transportation in solar active regions (ARs). To better understand the physical mechanism of the entire decay process of a sunspot, an α-configuration sunspot in AR NOAA 12411 was studied. Based on the continuum intensity images and vector magnetic field data with stray light correction from Solar Dynamics Observatory/Helioseismic and Magnetic… ▽ More The decay of sunspot plays a key role in magnetic flux transportation in solar active regions (ARs). To better understand the physical mechanism of the entire decay process of a sunspot, an α-configuration sunspot in AR NOAA 12411 was studied. Based on the continuum intensity images and vector magnetic field data with stray light correction from Solar Dynamics Observatory/Helioseismic and Magnetic Imager, the area, vector magnetic field and magnetic flux in the umbra and penumbra are calculated with time, respectively. Our main results are as follows: (1) The decay curves of the sunspot area in its umbra, penumbra, and whole sunspot take the appearance of Gaussian profiles. The area decay rates of the umbra, penumbra and whole sunspot are -1.56 MSH/day, -12.61 MSH/day and -14.04 MSH/day, respectively; (2) With the decay of the sunspot, the total magnetic field strength and the vertical component of the penumbra increase, and the magnetic field of the penumbra becomes more vertical. Meanwhile, the total magnetic field strength and vertical magnetic field strength for the umbra decrease, and the inclination angle changes slightly with an average value of about 20°; (3) The magnetic flux decay curves of the sunspot in its umbra, penumbra, and whole sunspot exhibit quadratic patterns, their magnetic flux decay rates of the umbra, penumbra and whole sunspot are -9.84 * 10^19 Mx/day, -1.59 * 10^20 Mx/day and -2.60 * 10^20 Mx/day , respectively. The observation suggests that the penumbra may be transformed into the umbra, resulting in the increase of the average vertical magnetic field strength and the reduction of the inclination angle in the penumbra during the decay of the sunspot. △ Less

Submitted 11 January, 2024; originally announced January 2024.

The Intrinsic Energy Resolution of LaBr$_3$(Ce) Crystal for GECAM

Authors: Pei-Yi Feng, Xi-Lei Sun, Cheng-Er Wang, Yong Deng, Zheng-Hua An, Da-Li Zhang, Chao Zheng, Xin-Qiao Li, Shao-Lin Xiong, Hong Lu

Abstract: The intrinsic resolution is the primary limitation on the total energy resolution of LaBr$_3$(Ce) crystal. This intrinsic resolution arises from two effects: fluctuations occurring in the process of energy transfer to luminescent centers within the LaBr$_3$(Ce) crystal and the LaBr$_3$(Ce) crystal's non-proportional luminescence. Presently, experimental measurements regarding the intrinsic resolut… ▽ More The intrinsic resolution is the primary limitation on the total energy resolution of LaBr$_3$(Ce) crystal. This intrinsic resolution arises from two effects: fluctuations occurring in the process of energy transfer to luminescent centers within the LaBr$_3$(Ce) crystal and the LaBr$_3$(Ce) crystal's non-proportional luminescence. Presently, experimental measurements regarding the intrinsic resolution of LaBr$_3$(Ce) crystal are scarce, and the underlying physical mechanisms remain incompletely understood. In this paper, we aim to elucidate the concept of intrinsic resolution. We investigated the entire physical process of luminescence following energy deposition in the LaBr$_3$(Ce) crystal, quantifying the various components in the total energy resolution. We conducted a series of experimental measurements and Geant4 simulations, determining the intrinsic resolution of LaBr$_3$(Ce) crystal to 100 keV electrons as 2.12%. The non-proportionality contributes significantly at 1.43%, while fluctuations in the energy transfer process accounted for 0.27%. It is evident that non-proportionality in light output constitutes the primary source of intrinsic resolution. Horizontal and vertical unevenness in light collection contributed 0.25% and 0.07%, respectively. Statistical fluctuations showed the largest impact on the total energy resolution, at 2.86%. The contribution from fluctuations in single-photoelectron events was 0.77%. Furthermore, we reconstructed the photon response using Geant4, and the consistency between the simulated relative light yield and the experimentally measured one confirmed the reliability of the LaBr$_3$(Ce) detector mass model employed in the simulation. △ Less

Submitted 30 December, 2023; originally announced January 2024.

Comments: 11 pages, 16 figures

Onset mechanism of an inverted U-shaped solar filament eruption revealed by NVST, SDO, and STEREO-A observations

Authors: Jincheng Wang, Xiaoli Yan, Qiangwei Cai, Zhike Xue, Liheng Yang, Qiaoling Li, Zhe Xu, Yunfang Cai, Liping Yang, Yang Peng, Xia Sun, Xinsheng Zhang, Yian Zhou

Abstract: Utilizing observations from the New Vacuum Solar Telescope (NVST), Solar Dynamics Observatory (SDO), and Solar Terrestrial Relations Observatory-Ahead (STEREO-A), we investigate the event from two distinct observational perspectives: on the solar disk using NVST and SDO, and on the solar limb using STEREO-A. We employ both a non-linear force-free field model and a potential field model to reconstr… ▽ More Utilizing observations from the New Vacuum Solar Telescope (NVST), Solar Dynamics Observatory (SDO), and Solar Terrestrial Relations Observatory-Ahead (STEREO-A), we investigate the event from two distinct observational perspectives: on the solar disk using NVST and SDO, and on the solar limb using STEREO-A. We employ both a non-linear force-free field model and a potential field model to reconstruct the coronal magnetic field, aiming to understand its magnetic properties. Two precursor jet-like activities were observed before the eruption, displaying an untwisted rotation. The second activity released an estimated twist of over two turns. During these two jet-like activities, Y-shaped brightenings, newly emerging magnetic flux accompanied by magnetic cancellation, and the formation of newly moving fibrils were identified. Combining these observational features, it can be inferred that these two precursor jet-like activities released the magnetic field constraining the filament and were triggered by newly emerging magnetic flux. Before the filament eruption, it was observed that some moving flows had been ejected from the site as the onset of two jet-like activities, indicating the same physical process as two jet-like activities. Extrapolations revealed that the filament laid under the height of the decay index of 1.0 and had strong magnetic field (540 Gauss) and a high twisted number (2.4 turns) before the eruption. An apparent rotational motion was observed during the filament eruption. We deduce that the solar filament, exhibiting an inverted U-shape, is a significantly twisted flux rope. The eruption of the filament was initiated by the release of constraining magnetic fields through continuous magnetic reconnection. This reconnection process was triggered by the emergence of newly magnetic flux. △ Less

Submitted 30 December, 2023; originally announced January 2024.

Comments: 18 pages, 12 figures, accepted for publication in Astronomy & Astrophysics

The Energy Response of LaBr3(Ce), LaBr3(Ce,Sr) and NaI(Tl) Crystals for GECAM

Authors: Pei-Yi Feng, Xi-Lei Sun, Zheng-Hua An, Yong Deng, Cheng-Er Wang, Huang Jiang, Jun-Jie Li, Da-Li Zhang, Xin-Qiao Li, Shao-Lin Xiong, Chao Zheng, Ke Gong, Sheng Yang, Xiao-Jing Liu, Min Gao, Xiang-Yang Wen, Ya-Qing Liu, Yan-Bing Xu, Xiao-Yun Zhao, Jia-Cong Liu, Fan Zhang, Hong Lu

Abstract: The GECAM series of satellites utilize LaBr3(Ce), LaBr3(Ce,Sr), and NaI(Tl) crystals as sensitive materials for gamma-ray detectors (GRDs). To investigate the non-linearity in the detection of low-energy gamma rays and address errors in the E-C relationship calibration, comprehensive tests and comparative studies of the non-linearity of these three crystals were conducted using Compton electrons,… ▽ More The GECAM series of satellites utilize LaBr3(Ce), LaBr3(Ce,Sr), and NaI(Tl) crystals as sensitive materials for gamma-ray detectors (GRDs). To investigate the non-linearity in the detection of low-energy gamma rays and address errors in the E-C relationship calibration, comprehensive tests and comparative studies of the non-linearity of these three crystals were conducted using Compton electrons, radioactive sources, and mono-energetic X-rays. The non-linearity test results for Compton electrons and X-rays displayed substantial differences, with all three crystals showing higher non-linearity for X-rays and gamma-rays than for Compton electrons. Despite LaBr3(Ce) and LaBr3(Ce,Sr) crystals having higher absolute light yields, they exhibited a noticeable non-linear decrease in light yield, especially at energies below 400 keV. The NaI(Tl) crystal demonstrated excess light output in the 6~200 keV range, reaching a maximum excess of 9.2% at 30 keV in X-ray testing and up to 15.5% at 14 keV during Compton electron testing, indicating a significant advantage in the detection of low-energy gamma rays. Furthermore, this paper explores the underlying causes of the observed non-linearity in these crystals. This study not only elucidates the detector responses of GECAM, but also marks the inaugural comprehensive investigation into the non-linearity of domestically produced lanthanum bromide and sodium iodide crystals. △ Less

Submitted 27 December, 2023; originally announced December 2023.

Comments: 12pages, 16 figures

Stellar Flares Are Far-Ultraviolet Luminous

Authors: Vera L. Berger, Jason T. Hinkle, Michael A. Tucker, Benjamin J. Shappee, Jennifer L. van Saders, Daniel Huber, Jeffrey W. Reep, Xudong Sun, Kai E. Yang

Abstract: We identify 182 flares on 158 stars within 100 pc of the Sun in both the near-ultraviolet (NUV: 1750-2750 Å) and far-ultraviolet (FUV: 1350-1750 Å) using high-cadence light curves from the Galaxy Evolution Explorer (GALEX). Ultraviolet (UV) emission from stellar flares plays a crucial role in determining the habitability of exoplanetary systems. However, whether such UV emission promotes or threat… ▽ More We identify 182 flares on 158 stars within 100 pc of the Sun in both the near-ultraviolet (NUV: 1750-2750 Å) and far-ultraviolet (FUV: 1350-1750 Å) using high-cadence light curves from the Galaxy Evolution Explorer (GALEX). Ultraviolet (UV) emission from stellar flares plays a crucial role in determining the habitability of exoplanetary systems. However, whether such UV emission promotes or threatens such life depends strongly on the energetics of these flares. Most studies assessing the effect of flares on planetary habitability assume a 9000 K blackbody spectral energy distribution that produces more NUV flux than FUV flux ($R \equiv F_{\rm FUV} / F_{\rm NUV} \approx \frac{1}{6}$). Instead, we observe the opposite with the excess FUV reaching $R \approx \frac{1}{2} - 2$, roughly $3-12$ times the expectation of a 9000 K blackbody. The ratio of FUV to NUV time-integrated flare energies is 3.0 times higher on average than would be predicted by a constant 9000 K blackbody during the flare. Finally, we find that the FUV/NUV ratio at peak tentatively correlates ($\sim 2 σ$ significance) both with total UV flare energy and with the G - RP color of the host star. On average, we observe higher FUV/NUV ratios at peak in $E_{\text{UV}}>10^{32}$ erg flares and in flares on fully convective stars. △ Less

Submitted 19 December, 2023; originally announced December 2023.

Comments: Submitted to MNRAS, comments welcome

Ground Calibration Result of the Lobster Eye Imager for Astronomy

Authors: Huaqing Cheng, Zhixing Ling, Chen Zhang, Xiaojin Sun, Shengli Sun, Yuan Liu, Yanfeng Dai, Zhenqing Jia, Haiwu Pan, Wenxin Wang, Donghua Zhao, Yifan Chen, Zhiwei Cheng, Wei Fu, Yixiao Han, Junfei Li, Zhengda Li, Xiaohao Ma, Yulong Xue, Ailiang Yan, Qiang Zhang, Yusa Wang, Xiongtao Yang, Zijian Zhao, Weimin Yuan

Abstract: We report on results of the on-ground X-ray calibration of the Lobster Eye Imager for Astronomy (LEIA), an experimental space wide-field (18.6*18.6 square degrees) X-ray telescope built from novel lobster eye mirco-pore optics. LEIA was successfully launched on July 27, 2022 onboard the SATech-01 satellite. To achieve full characterisation of its performance before launch, a series of tests and ca… ▽ More We report on results of the on-ground X-ray calibration of the Lobster Eye Imager for Astronomy (LEIA), an experimental space wide-field (18.6*18.6 square degrees) X-ray telescope built from novel lobster eye mirco-pore optics. LEIA was successfully launched on July 27, 2022 onboard the SATech-01 satellite. To achieve full characterisation of its performance before launch, a series of tests and calibrations have been carried out at different levels of devices, assemblies and the complete module. In this paper, we present the results of the end-to-end calibration campaign of the complete module carried out at the 100-m X-ray Test Facility at IHEP. The PSF, effective area and energy response of the detectors were measured in a wide range of incident directions at several X-ray line energies. The distributions of the PSF and effective areas are roughly uniform across the FoV, in large agreement with the prediction of lobster-eye optics. The mild variations and deviations from the prediction of idealized lobster-eye optics can be understood to be caused by the imperfect shapes and alignment of the micro-pores as well as the obscuration by the supporting frames, which can be well reproduced by MC simulations. The spatial resolution of LEIA defined by the FWHM of the focal spot ranges from 4-8 arcmin with a median of 5.7. The measured effective areas are in range of 2-3 $cm^2$ at ~1.25 keV across the entire FoV, and its dependence on photon energy is in large agreement with simulations. The gains of the CMOS sensors are in range of 6.5-6.9 eV/DN, and the energy resolutions in the range of ~120-140 eV at 1.25 keV and ~170-190 eV at 4.5 keV. These results have been ingested into the calibration database and applied to the analysis of the scientific data acquired by LEIA. This work paves the way for the calibration of the Wide-field X-Ray Telescope modules of the Einstein Probe mission. △ Less

Submitted 11 December, 2023; originally announced December 2023.

Comments: 24 pages, 13 figures. Submitted to Experimental Astronomy

New Continuum Observations of the Andromeda galaxy M31 with FAST

Authors: Wenjun Zhang, Xiaohui Sun, Jie Wang

Abstract: We present a new total intensity image of M31 at 1.248 GHz, observed with the Five-hundred-meter Aperture Spherical radio telescope (FAST) with an angular resolution of 4 arcmin and a sensitivity of about 16 mK. The new FAST image clearly reveals weak emission outside the ring due to its high sensitivity on large-scale structures. We derive a scale length of 2.7 kpc for the cosmic ray electrons an… ▽ More We present a new total intensity image of M31 at 1.248 GHz, observed with the Five-hundred-meter Aperture Spherical radio telescope (FAST) with an angular resolution of 4 arcmin and a sensitivity of about 16 mK. The new FAST image clearly reveals weak emission outside the ring due to its high sensitivity on large-scale structures. We derive a scale length of 2.7 kpc for the cosmic ray electrons and find that the cosmic ray electrons propagate mainly through diffusion by comparing the scale length at 4.8 GHz. The spectral index of the total intensity varies along the ring, which can be attributed to the variation of the spectra of synchrotron emission. This variation is likely caused by the change of star formation rates along the ring. We find that the azimuthal profile of the non-thermal emission can be interpreted by an axisymmetric large-scale magnetic field with varying pitch angle along the ring, indicating a complicated magnetic field configuration in M31. △ Less

Submitted 1 December, 2023; originally announced December 2023.

Comments: 9 papges

G213.0$-$0.6, a true supernova remnant or just an HII region?

Authors: X. Y. Gao, C. J. Wu, X. H. Sun, W. Reich, J. L. Han

Abstract: G213.0$-$0.6 is a faint extended source situated in the anti-center region of the Galactic plane. It has been classified as a shell-type supernova remnant (SNR) based on its shell-like morphology, steep radio continuum spectrum, and high ratio of [S II]/H$α$. With new optical emission line data of H$α$, [S II], and [N II] recently observed by the Large Sky Area Multi-Object Fiber Spectroscopic Tel… ▽ More G213.0$-$0.6 is a faint extended source situated in the anti-center region of the Galactic plane. It has been classified as a shell-type supernova remnant (SNR) based on its shell-like morphology, steep radio continuum spectrum, and high ratio of [S II]/H$α$. With new optical emission line data of H$α$, [S II], and [N II] recently observed by the Large Sky Area Multi-Object Fiber Spectroscopic Telescope, the ratios of [S II]/H$α$ and [N II]/H$α$ are re-assessed. The lower values than those previously reported put G213.0$-$0.6 around the borderline of SNR-HII region classification. We decompose the steep-spectrum synchrotron and the flat-spectrum thermal free-free emission in the area of G213.0$-$0.6 with multi-frequency radio continuum data. G213.0$-$0.6 is found to show a flat spectrum, in conflict with the properties of a shell-type SNR. Such a result is further confirmed by TT-plots made between the 863-MHz, 1.4-GHz, and 4.8-GHz data. Combining the evidence extracted in both optical and radio continuum, we argue that G213.0$-$0.6 is possibly not an SNR, but an HII region instead. The $V_{LSR}$ pertaining to the H$α$ filaments places G213.0$-$0.6 approximately 1.9 kpc away in the Perseus Arm. △ Less

Submitted 19 November, 2023; originally announced November 2023.

Comments: 9 pages, 4 figures, accepted for publication in MNRAS

A Possible Mechanism for "Late Phase" in Stellar White-Light Flares

Authors: Kai E. Yang, Xudong Sun, Graham S. Kerr, Hugh S. Hudson

Abstract: M-dwarf flares observed by the \textit{Transiting Exoplanet Survey Satellite} (\textit{TESS}) sometimes exhibit a "peak-bump" light-curve morphology, characterized by a secondary, gradual peak well after the main, impulsive peak. A similar "late phase" is frequently detected in solar flares observed in the extreme-ultraviolet from longer hot coronal loops distinct from the impulsive flare structur… ▽ More M-dwarf flares observed by the \textit{Transiting Exoplanet Survey Satellite} (\textit{TESS}) sometimes exhibit a "peak-bump" light-curve morphology, characterized by a secondary, gradual peak well after the main, impulsive peak. A similar "late phase" is frequently detected in solar flares observed in the extreme-ultraviolet from longer hot coronal loops distinct from the impulsive flare structures. White-light emission has also been observed in off-limb solar flare loops. Here, we perform a suite of one-dimensional hydrodynamic loop simulations for M-dwarf flares inspired by these solar examples. Our results suggest that coronal plasma condensation following impulsive flare heating can yield high electron number density in the loop, allowing it to contribute significantly to the optical light curves via free-bound and free-free emission mechanisms. Our simulation results qualitatively agree with \textit{TESS} observations: the longer evolutionary time scale of coronal loops produces a distinct, secondary emission peak; its intensity increases with the injected flare energy. We argue that coronal plasma condensation is a possible mechanism for the \textit{TESS} late-phase flares. △ Less

Submitted 30 October, 2023; originally announced October 2023.

Comments: 31 pages, 13 figures, accepted for publication in ApJ

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

KIC 9845907: A $δ$ Scuti star with the first overtone as the dominant frequency and with many equidistant structures in its spectrum

Authors: Xiao-Ya Sun, Zhao-Yu Zuo, Tao-Zhi Yang, Antonio García Hernández

Abstract: In this paper, we present an analysis of the pulsating behavior of Kepler target KIC 9845907. Using the data from Kepler, we detected 85 significant frequencies, including the first overtone $f_{1}$ = 17.597 day$^{-1}$ as the dominant frequency, the non-radial independent frequency $f_{3}$ = 31.428 day$^{-1}$ ($\ell$=1), as well as two modulation terms $f_{m1}$ = 0.065 day$^{-1}$ and $f_{m2}$ = 1.… ▽ More In this paper, we present an analysis of the pulsating behavior of Kepler target KIC 9845907. Using the data from Kepler, we detected 85 significant frequencies, including the first overtone $f_{1}$ = 17.597 day$^{-1}$ as the dominant frequency, the non-radial independent frequency $f_{3}$ = 31.428 day$^{-1}$ ($\ell$=1), as well as two modulation terms $f_{m1}$ = 0.065 day$^{-1}$ and $f_{m2}$ = 1.693 day$^{-1}$. We found fourteen pairs of triplet structures with $f_{m1}$ or $f_{m2}$, four pairs of which can further form quintuplet structures. We note these are the most intriguing features discovered in this study and they were recognized for the first time in $δ$ Scuti stars. We discussed several possible explanations, i.e., beating, the Blazhko effect, combination mode hypothesis, nonlinear mode coupling, large separation, and stellar rotational splitting for these equidistant structures. Our asteroseismic models indicate this modulation with $f_{m1}$ might be related to the rotational splitting. The study of more $δ$ Scuti stars with triplet and/or quintuplet structures using high-precision space photometry would be helpful to further explore its origin. △ Less

Submitted 14 September, 2023; originally announced September 2023.

Comments: 28 pages, 10 figures, accepted for publication in ApJ

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

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

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

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

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

Studying X-ray spectra from large-scale jets of FR II radio galaxies: application of shear particle acceleration

Authors: Jia-Chun He, Xiao-Na Sun, Jie-Shuang Wang, Frank M. Rieger, Ruo-Yu Liu, En-Wei Liang

Abstract: Shear particle acceleration is a promising candidate for the origin of extended high-energy emission in extra-galactic jets. In this paper, we explore the applicability of a shear model to 24 X-ray knots in the large-scale jets of FR II radio galaxies, and study the jet properties by modeling the multi-wavelength spectral energy distributions (SEDs) in a leptonic framework including synchrotron an… ▽ More Shear particle acceleration is a promising candidate for the origin of extended high-energy emission in extra-galactic jets. In this paper, we explore the applicability of a shear model to 24 X-ray knots in the large-scale jets of FR II radio galaxies, and study the jet properties by modeling the multi-wavelength spectral energy distributions (SEDs) in a leptonic framework including synchrotron and inverse Compton - CMB processes. In order to improve spectral modelling, we analyze Fermi-LAT data for five sources and reanalyzed archival data of Chandra on 15 knots, exploring the radio to X-ray connection. We show that the X-ray SEDs of these knots can be satisfactorily modelled by synchrotron radiation from a second, shear-accelerated electron population reaching multi-TeV energies. The inferred flow speeds are compatible with large-scale jets being mildly relativistic. We explore two different shear flow profiles (i.e., linearly decreasing and power-law) and find that the required spine speeds differ only slightly, supporting the notion that for higher flow speeds the variations in particle spectral indices are less dependent on the presumed velocity profile. The derived magnetic field strengths are in the range of a few to ten microGauss, and the required power in non-thermal particles typically well below the Eddington constraint. Finally, the inferred parameters are used to constrain the potential of FR II jets as possible UHECR accelerators. △ Less

Submitted 22 August, 2023; originally announced August 2023.

Comments: 14 pages, 4 figures, 6 table, accepted for publication in MNRAS

Sub-second periodic radio oscillations in a microquasar

Authors: Pengfu Tian, Ping Zhang, Wei Wang, Pei Wang, Xiaohui Sun, Jifeng Liu, Bing Zhang, Zigao Dai, Feng Yuan, Shuangnan Zhang, Qingzhong Liu, Peng Jiang, Xuefeng Wu, Zheng Zheng, Jiashi Chen, Di Li, Zonghong Zhu, Zhichen Pan, Hengqian Gan, Xiao Chen, Na Sai

Abstract: Powerful relativistic jets are one of the ubiquitous features of accreting black holes in all scales. GRS 1915+105 is a well-known fast-spinning black-hole X-ray binary with a relativistic jet, termed as a ``microquasar'', as indicated by its superluminal motion of radio emission. It exhibits persistent x-ray activity over the last 30 years, with quasi-periodic oscillations of $\sim 1-10$ Hz and 3… ▽ More Powerful relativistic jets are one of the ubiquitous features of accreting black holes in all scales. GRS 1915+105 is a well-known fast-spinning black-hole X-ray binary with a relativistic jet, termed as a ``microquasar'', as indicated by its superluminal motion of radio emission. It exhibits persistent x-ray activity over the last 30 years, with quasi-periodic oscillations of $\sim 1-10$ Hz and 34 and 67 Hz in the x-ray band. These oscillations likely originate in the inner accretion disk, but other origins have been considered. Radio observations found variable light curves with quasi-periodic flares or oscillations with periods of $\sim 20-50$ minutes. Here we report two instances of $\sim$5 Hz transient periodic oscillation features from the source detected in the 1.05-1.45 GHz radio band that occurred in January 2021 and June 2022, respectively. Circular polarization was also observed during the oscillation phase. △ Less

Submitted 26 July, 2023; originally announced July 2023.

Comments: The author version of the article which will appear in Nature on 26 July 2023, 32 pages including the extended data. The online publication version can be found at the following URL: https://www.nature.com/articles/s41586-023-06336-6

Journal ref: 2023, Nature, 621, 271-275

Large Photospheric Doppler Shift in Solar Active Region 12673: I. Field-Aligned Flows

Authors: Jiayi Liu, Xudong Sun, Peter W. Schuck, Sarah A. Jaeggli, Brian T. Welsch, Carlos Quintero Noda

Abstract: Delta ($δ$) sunspots sometimes host fast photospheric flows along the central magnetic polarity inversion line (PIL). Here we study the strong Doppler shift signature in the central penumbral light bridge of solar active region NOAA 12673. Observations from the Helioseismic and Magnetic Imager (HMI) indicate highly sheared, strong magnetic fields. Large Doppler shifts up to 3.2 km s$^{-1}$ appeare… ▽ More Delta ($δ$) sunspots sometimes host fast photospheric flows along the central magnetic polarity inversion line (PIL). Here we study the strong Doppler shift signature in the central penumbral light bridge of solar active region NOAA 12673. Observations from the Helioseismic and Magnetic Imager (HMI) indicate highly sheared, strong magnetic fields. Large Doppler shifts up to 3.2 km s$^{-1}$ appeared during the formation of the light bridge and persisted for about 16 hours. A new velocity estimator, called DAVE4VMwDV, reveals fast converging and shearing motion along the PIL from HMI vector magnetograms, and recovers the observed Doppler signal much better than an old version of the algorithm. The inferred velocity vectors are largely (anti-)parallel to the inclined magnetic fields, suggesting that the observed Doppler shift contains significant contribution from the projected, field-aligned flows. High-resolution observations from the Hinode/Spectro-Polarimeter (SP) further exhibit a clear correlation between the Doppler velocity and the cosine of the magnetic inclination, which is in agreement with HMI results and consistent with a field-aligned flow of about 9.6 km s$^{-1}$. The complex Stokes profiles suggest significant gradients of physical variables along the line of sight. We discuss the implications on the $δ$-spot magnetic structure and the flow-driving mechanism. △ Less

Submitted 18 July, 2023; originally announced July 2023.

Comments: 17 pages, 14 figures. Accepted for publication in ApJ

Magnetar emergence in a peculiar gamma-ray burst from a compact star merger

Authors: H. Sun, C. -W. Wang, J. Yang, B. -B. Zhang, S. -L. Xiong, Y. -H. I. Yin, Y. Liu, Y. Li, W. -C. Xue, Z. Yan, C. Zhang, W. -J. Tan, H. -W. Pan, J. -C. Liu, H. -Q. Cheng, Y. -Q. Zhang, J. -W. Hu, C. Zheng, Z. -H. An, C. Cai, L. Hu, C. Jin, D. -Y. Li, X. -Q. Li, H. -Y. Liu , et al. (19 additional authors not shown)

Abstract: The central engine that powers gamma-ray bursts (GRBs), the most powerful explosions in the universe, is still not identified. Besides hyper-accreting black holes, rapidly spinning and highly magnetized neutron stars, known as millisecond magnetars, have been suggested to power both long and short GRBs. The presence of a magnetar engine following compact star mergers is of particular interest as i… ▽ More The central engine that powers gamma-ray bursts (GRBs), the most powerful explosions in the universe, is still not identified. Besides hyper-accreting black holes, rapidly spinning and highly magnetized neutron stars, known as millisecond magnetars, have been suggested to power both long and short GRBs. The presence of a magnetar engine following compact star mergers is of particular interest as it would provide essential constraints on the poorly understood equation of state for neutron stars. Indirect indications of a magnetar engine in these merger sources have been observed in the form of plateau features present in the X-ray afterglow light curves of some short GRBs. Additionally, some X-ray transients lacking gamma-ray bursts (GRB-less) have been identified as potential magnetar candidates originating from compact star mergers. Nevertheless, smoking gun evidence is still lacking for a magnetar engine in short GRBs, and the associated theoretical challenges have been addressed. Here we present a comprehensive analysis of the broad-band prompt emission data of a peculiar, very bright GRB 230307A. Despite its apparently long duration, the prompt emission and host galaxy properties point toward a compact star merger origin, being consistent with its association with a kilonova. More intriguingly, an extended X-ray emission component emerges as the $γ$-ray emission dies out, signifying the emergence of a magnetar central engine. We also identify an achromatic temporal break in the high-energy band during the prompt emission phase, which was never observed in previous bursts and reveals a narrow jet with half opening angle of approximately $3.4^\circ$. △ Less

Submitted 11 July, 2023; originally announced July 2023.

Comments: 44 pages, 10 figures, 5 tables

Beyond spectroscopy. II. Stellar parameters for over twenty million stars in the northern sky from SAGES DR1 and Gaia DR3

Authors: Yang Huang, Timothy C. Beers, Hai-Bo Yuan, Ke-Feng Tan, Wei Wang, Jie Zheng, Chun Li, Young Sun Lee, Hai-Ning Li, Jing-Kun Zhao, Xiang-Xiang Xue, Yu-Juan Liu, Hua-Wei Zhang, Xue-Ang Sun, Ji Li, Hong-Rui Gu, Christian Wolf, Christopher A. Onken, Ji-Feng Liu, Zhou Fan, Gang Zhao

Abstract: We present precise photometric estimates of stellar parameters, including effective temperature, metallicity, luminosity classification, distance, and stellar age, for nearly 26 million stars using the methodology developed in the first paper of this series, based on the stellar colors from the Stellar Abundances and Galactic Evolution Survey (SAGES) DR1 and Gaia EDR3. The optimal design of stella… ▽ More We present precise photometric estimates of stellar parameters, including effective temperature, metallicity, luminosity classification, distance, and stellar age, for nearly 26 million stars using the methodology developed in the first paper of this series, based on the stellar colors from the Stellar Abundances and Galactic Evolution Survey (SAGES) DR1 and Gaia EDR3. The optimal design of stellar-parameter sensitive $uv$ filters by SAGES has enabled us to determine photometric-metallicity estimates down to $-3.5$, similar to our previous results with the SkyMapper Southern Survey (SMSS), yielding a large sample of over five million metal-poor (MP; [Fe/H]$\le -1.0$) stars and nearly one million very metal-poor (VMP; [Fe/H]$\le -2.0$) stars. The typical precision is around $0.1$ dex for both dwarf and giant stars with [Fe/H]$>-1.0$, and 0.15-0.25/0.3-0.4 dex for dwarf/giant stars with [Fe/H]$<-1.0$. Using the precise parallax measurements and stellar colors from Gaia, effective temperature, luminosity classification, distance and stellar age are further derived for our sample stars. This huge data set in the Northern sky from SAGES, together with similar data in the Southern sky from SMSS, will greatly advance our understanding of the Milky Way, in particular its formation and evolution. △ Less

Submitted 10 July, 2023; originally announced July 2023.

Comments: 14 pages, 14 figures, 3 tables, accepted by ApJ. arXiv admin note: text overlap with arXiv:2104.14154

A Catalogue of Radio Supernova Remnants and Candidate Supernova Remnants in the EMU/POSSUM Galactic Pilot Field

Authors: Brianna D. Ball, Roland Kothes, Erik Rosolowsky, Jennifer West, Werner Becker, Miroslav D. Filipović, B. M. Gaensler, Andrew M. Hopkins, Bärbel Koribalski, Tom Landecker, Denis Leahy, Joshua Marvil, Xiaohui Sun, Filomena Bufano, Ettore Carretti, Adriano Ingallinera, Cameron L. Van Eck, Tony Willis

Abstract: We use data from the pilot observations of the EMU/POSSUM surveys to study the "missing supernova remnant (SNR) problem", the discrepancy between the number of Galactic SNRs that have been observed and the number that are estimated to exist. The Evolutionary Map of the Universe (EMU) and the Polarization Sky Survey of the Universe's Magnetism (POSSUM) are radio sky surveys that are conducted using… ▽ More We use data from the pilot observations of the EMU/POSSUM surveys to study the "missing supernova remnant (SNR) problem", the discrepancy between the number of Galactic SNRs that have been observed and the number that are estimated to exist. The Evolutionary Map of the Universe (EMU) and the Polarization Sky Survey of the Universe's Magnetism (POSSUM) are radio sky surveys that are conducted using the Australian Square Kilometre Array Pathfinder (ASKAP). We report on the properties of 7 known SNRs in the joint Galactic pilot field, with an approximate longitude and latitude of 323$^\circ\leq$ l $\leq$ 330$^\circ$ and -4$^\circ\leq$ b $\leq$ 2$^\circ$ respectively, and identify 21 SNR candidates. Of these, 4 have been previously identified as SNR candidates, 3 were previously listed as a single SNR, 13 have not been previously studied, and 1 has been studied in the infrared. These are the first discoveries of Galactic SNR candidates with EMU/POSSUM and, if confirmed, they will increase the SNR density in this field by a factor of 4. By comparing our SNR candidates to the known Galactic SNR population, we demonstrate that many of these sources were likely missed in previous surveys due to their small angular size and/or low surface brightness. We suspect that there are SNRs in this field that remain undetected due to limitations set by the local background and confusion with other radio sources. The results of this paper demonstrate the potential of the full EMU/POSSUM surveys to uncover more of the missing Galactic SNR population. △ Less

Submitted 4 July, 2023; originally announced July 2023.

The First GECAM Observation Results on Terrestrial Gamma-ray Flashes and Terrestrial Electron Beams

Authors: Y. Zhao, J. C. Liu, S. L. Xiong, W. C. Xue, Q. B. Yi, G. P. Lu, W. Xu, F. C. Lyu, J. C. Sun, W. X. Peng, C. Zheng, Y. Q. Zhang, C. Cai, S. Xiao, S. L. Xie, C. W. Wang, W. J. Tan, Z. H. An, G. Chen, Y. Q. Du, Y. Huang, M. Gao, K. Gong, D. Y. Guo, J. J. He , et al. (37 additional authors not shown)

Abstract: Gravitational-wave high-energy Electromagnetic Counterpart All-sky Monitor (GECAM) is a space-borne instrument dedicated to monitoring high-energy transients, including Terrestrial Gamma-ray Flashes (TGFs) and Terrestrial Electron Beams (TEBs). We implemented a TGF/TEB search algorithm for GECAM, with which 147 bright TGFs, 2 typical TEBs and 2 special TEB-like events are identified during an effe… ▽ More Gravitational-wave high-energy Electromagnetic Counterpart All-sky Monitor (GECAM) is a space-borne instrument dedicated to monitoring high-energy transients, including Terrestrial Gamma-ray Flashes (TGFs) and Terrestrial Electron Beams (TEBs). We implemented a TGF/TEB search algorithm for GECAM, with which 147 bright TGFs, 2 typical TEBs and 2 special TEB-like events are identified during an effective observation time of $\sim$9 months. We show that, with gamma-ray and charged particle detectors, GECAM can effectively identify and distinguish TGFs and TEBs, and measure their temporal and spectral properties in detail. A very high TGF-lightning association rate of $\sim$80\% is obtained between GECAM and GLD360 in east Asia region. △ Less

Submitted 17 June, 2023; originally announced June 2023.

Comments: The paper was accepted by Geophysical Research Letters on June 16th, 2023

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

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

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

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

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

Journal ref: JCAP 09 (2023) 001

Radio Sources Segmentation and Classification with Deep Learning

Authors: Baoqiang Lao, Sumit Jaiswal, Zhen Zhao, Leping Lin, Junyi Wang, Xiaohui Sun, Shengli Qin

Abstract: Modern large radio continuum surveys have high sensitivity and resolution, and can resolve previously undetected extended and diffuse emissions, which brings great challenges for the detection and morphological classification of extended sources. We present HeTu-v2, a deep learning-based source detector that uses the combined networks of Mask Region-based Convolutional Neural Networks (Mask R-CNN)… ▽ More Modern large radio continuum surveys have high sensitivity and resolution, and can resolve previously undetected extended and diffuse emissions, which brings great challenges for the detection and morphological classification of extended sources. We present HeTu-v2, a deep learning-based source detector that uses the combined networks of Mask Region-based Convolutional Neural Networks (Mask R-CNN) and a Transformer block to achieve high-quality radio sources segmentation and classification. The sources are classified into 5 categories: Compact or point-like sources (CS), Fanaroff-Riley Type I (FRI), Fanaroff-Riley Type II (FRII), Head-Tail (HT), and Core-Jet (CJ) sources. HeTu-v2 has been trained and validated with the data from the Faint Images of the Radio Sky at Twenty-one centimeters (FIRST). We found that HeTu-v2 has a high accuracy with a mean average precision ($AP_{\rm @50:5:95}$) of 77.8%, which is 15.6 points and 11.3 points higher than that of HeTu-v1 and the original Mask R-CNN respectively. We produced a FIRST morphological catalog (FIRST-HeTu) using HeTu-v2, which contains 835,435 sources and achieves 98.6% of completeness and up to 98.5% of accuracy compared to the latest 2014 data release of the FIRST survey. HeTu-v2 could also be employed for other astronomical tasks like building sky models, associating radio components, and classifying radio galaxies. △ Less

Submitted 5 June, 2023; v1 submitted 2 June, 2023; originally announced June 2023.

Comments: 20 pages, 15 figures, Accepted by Astronomy and Computing

Journal ref: 2023, Astronomy and Computing Volume 44, Pages 100728

FAST polarization mapping of the SNR VRO 42.05.01

Authors: Li Xiao, Ming Zhu, Xiao-Hui Sun, Peng Jiang, Chun Sun

Abstract: We have obtained the polarization data cube of the VRO 42.05.01 supernova remnant at 1240 MHz using the Five-hundred-meter Aperture Spherical radio Telescope (FAST). Three-dimensional Faraday Synthesis is applied to the FAST data to derive the Faraday depth spectrum. The peak Faraday depth map shows a large area of enhanced foreground RM of ~60 rad m-2 extending along the remnant's "wing" section,… ▽ More We have obtained the polarization data cube of the VRO 42.05.01 supernova remnant at 1240 MHz using the Five-hundred-meter Aperture Spherical radio Telescope (FAST). Three-dimensional Faraday Synthesis is applied to the FAST data to derive the Faraday depth spectrum. The peak Faraday depth map shows a large area of enhanced foreground RM of ~60 rad m-2 extending along the remnant's "wing" section, which coincides with a large-scale HI shell at -20 km/s. The two depolarization patches within the "wing" region with RM of 97 rad m-2 and 55 rad m-2 coincide with two HI structures in the HI shell. Faraday screen model fitting on the Canadian Galactic Plane Survey (CGPS) 1420 MHz full-scale polarization data reveals a distance of 0.7-0.8d_{SNR} in front of the SNR with enhanced regular magnetic field there. The highly piled-up magnetic field indicates that the HI shell at -20 km/s could originate from an old evolved SNR. △ Less

Submitted 2 June, 2023; originally announced June 2023.

Comments: 9 pages, 8 figures, accepted by ApJ

Fine structures of radio bursts from flare star AD Leo with FAST observations

Authors: Jiale Zhang, Hui Tian, Philippe Zarka, Corentin K. Louis, Hongpeng Lu, Dongyang Gao, Xiaohui Sun, Sijie Yu, Bin Chen, Xin Cheng, Ke Wang

Abstract: Radio bursts from nearby active M-dwarfs have been frequently reported and extensively studied in solar or planetary paradigms. Whereas, their sub-structures or fine structures remain rarely explored despite their potential significance in diagnosing the plasma and magnetic field properties of the star. Such studies in the past have been limited by the sensitivity of radio telescopes. Here we repo… ▽ More Radio bursts from nearby active M-dwarfs have been frequently reported and extensively studied in solar or planetary paradigms. Whereas, their sub-structures or fine structures remain rarely explored despite their potential significance in diagnosing the plasma and magnetic field properties of the star. Such studies in the past have been limited by the sensitivity of radio telescopes. Here we report the inspiring results from the high time-resolution observations of a known flare star AD Leo with the Five-hundred-meter Aperture Spherical radio Telescope (FAST). We detected many radio bursts in the two days of observations with fine structures in the form of numerous millisecond-scale sub-bursts. Sub-bursts on the first day display stripe-like shapes with nearly uniform frequency drift rates, which are possibly stellar analogs to Jovian S-bursts. Sub-bursts on the second day, however, reveal a different blob-like shape with random occurrence patterns and are akin to solar radio spikes. The new observational results suggest that the intense emission from AD Leo is driven by electron cyclotron maser instability which may be related to stellar flares or interactions with a planetary companion. △ Less

Submitted 1 June, 2023; originally announced June 2023.

Comments: 25 pages, 12 figures, accepted for publication in ApJ

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

Polarization Criterion in Targeted SETI Observation

Authors: Jian-Kang Li, Yu Chen, Zhen-Zhao Tao, Xiao-Hang Luan, Tong-Jie Zhang, Bo-Lun Huang, Xiao-Hui Sun, Vishal Gajjar

Abstract: In this paper, we propose a novel method for distinguishing extraterrestrial intelligence (ETI) signals from radio frequency interference (RFI) by leveraging polarization features. We exploit the sinusoidal variation of the linearly polarized components of Stokes parameters with the parallactic angle as a characteristic signature of ETI signals, while such linearly polarized components remain rela… ▽ More In this paper, we propose a novel method for distinguishing extraterrestrial intelligence (ETI) signals from radio frequency interference (RFI) by leveraging polarization features. We exploit the sinusoidal variation of the linearly polarized components of Stokes parameters with the parallactic angle as a characteristic signature of ETI signals, while such linearly polarized components remain relatively stable for terrestrial RFI. Typically, a minimum of 4-8 hours of observation time is required to detect these sinusoidal variations. The polarization approach in the search for extraterrestrial intelligence (SETI) also enables us to study the radio stellar bursts emitted by M-type stars as ancillary science, which is relevant to assessing the habitability of exoplanets. Compared to the frequency drift method, the polarization method effectively reduces the required observation time for signal identification while improving the signal identification process. △ Less

Submitted 7 December, 2023; v1 submitted 25 May, 2023; originally announced May 2023.

Comments: 8 pages, 2 figures, published in The Astronomical Journal. Using polarization features as the criterion in SETI research

Journal ref: The Astronomical Journal, 167:8 (5pp), 2024

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)

Limitations and biases in the retrieval of the polar magnetic field I: the role of the magnetic filling factor in Milne-Eddington inversions of simulated Hinode/SP data

Authors: Rebecca Centeno, Ivan Milić, Matthias Rempel, Nariaki V. Nitta, Xudong Sun

Abstract: We study the extent to which Milne-Eddington inversions are able to retrieve and characterize the magnetic landscape of the solar poles from observations by the spectropolarimeter onboard Hinode. In particular, we evaluate whether a variable magnetic filling factor is an adequate modeling technique for retrieving the intrinsic magnetic properties from every pixel in the polar field of view. We fir… ▽ More We study the extent to which Milne-Eddington inversions are able to retrieve and characterize the magnetic landscape of the solar poles from observations by the spectropolarimeter onboard Hinode. In particular, we evaluate whether a variable magnetic filling factor is an adequate modeling technique for retrieving the intrinsic magnetic properties from every pixel in the polar field of view. We first generate synthetic spectra emerging from a numerical simulation of a "plage" region at an inclined line of sight of 65$^{\circ}$, and degrade the data to emulate real observations. Then, we invert the synthetic spectra with two Milne-Eddington inversion codes that feature different treatments of the magnetic filling factor, and relate the retrieved magnetic quantities back to their original values in the simulation cube. We find that while the apparent retrieved magnetic properties map well the spatially-degraded simulation, the intrinsic magnetic quantities bear little relation to the magnetic field at the native resolution of the simulation. We discuss the systematic biases caused by line-of-sight foreshortening, spatial degradation, photon noise and modeling assumptions embedded in the inversion algorithm. △ Less

Submitted 1 May, 2023; originally announced May 2023.

Comments: Accepted for publication in ApJ

Plasma lensing with magnetic field and a small correction to the Faraday rotation measurement

Authors: Xinzhong Er, Ue-Li Pen, Xiaohui Sun, Dongzi Li

Abstract: Plasma lensing displays interesting characteristics that set it apart from gravitational lensing. The magnetised medium induces birefringence in the two polarisation modes. As the lensing deflection grows stronger, e.g. when images form near the critical curve, the geometric delay of the signal can cause rotation in linear polarisation, in addition to Faraday rotation. This rotation has a frequenc… ▽ More Plasma lensing displays interesting characteristics that set it apart from gravitational lensing. The magnetised medium induces birefringence in the two polarisation modes. As the lensing deflection grows stronger, e.g. when images form near the critical curve, the geometric delay of the signal can cause rotation in linear polarisation, in addition to Faraday rotation. This rotation has a frequency dependence to the power of four. We study the geometric rotation of the lensed image in a Gaussian density model and find that it is necessary to take into account the geometric rotation when estimating magnetised media, especially in the under-dense lens. At frequencies of $\sim 1$ GHz or lower, the geometric rotation can dominate. We simulate the flux of lensed images and find that when the image forms near the lensing critical curve, the birefringence can convert the linear polarisation and un-polarisation pulse into a circular mode. The lensing magnification has the potential to increase the probability of detecting such events. △ Less

Submitted 25 April, 2023; originally announced April 2023.

Comments: MNRAS, 7 pages, comments welcome

The Magnetohydrodynamic-Particle-In-Cell Module in Athena++: Implementation and Code Tests

Authors: Xiaochen Sun, Xue-Ning Bai

Abstract: We present a new magnetohydrodynamic-particle-in-cell (MHD-PIC) code integrated into the Athena++ framework. It treats energetic particles as in conventional PIC codes while the rest of thermal plasmas are treated as background fluid described by MHD, thus primarily targeting at multi-scale astrophysical problems involving the kinetic physics of the cosmic-rays (CRs). The code is optimized toward… ▽ More We present a new magnetohydrodynamic-particle-in-cell (MHD-PIC) code integrated into the Athena++ framework. It treats energetic particles as in conventional PIC codes while the rest of thermal plasmas are treated as background fluid described by MHD, thus primarily targeting at multi-scale astrophysical problems involving the kinetic physics of the cosmic-rays (CRs). The code is optimized toward efficient vectorization in interpolation and particle deposits, with excellent parallel scaling. The code is also compatible with static/adaptive mesh refinement, with dynamic load balancing to further enhance multi-scale simulations. In addition, we have implemented a compressing/expanding box framework which allows adiabatic driving of CR pressure anisotropy, as well as the $δf$ method that can dramatically reduce Poisson noise in problems where distribution function $f$ is only expected to slightly deviate from the background. The code performance is demonstrated over a series of benchmark test problems including particle acceleration in non-relativistic parallel shocks. In particular, we reproduce the linear growth of the CR gyro-resonant (streaming and pressure anisotropy) instabilities, under both the periodic and expanding/compressing box setting. We anticipate the code to open up the avenue for a wide range of astrophysical and plasma physics applications. △ Less

Submitted 20 April, 2023; originally announced April 2023.

Comments: 20 pages, 19 figures, submitted to MNRAS

Removal Of Active Region Inflows Reveals a Weak Solar Cycle Scale Trend In Near-surface Meridional Flow

Authors: Sushant S. Mahajan, Xudong Sun, Junwei Zhao

Abstract: Using time-distance local helioseismology flow maps within 1 Mm of the solar photosphere, we detect inflows toward activity belts that contribute to solar cycle scale variations in near-surface meridional flow. These inflows stretch out as far as 30 degrees away from active region centroids. If active region neighborhoods are excluded, the solar cycle scale variation in background meridional flow… ▽ More Using time-distance local helioseismology flow maps within 1 Mm of the solar photosphere, we detect inflows toward activity belts that contribute to solar cycle scale variations in near-surface meridional flow. These inflows stretch out as far as 30 degrees away from active region centroids. If active region neighborhoods are excluded, the solar cycle scale variation in background meridional flow diminishes to below 2~m~s$^{-1}$, but still shows systematic variations in the absence of active regions between Sunspot Cycles 24 and 25. We, therefore, propose that the near-surface meridional flow is a three component flow made up of: a constant baseline flow profile that can be derived from quiet Sun regions, variations due to inflows around active regions, and solar cycle scale variation of the order of 2~m~s$^{-1}$. Torsional oscillation, on the other hand, is found to be a global phenomenon i.e. exclusion of active region neighborhoods does not affect its magnitude or phase significantly. This non-variation of torsional oscillation with distance away from active regions and the three-component breakdown of the near-surface meridional flow serve as vital constraints for solar dynamo models and surface flux transport simulations. △ Less

Submitted 4 April, 2023; originally announced April 2023.

Comments: 14 pages, 9 figures, accepted for publication in the Astrophysical Journal

Probing darK Matter Using free leptONs: PKMUON

Authors: Alim Ruzi, Chen Zhou, Xiaohu Sun, Dayong Wang, Siguang Wang, Yong Ban, Yajun Mao, Qite Li, Qiang Li

Abstract: We propose a new method to detect sub-GeV dark matter, through their scatterings from free leptons and the resulting kinematic shifts. Specially, such an experiment can detect dark matter interacting solely with muons. The experiment proposed here is to directly probe muon-philic dark matter, in a model-independent way. Its complementarity with the muon on target proposal, is similar to, e.g. XENO… ▽ More We propose a new method to detect sub-GeV dark matter, through their scatterings from free leptons and the resulting kinematic shifts. Specially, such an experiment can detect dark matter interacting solely with muons. The experiment proposed here is to directly probe muon-philic dark matter, in a model-independent way. Its complementarity with the muon on target proposal, is similar to, e.g. XENON/PandaX and ATLAS/CMS on dark matter searches. Moreover, our proposal can work better for relatively heavy dark matter such as in the sub-GeV region. We start with a small device of a size around 0.1 to 1 meter, using atmospheric muons to set up a prototype. Within only one year of operation, the sensitivity on cross section of dark matter scattering with muons can already reach $σ_D\sim 10^{-19 (-20,\,-18)}\rm{cm}^{2}$ for a dark mater $\rm{M_D}=100\, (10,\,1000)$ MeV. We can then interface the device with a high intensity muon beam of $10^{12}$/bunch. Within one year, the sensitivity can reach $σ_D\sim 10^{-27 (-28,\,-26)}\rm{cm}^{2}$ for $\rm{M_D}=100\, (10,\,1000)$ MeV. △ Less

Submitted 10 October, 2023; v1 submitted 31 March, 2023; originally announced March 2023.

Comments: Accepted version