Study on the Performance of the GRANDProto300 Particle Detector Array by Simulation

The Giant Radio Array for Neutrino Detection(GRAND)is a proposed large-scale observatory designed to detect cosmic rays,gamma-rays,and neutrinos with energies exceeding 100 Pe V.The GRANDProto300 experiment is proposed as the early stage of the GRAND project,consisting of a hybrid array of radio antennas and scintillator detectors.The latter,as a mature and traditional detector,is used to cross-check the nature of the candidate events selected from radio observations.In this study,we developed a simulation software called G4GRANDProto300,based on the Geant4 software package,to optimize the spacing of the scintillator detector array and to investigate its effective area.The analysis was conducted at various zenith angles under different detector spacings,including 300,500,600,700,and 900 m.Our results indicate that,for large zenith angles used to search for cosmic-ray in the GRAND project,the optimized effective area is with a detector spacing of 500 m.The G4GRANDProto300 software that we developed could be used to further optimize the layout of the particle detector array in future work.

Analysis of Bright Source Hardness Ratios in the 4 yr Insight-HXMT Galactic Plane Scanning Survey Catalog

We conduct a statistical analysis of the hardness ratio(HR)for bright sources in the 4 yr Galactic Plane Scanning Survey catalog of Insight-HXMT.Depending on the stable(variable)flux F_(s)(F_(v))or spectrum S_(s)(S_(v))of each source,the bright sources are classified into three groups:F_(v)&S_(v),F_(v)&S_(s),and F_(s)&_(s).Our study of the HR characteristics in different types of sources reveals that accretion-powered neutron star(NS)low-mass X-ray binaries(LMXBs)exhibit softer energy spectra than NS high-mass X-ray binaries(HMXBs),but harder energy spectra than black hole binaries in most cases.This difference is probably due to their different magnetic field strengths.Additionally,Fv&Sv LMXBs tend to be harder than Fv&Ss LMXBs below 7 keV,while the opposite is true for HMXBs.Our results suggest that LMXBs may dominate unclassified sources,and NS binaries are likely to be the primary type of X-ray binaries with ambiguous compact stars.By comparing the HR of transient sources in their outburst and low-flux states,it is found that the averaged HR of four sources in the two states are roughly comparable within uncertainties.We also investigate the spatial properties of the three groups and find that Fv&Sv sources are mainly located in the longitude of-20°<l<9°,Fv&Ss sources cross the Galactic Plane,and Fs&Ss sources are predominantly concentrated in 19°<l<42°.In addition,analyzing the HR spatial distributions shows the absorption of soft X-rays(primarily below 2 keV)in the Galactic Plane.

The Precursor of GRB211211A:A Tide-induced Giant Quake?

The equilibrium configuration of a solid strange star in the final inspiral phase with another compact object is generally discussed,and the starquake-related issue is revisited,for a special purpose to understand the precursor emission of binary compact star merger events(e.g.,that of GRB211211A).As the binary system inspirals inward due to gravitational wave radiation,the ellipticity of the solid strangeon star increases due to the growing tidal field of its compact companion.Elastic energy is hence accumulated during the inspiral stage which might trigger a starquake before the merger when the energy exceeds a critical value.The energy released during such starquakes is calculated and compared to the precursor observation of GRB211211 A.The result shows that the energy might be insufficient for binary strangeon-star case unless the entire solid strangeon star shatters,and hence favors a black hole-strangeon star scenario for GRB211211A.The timescale of the precursor as well as the frequency of the observed quasi-periodic-oscillation have also been discussed in the starquake model.

The Influence of the Sun and Moon on the Observation of Very High Energy Gamma-ray Sources Using EAS Arrays

With great advance of ground-based extensive air shower arrays,such as LHAASO and HAWC,many very high energy(VHE)gamma-ray sources have been discovered and are being monitored regardless of the day and the night.Hence,the Sun and Moon would have some impacts on the observation of gamma-ray sources,which have not been taken into account in previous analysis.In this paper,the influence of the Sun and Moon on the observation of very high energy gamma-ray sources when they are near the line of sight of the Sun or Moon is estimated.The tracks of all the known VHE sources are scanned and several VHE sources are found to be very close to the line of sight of the Sun or Moon during some period.The absorption of very high energy gamma rays by sunlight is estimated with detailed method and some useful conclusions are achieved.The main influence is the block of the Sun and Moon on gamma rays and the shadow on the cosmic ray background.The influence is investigated considering the detector angular resolution and some strategies on data analysis are proposed to avoid the underestimation of the gamma-ray emission.

The energy response of LaBr_(3)(Ce),LaBr_(3)(Ce,Sr),and NaI(Tl)crystals for GECAM

The GECAM series of satellites utilizes LaBr_(3)(Ce),LaBr_(3)(Ce,Sr),and NaI(Tl)crystals as sensitive materials for gamma-ray detectors(GRDs).To investigate the nonlinearity in the detection of low-energy gamma rays and address the errors in the calibration of the E-C relationship,comprehensive tests and comparative studies of the three aforementioned crystals were conducted using Compton electrons,radioactive sources,and mono-energetic X-rays.The nonlinearity test results of the Compton electrons and X-rays demonstrated substantial differences,with all three crystals presenting a higher nonlinearity for X/-rays than for Compton electrons.Despite the LaBr_(3)(Ce)and LaBr_(3)(Ce,Sr)crystals having higher absolute light yields,they exhibited a noticeable nonlinear decrease in the light yield,especially at energies below 400 keV.The NaI(Tl)crystal demonstrated an"excess"light output in the 6-200 keV range,reaching a maximum"excess"of 9.2%at 30 keV in the 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,we explored the underlying causes of the observed nonlinearity in these crystals.This study not only elucidates the detector responses of GECAM,but also initiates a comprehensive investigation of the nonlinearity of domestically produced lanthanum bromide and sodium iodide crystals.

Weak Merging Scenario of CLASH Cluster A209

We study the structural and dynamical properties of A209 based on Chandra and XMM-Newton observations.We obtain detailed temperature,pressure,and entropy maps with the contour binning method,and find a hot region in the NW direction.The X-ray brightness residual map and corresponding temperature profiles reveal a possible shock front in the NW direction and a cold front feature in the SE direction.Combined with the galaxy luminosity density map we propose a weak merger scenario.A young sub-cluster passing from the SE to NW direction could explain the optical subpeak,the intracluster medium temperature map,the X-ray surface brightness excess,and the X-ray peak offset together.

Temporal and Spectral Characteristics of Persistent Emission and Special Bursts of Magnetar SGR J1935+2154 Based on Insight-HXMT

In October 2022,the magnetar SGR J1935+2154 entered the active outburst state.During the episode,the InsightHXMT satellite carried out a long observation that lasted for 20 days.More than 300 bursts were detected,and a certain amount of persistent radiation signals were also accumulated.This paper mainly introduces the results of persistent radiation profile folding and period search based on Insight-HXMT data.At the same time,the burst phase distribution characteristics,spectral lag results of burst,the spectral characteristics of zero-lag bursts and the time-resolved spectral evolution characteristics of high-flux bursts are reported.We found that there is no significant delay feature during different energy bands for the bursts of SGR J1935+2154.The observed zero-lag burst does not have a unique spectrum.The time-resolved spectrum of the individual burst has consistent spectral types and spectral parameters at different time periods of the burst.We also find that the burst number phase distribution and the burst photon phase distribution have the same tendency to concentrate in specific regions of the persistent emission profile.

The Design of GECAM Scientific Ground Segment

The Gravitational wave burst high-energy Electromagnetic Counterpart All-sky Monitor(GECAM)is a dedicated mission for monitoring high-energy transients.Here we report the design of the GECAM Scientific Ground Segment(GSGS)in terms of the scientific requirements,including the architecture,the external interfaces,the main function,and workflow.Judging from the analysis and verification results during the commissioning phase,the GSGS functions well and is able to monitor the status of the payloads,adjust the parameters,develop the scientific observation plans,generate the scientific data products,analyze the data,etc.Thus,the on-orbit operation and scientific researches of GECAM are guaranteed.

Timing and Spectral Analysis of the Black Hole X-Ray Binary MAXI J1803-298 with Insight-HXMT Data

We present a comprehensive analysis of the 2021 outburst of MAXI J1803–298 utilizing observations of the Insight-Hard X-ray Modulation Telescope(Insight-HXMT)spanning from the low hard state to the high soft state.Within the Insight-HXMT data set,compared to the previous work,we identify a more prolonged presence of typeC quasi-periodic oscillations(QPOs)with centroid frequencies ranging from~0.16 to 6.3 Hz,which present correlations with the hardness ratio and the photon index of the Comptonized component.For QPO frequencies less than~2 Hz,the QPO phase lags are hard(photons of 10–19 keV arrive later than those of 1–4 keV),while at higher frequencies,the lags become soft at and above~4 Hz.Furthermore,the spectra in all Insight-HXMT observations consist of a multi-color blackbody component and a Comptonized component,as commonly observed in classical black hole X-ray binaries.We analyze state transitions and the evolution of accretion geometry in this work.The fitted inner disk radius increases abnormally during the low hard state,hypothesized to result from the corona condensing onto the inner disk.Additionally,two significant drops in flux are observed during the soft intermediate state,maybe implying changes in the corona/jet and the disk,respectively.

Progress of GECAM Observation Research

Gravitational wave high-energy Electromagnetic Counterpart All-sky Monitor(GECAM)is a constellation with four instruments(launch date):GECAM-A/B(10 December 2020),GECAM-C(27 July 2022)and GECAM-D(13 March 2024),which are dedicated to monitoring gamma-ray transients in all-sky.The primary science objectives of GECAM include Gamma-Ray Bursts(GRBs),Soft Gamma-ray Repeaters(SGRs),high energy counterparts of Gravitation Wave(GW)and Fast Radio Burst(FRB),Solar Flares(SFLs),as well as Terrestrial Gamma-ray Flashes(TGFs)and Terrestrial Electron Beams(TEBs).A series of observations and research have been made since the launch of GECAM-A/B.GECAM observations provide new insights into these highenergy transients,demonstrating the unique role of GECAM in the“multi-wavelength,multi-messenger”era.

Searching for the Highest Energy of Pulsation and Critical Luminosity of Swift J0243.6+6124 Observed by Insight-HXMT

Owing to the broad energy coverage of Insight-HXMT in the hard X-ray band,we detected the highest energy of pulsation exceeding 200 keV around the 2017–2018 outburst peak of the first Galactic pulsating ultraluminous X-ray source (PULX) Swift J0243.6+6124,which is the highest energy detected from PULXs to date.We also obtained the highest energy of pulsation of every exposure during the outburst in 2017–2018,and found the highest energy is roughly positively correlated with luminosity.Using our newly developed method,we identified the critical luminosity being 4×10^(38)erg s^(-1) when the main peaks of the low and high energy pulse profiles became aligned,which separates the fan-beam dominated and pencil-beam dominated accretion regimes.Above the critical luminosity,the phase of the main peak shifted gradually from 0.5 to 0.8 until the outburst peak in al energy bands is reached,which is in agreement with the phase shift found previously at low energies.Our result is consistent with what is derived from spectral analysis.

Polarization Study of Swift J151857.0–572147 with IXPE Observation

We present an analysis of the Imaging X-ray Polarimetry Explorer observation from a newly discovered transient source:Swift J151857.0-572147.The obtained polarization degree(PD)and angle are 0.3%±0.3% and -24°±26°respectively in 2-8 keV within 68%confidence level errors,and polarization results are below MDP99in all energy bins,with the upper limit on PD of 0.8% in the 2-8 keV energy range.No quasi-periodic oscillations(QPOs)are detected in this observation.The polarization and QPO analyses support the hypothesis that the source was in the high soft state,and the results are consistent with predictions for a thin accretion disk model.

Application of deep learning methods combined with physical background in wide field of view imaging atmospheric Cherenkov telescopes

The High Altitude Detection of Astronomical Radiation(HADAR)experiment,which was constructed in Tibet,China,combines the wide-angle advantages of traditional EAS array detectors with the high-sensitivity advantages of focused Cherenkov detectors.Its objective is to observe transient sources such as gamma-ray bursts and the counterparts of gravitational waves.This study aims to utilize the latest AI technology to enhance the sensitivity of HADAR experiments.Training datasets and models with distinctive creativity were constructed by incorporating the relevant physical theories for various applications.These models can determine the type,energy,and direction of the incident particles after careful design.We obtained a background identification accuracy of 98.6%,a relative energy reconstruction error of 10.0%,and an angular resolution of 0.22°in a test dataset at 10 TeV.These findings demonstrate the significant potential for enhancing the precision and dependability of detector data analysis in astrophysical research.By using deep learning techniques,the HADAR experiment’s observational sensitivity to the Crab Nebula has surpassed that of MAGIC and H.E.S.S.at energies below 0.5 TeV and remains competitive with conventional narrow-field Cherenkov telescopes at higher energies.In addition,our experiment offers a new approach for dealing with strongly connected,scattered data.

The In-Flight Realtime Trigger and Localization Software of GECAM

Realtime trigger and localization of bursts are the key functions of GECAM,an all-sky gamma-ray monitor launched on 2020 December 10.We developed a multifunctional trigger and localization software operating in the CPU of the GECAM Electronic Box.This onboard software has the following features:high trigger efficiency for real celestial bursts with a suppression of false triggers caused by charged particle bursts and background fluctuation,dedicated localization algorithm optimized for both short and long bursts,and low time latency of the trigger information which is downlinked through the Global Short Message Communication service of the global BeiDou navigation system.This paper provides a detailed description of the design and development of the trigger and localization software system for GECAM.It covers the general design,workflow,the main functions,and the algorithms used in the system.The paper also includes on-ground trigger tests using simulated gamma-ray bursts generated by a dedicated X-ray tube,as well as an overview of the performance for real celestial bursts during its in-orbit operation.

Spin Evolution of the Magnetar SGR J1935+2154

Fast radio bursts(FRBs)are short pulses observed in radio frequencies usually originating from cosmological distances.The discovery of FRB 200428 and its X-ray counterpart from the Galactic magnetar SGR J1935+2154suggests that at least some FRBs can be generated by magnetars.However,the majority of X-ray bursts from magnetars are not associated with radio emission.The fact that only in rare cases can an FRB be generated raises the question regarding the special triggering mechanism of FRBs.Here we report long time spin evolution of SGR J1935+2154 until the end of 2022.According to v and v,the spin evolution of SGR J1935+2154 could be divided into two stages.The first stage evolves relatively steady evolution until 2020 April 27.After the burst activity in2020,the spin of SGR J1935+2154 shows strong variations,especially for v.After the burst activity in 2022October,a new spin-down glitch with△v/v=(-7.2±0.6)×10^(-6)is detected around MJD 59876,which is the second event in SGR J1935+2154.At the end,spin frequency and pulse profile do not show variations around the time of FRB 200428 and radio bursts 221014 and 221021,which supply strong clues to constrain the trigger mechanism of FRBs or radio bursts.

A Localization Method of High Energy Transients for All-sky Gamma-ray Monitor

Fast and reliable localization of high-energy transients is crucial for characterizing the burst properties and guiding the follow-up observations.Localization based on the relative counts of different detectors has been widely used for all-sky gamma-ray monitors.There are two major methods for this count distribution localization:χ^(2)minimization method and the Bayesian method.Here we propose a modified Bayesian method that could take advantage of both the accuracy of the Bayesian method and the simplicity of the χ^(2)method.With comprehensive simulations,we find that our Bayesian method with Poisson likelihood is generally more applicable for various bursts than the χ^(2)method,especially for weak bursts.We further proposed a location-spectrum iteration approach based on the Bayesian inference,which could alleviate the problems caused by the spectral difference between the burst and location templates.Our method is very suitable for scenarios with limited computation resources or timesensitive applications,such as in-flight localization software,and low-latency localization for rapidly follow-up observations.

Burst Phase Distribution of SGR J1935+2154 Based on Insight-HXMT

On 2020 April 27,the soft gamma-ray repeater SGR J1935+2154 entered its intense outburst episode again.Insight-HXMT carried out about one month observation of the source.A total number of 75 bursts were detected during this activity episode by Insight-HXMT,and persistent emission data were also accumulated.We report on the spin period search result and the phase distribution of burst start times and burst photon arrival times of the Insight-HXMT high energy detectors and Fermi/Gamma-ray Burst Monitor(GBM).We find that the distribution of burst start times is uniform within its spin phase for both Insight-HXMT and Fermi/GBM observations,whereas the phase distribution of burst photons is related to the type of a burst’s energy spectrum.The bursts with the same spectrum have different distribution characteristics in the initial and decay episodes for the activity of magnetar SGR J1935+2154.

Statistical Properties of X-Ray Bursts from SGR J1935+2154 Detected by Insight-HXMT

As one class of the most important objects in the universe,magnetars can produce a lot of different frequency bursts including X-ray bursts.In Cai et al.,75 X-ray bursts produced by magnetar SGR J1935+2154 during an active period in 2020 are published,including the duration and net photon counts of each burst,and waiting time based on the trigger time difference.In this paper,we utilize the power-law model,dN(x)/dx∝(x+x_0)~((-α)_x),to fit the cumulative distributions of these parameters.It can be found that all the cumulative distributions can be well fitted,which can be interpreted by a self-organizing criticality theory.Furthermore,we check whether this phenomenon still exists in different energy bands and find that there is no obvious evolution.These findings further confirm that the X-ray bursts from magnetars are likely to be generated by some self-organizing critical process,which can be explained by a possible magnetic reconnection scenario in magnetars.

Potential for Constraining Propagation Parameters of Galactic Cosmic Rays with the High Energy Cosmic-radiation Detection Facility on Board China's Space Station

Precise measurements of the spectra of secondary and primary cosmic rays are crucial for understanding the origin and propagation of those energetic particles.The High Energy Cosmic-radiation Detection(HERD)facility on board China's Space Station,which is expected to operate in 2027,will push the direct and precise measurements of cosmic-ray fluxes up to PeV energies.In this work,we investigate the potential of HERD for studying the propagation of cosmic rays using measurements of boron,carbon,and oxygen spectra.We find that,compared with the current results,the new HERD measurements can improve the accuracy of the propagation parameters by 8%–40%.The constraints on the injection spectra at high energies will also be improved.

The Year-scale X-Ray Variations in the Core of M87

The analysis of light variation of M87 can help us understand the disk evolution.In the past decade,M87 has experienced several short-term light variabilities related to flares.We also find that there are year-scale X-ray variations in the core of M87.Their light variability properties are similar to clumpy-ADAF.By re-analyzing 56Chandra observations from 2007 to 2019,we distinguish the“non-flaring state”from“flaring state”in the light variability.After removing flaring state data,we identify four gas clumps in the nucleus and all of them can be well fitted by the clumpy-ADAF model.The average mass accretion rate is~0.16M⊙yr^(-1).We analyze the photon index(Γ)-flux(2-10 keV)correlation between the non-flaring state and flaring state.For the non-flaring states,the flux is inversely proportional to the photon index.For the flaring states,we find no obvious correlation between the two parameters.In addition,we find that the flare always occurs at a high mass accretion rate,and after the luminosity of the flare reaches the peak,it will be accompanied by a sudden decrease in luminosity.Our results can be explained as that the energy released by magnetic reconnection destroys the structure of the accretion disk,thus the luminosity decreases rapidly and returns to normal levels thereafter.