Low Energy Threshold Research Articles

Layout optimization and Performance of Large Array of imaging atmospheric Cherenkov Telescope (LACT)

Abstract Large Array of imaging atmospheric Cherenkov Telescope (LACT) is an array of 32 Cherenkov telescopes with 6-meter diameter mirrors to be constructed at the LHAASO site. In this work, we present a study on the layout optimization and performance analysis of LACT. We investigate two observation modes: large zenith angle observations for ultra-high energy events and small zenith angle observations for lower energy thresholds. For large zenith angles (60°), simulations show that an 8-telescope subarray can achieve an effective area of $3 ~\rm km^2$ and excellent angular resolution. For small zenith angles, we optimize the layout of 4-telescope cells and the full 32-telescope array. The threshold of the full array is about $200~\rm GeV$, which is particularly crucial for studying transient phenomena, including gamma-ray bursts (GRBs) and active galactic nuclei (AGNs). This study provides important guidance for the final LACT layout design and performance estimates under different observational conditions, demonstrating LACT's potential for deep observations of ultra-high energy \gray sources and morphological studies of PeVatrons, as well as time-domain \gray astronomy.

Reactor neutrino liquid xenon coherent elastic scattering experiment

Coherent elastic neutrino-nucleus scattering (CEνNS) provides a unique probe for neutrino properties beyond the Standard Model physics. The reactor neutrino liquid xenon coherent scattering experiment (RELICS), a proposed reactor neutrino program using liquid xenon time projection chamber technology, aims to investigate the CEνNS process of antineutrinos off xenon atomic nuclei. In this work, the design of the experiment is studied and optimized based on Monte Carlo simulations. To achieve a sufficiently low-energy threshold for CEνNS detection, an ionization-only analysis channel is adopted for RELICS. A high emission rate of delayed electrons after a big ionization signal is the major background, leading to an analysis threshold of 120 photoelectrons in the CEνNS search. The second largest background, nuclear recoils induced by cosmic-ray neutrons, is suppressed via a passive water shield. The physics potential of RELICS is explored with a 30 kg·yr exposure at a baseline of 25 m from a reactor core with a 3 GW thermal power. In an energy range of 120–300 photoelectrons, corresponding to an average nuclear recoil from 0.63 to 1.36 keV considering the liquid xenon response and detector-related effect, we expect 4639.7 CEνNS and 1687.8 background events. The sensitivity of RELICS to the weak mixing angle is investigated at a low momentum transfer. Our study shows that RELICS can further improve the constraints on the nonstandard neutrino interaction compared to the current best results. Published by the American Physical Society 2024

Programmable multi-wavelength distributed feedback laser array integrated in a liquid crystal polymer waveguide.

Liquid crystal (LC) distributed feedback (DFB) lasers hold significant potential for integrated photonics applications. However, limitations in wavelength spacing for wavelength switching, device size, and compatibility with other technologies have impeded advancements of the LC DFB laser in integration and responsiveness. Herein, we propose a thin-film multi-wavelength DFB laser array utilizing high-resolution patterned programmable nematic LC polymers, enabling rapid switching with high-resolution wavelength spacing between wavelength division multiplexing channels while maintaining a stable single longitudinal mode (SLM) for each laser. The underlying physical mechanism involves modulating the effective refractive index of the DFB laser by varying the LC molecules' orientation angles between adjacent regions of the LC grating to achieve wavelength modulation. Additionally, a specialized LC waveguide design connects the DFB lasers, facilitating wavelength modulation as well as straight-line and bending propagation of the laser. Furthermore, the laser array demonstrates a relatively low energy threshold, facilitating its applications in high-integration scenarios.

LANTERN: A multichannel light calibration system for cryogenic detectors

Recent advancements in low-energy rare-event searches rely on the use of cryogenic calorimeters in order to reach extremely low energy thresholds. When attempting to accurately characterize their responses within the region of interest (ROI) these detectors face severe challenges due to the fact that traditional radioactive calibration sources produce signals above the ROI and cannot be used in low background environments.LANTERN, an optical calibration system, can instead be used to accurately study the ROI of cryogenic detectors. It utilizes LED pulses to generate photon signals to perform detector calibration. LANTERN’s current LED matrix allows for a fast characterization of up to 64 calorimeters independently, spanning energies from a few eV to hundreds of keV. Its minimal electronics and optics make it a cost-effective and customizable solution and it will be employed by the BULLKID, NUCLEUS and CRAB experiments.

Comprehensive study of Kinetic Trajectory Simulation method for multi-component magnetized plasma-wall interaction process

For a wide range of plasma applications across diverse fields, a comprehensive understanding of the plasma-wall interaction mechanism is indispensable due to its inherent connection with confined plasma. This work compilation delves into the Kinetic Trajectory Simulation (KTS) method for the interaction of multi-component magnetized plasma with wall, specifically focusing on its implications for the tungsten wall sputtering model. In the evolution of the 1d3v (one dimensional spatial coordinates and three dimensional velocity coordinates) KTS method, the coupled set of kinetic equations has been solved under specified boundary conditions which yields results of higher accuracy. At the particle injection boundary, we have assumed the velocity distribution function of particle species to be cut-off Maxwellians, meeting essential requirements for plasma-wall transition processes: quasineutrality, sheath edge singularity, continuity of macroscopic fluid variables, and the kinetic Bohm sheath condition. The kinetic Bohm sheath condition, a fundamental criterion for plasma sheath formation, is extended for multi-component plasmas, accounting for the cut-off Maxwellian distribution of negatively charged particles. A comparative study of the kinetic Bohm sheath condition for cut-off and Boltzmann distributions reveals a deviation of less than 2.0% in magnitude. The concentration ratio of positive or negative ion species and the presheath side electron temperature influence various plasma-wall transition characteristics, including wall potential, Debye sheath thickness, particle densities, potential distribution, particle fluxes towards the surface, particle drift velocity, phase-space trajectory evolution, and physical sputtering of the tungsten surface. Although lighter ions possess higher energy when striking the surface, the physical sputtering yield of the tungsten surface is greater for heavier ions due to their lower threshold energy and larger collision cross-section. Furthermore, a comparative study of plasma-wall transition properties using kinetic and fluid approaches demonstrates qualitative similarities, with a notable deviation of approximately 4.0% in the magnitude in the vicinity of the material surface.

Production of Positrons by Cosmic Rays

It is shown that electromagnetic interactions of cosmic ray protons provide a noticeable contribution to positrons production. This is due to a combination of low energy threshold of electron–positron pair creation compared to the thresholds of pion creation in strong interactions and rapid decrease of the cosmic rays energy spectrum. Moreover, the electromagnetically produced positrons are very soft, therefore their annihilation with background electrons directly produces the observed 511 keV gamma-line.

Water Cherenkov muon veto for the COSINUS experiment: design and simulation optimization

COSINUS is a dark matter (DM) direct search experiment that uses sodium iodide (NaI) crystals as cryogenic calorimeters. Thanks to the low nuclear recoil energy threshold and event-by-event discrimination capability, COSINUS will address the long-standing DM claim made by the DAMA/LIBRA collaboration. The experiment is currently under construction at the Laboratori Nazionali del Gran Sasso, Italy, and employs a large cylindrical water tank as a passive shield to meet the required background rate. However, muon-induced neutrons can mimic a DM signal therefore requiring an active veto system, which is achieved by instrumenting the water tank with an array of photomultiplier tubes (PMTs). This study optimizes the number, arrangement, and trigger conditions of the PMTs as well as the size of an optically invisible region. The objective was to maximize the muon veto efficiency while minimizing the accidental trigger rate due to the ambient and instrumental background. The final configuration predicts a veto efficiency of 99.63 ± 0.16% and 44.4 ± 5.6% in the tagging of muon events and showers of secondary particles, respectively. The active veto will reduce the cosmogenic neutron background rate to 0.11 ± 0.02 cts·\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\cdot $$\\end{document}kg-1·\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$^{-1}\\cdot $$\\end{document}year-1,\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$^{-1},$$\\end{document} corresponding to less than one background event in the region of interest for the whole COSINUS-1π\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\pi $$\\end{document} exposure of 1000 kg·\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\cdot $$\\end{document}days.

Low-contrast detectability of photon-counting-detector CT at different scan modes and image types in comparison with energy-integrating-detector CT.

We aim to compare the low-contrast detectability of a clinical whole-body photon-counting-detector (PCD)-CT at different scan modes and image types with an energy-integrating-detector (EID)-CT. We used a channelized Hotelling observer (CHO) previously optimized for quality control purposes. An American College of Radiology CT accreditation phantom was scanned on both PCD-CT and EID-CT with 10 phantom positionings. For PCD-CT, images were generated using two scan modes, standard resolution (SR) and ultra-high-resolution (UHR); two image types, virtual monochromatic images at 70keV and low-energy threshold (T3D); both filtered-back-projection (FBP) and iterative reconstruction (IR) reconstruction methods; and three reconstruction kernels. For each positioning, three repeated scans were acquired for each scan mode, image type, and CTDIvol of 6, 12, and 24mGy. For EID-CT, images acquired from scans (10 positionings × 3 repeats × 3 doses) were reconstructed using the closest counterpart FBP and IR kernels. CHO was applied to calculate the index of detectability () on both scanners. With the smooth Br44 kernel, the of UHR was mostly comparable with that of the SR mode (difference: -11.4% to 8.3%, to 0.956), and the T3D images had a higher (difference: 0.7% to 25.6%) than 70keV images on PCD-CT. Compared with the EID-CT, UHR-T3D of PCD-CT had non-inferior (difference: -2.7% to 12.9%) with IR and non-superior (difference: 0.8% to 11.2%) with FBP using the Br44 kernel. PCD-CT produced higher than EID-CT by 61.8% to 247.1% with the sharper reconstruction kernels. The comparison between PCD-CT and EID-CT was significantly influenced by the reconstruction method and kernel. With a smooth kernel that is typically used in low-contrast detection tasks, the PCD-CT demonstrated low-contrast detectability that was comparable to EID-CT with IR and showed no superiority when using FBP. With the use of sharper kernels, the PCD-CT significantly outperformed EID-CT in low-contrast detectability.

Investigation of Electronic and Molecular Features of Zn3S3/PEG4000 Nano-Composite Using the DFT Method

Molecular geometry structures were accurately optimized to low convergence energy thresholds for the Zn3S3 cluster before and after adding Polyethylene Glycol (PEG4000). Density functional theory DFT/ B3LYP calculations with 6-113G (d, p) basis set were employed to investigate structural and electronic properties of Zn3S3/PEG4000 composite. The FTIR spectral lines were analyzed where an agreement of FTIR spectra of titled molecules was evaluated between experimental and theoretical findings of the active peaks of O–H, C–H, C=O, C–O–C, and Zn–S functional groups. The vibrational modes frequencies were systematically analyzed on the distribution basis of potential energy around the range 0–4000 cm-1 and observed 12 modes of vibrations for the Zn3S3 molecule, while 36 modes for the Zn3S3/PEG4000 compound. Frontier high occupied, and low unoccupied molecular orbitals (HOMO&LUMO) were calculated and plotted to obtain the energy gap (E𝒈) resulting from the difference between those orbitals. The promising indicator was obtained at increasing E𝒈 from (4.031 to 4.459) eV after adding PEG4000, pointing out the effect of polymer on the ZnS surface as a capping agent. Additionally, electronic features of the mentioned structures, such as IP, EA, Ef, E𝒈, 𝐶𝑝, χ, η, Ѕ, and ω, were calculated. Finally, the molecular electrostatic potential (MEP) diagram of Zn3S3 and Zn3S3/ PEG4000 and charge densities of isosurface and contour diagrams were estimated, showing the nucleophilic and electrophilic attack of these compounds.

DARWIN/XLZD: A future xenon observatory for dark matter and other rare interactions

The DARWIN/XLZD experiment is a next-generation dark matter detector with a multi-ten-ton liquid xenon time projection chamber at its core. Its principal goal will be to explore the experimentally accessible parameter space for Weakly Interacting Massive Particles (WIMPs) in a wide mass-range, until interactions of astrophysical neutrinos will become an irreducible background. The prompt scintillation light and the charge signals induced by particle interactions in the liquid xenon target will be observed by VUV-sensitive, ultra-low background photosensors. Besides its excellent sensitivity to WIMPs with masses above ∼5 GeV, such a detector with its large mass, low-energy threshold and ultra-low background level will also be sensitive to other rare interactions, and in particular also to bosonic dark matter candidates with masses at the keV-scale. We present the detector concept, discuss the main sources of backgrounds, the technological challenges and some of the ongoing detector design and R&D efforts, as well as the large-scale demonstrators. We end by discussing the sensitivity to particle dark matter interactions.

Searching for WIMPs with TREX-DM: achievements and challenges

The TREX-DM detector, a low background chamber with microbulk Micromegas readout, was commissioned in the Canfranc Underground Laboratory (LSC) in 2018. Since then, data taking campaigns have been carried out with argon and neon mixtures, at different pressures from 1 to 4 bar. By achieving a low energy threshold of 1 keV ee and a background level of 80 counts keV-1 kg-1 day-1 in the region from 1 to 7 keV ee , the experiment demonstrates its potential to search for low-mass WIMPs. Two of the most important challenges currently faced are the reduction of both, background level and energy threshold. With respect to the energy threshold, recently a new readout plane is being developed, based on the combination of Micromegas and GEM technologies, aiming to have a pre-amplification stage that would permit very low energy thresholds, close to the single-electron ionization energy. With respect to the background reduction, apart from studies to identify and minimize contamination population, a high sensitivity alpha detector is being developed in order to allow a proper material selection for the TREX-DM detector components. Both challenges, together with the optimization of the gas mixture used as target for the WIMP detection, will take TREX-DM to explore regions of WIMP's mass below 1 GeV c -2.

High-z gamma-ray burst detection by SVOM/ECLAIRs: Impact of instrumental biases on the bursts’ measured properties

Context. Gamma-ray bursts (GRBs) can be detected at cosmological distances, and therefore can be used to study the contents and phases of the early Universe. The 4−150 keV wide-field trigger camera ECLAIRs on board the Space-based multi-band Variable Object Monitor (SVOM) mission, dedicated to studying the high-energy transient sky in synergy with multi-messenger follow-up instruments, has been adapted to detect high-z GRBs. Aims. Investigating the detection capabilities of ECLAIRs for high-redshift GRBs and estimating the impacts of instrumental biases in reconstructing some of the source measured properties, focusing on GRB duration biases as a function of redshift. Methods. We simulated realistic detection scenarios for a sample of 162 already observed GRBs with known redshift values as they would have been seen by ECLAIRs. We simulated them at redshift values equal to and higher than their measured value. Then we assessed whether they would be detected with a trigger algorithm resembling that on board ECLAIRs, and derived quantities, such as T90, for those that would have been detected. Results. We find that ECLAIRs would be capable of detecting GRBs up to very high redshift values (e.g. 20 GRBs in our sample are detectable within more than 0.4 of the ECLAIRs field of view for zsim > 12). The ECLAIRs low-energy threshold of 4 keV, contributes to this great detection capability, as it may enhance it at high redshift (z > 10) by over 10% compared with a 15 keV low-energy threshold. We also show that the detection of GRBs at high-z values may imprint tip-of-the-iceberg biases on the GRB duration measurements, which can affect the reconstruction of other source properties.

Comparative study of the multiple wavelength lasing of nitrogen ions: the role of vibrational level-dependent photoionization

We report on an optical amplification and energy threshold of the two most prominent emission lines, 391.4 and 427.8 nm, of the cavity-less lasing of nitrogen ions pumped by femtosecond laser pulses. It was found that the two transitions both show optical amplification under a low gas pressure condition, while the 391.4 nm emission is barely amplified under high gas pressure. Moreover, the 427.8 nm emission presents a significant lower pump laser energy threshold and a larger gain factor than the 391.4 nm emission. Numerical simulations based on a three-state coupling model suggest that the smaller ionization Franck–Condon factor from the ground state of N2 to the vibrational level ν = 1 in X2Σ g + state of N2+ favors the formation of population inversion corresponding to the 427.8 nm emission. Meanwhile, the competition between the strong field ionization and excitation induced by the pumping laser requires higher laser intensity to acquire the population inversion for the 391.4 nm radiation, leading to a corresponding larger energy threshold.

Optimization of a BEGe Detector Setup for Testing Quantum Foundations in the Underground LNGS Laboratory

In this work, we report on tests performed with an experimental apparatus prototype based on a broad-energy germanium detector aimed at investigating topical, foundational issues in quantum mechanics: i.e., possible violations of the spin-statistics connection and models of dynamical wave function collapse. Our recent phenomenological analyses demonstrated the importance of pushing the research of new physics signal, predicted in these fields, to an energy range below 10 keV. We describe the development of the dedicated data acquisition system and of the pulse shape discrimination algorithm, which have already allowed us to get a factor two improvement in the lower energy threshold. Future plans are discussed to further improve the lower energy threshold to the level of a few keV.

Construction and initial performance evaluation of a newly developed plant PET system

In plant research, positron emission tomography (PET) is occasionally employed for physiological studies, offering valuable insights. However, the generally high cost of PET systems and their suboptimal design for plant research pose challenges to their application in this field. To address these issues, we have developed a new PET system optimized specifically for plant research. The PET detector ring was positioned vertically to enable measurements of plants in their normal upright position. The developed plant PET system features a transaxial field of view (FOV) of approximately 12 cm and an axial FOV slightly larger than 9 cm, allowing for the imaging of relatively small-sized plants. To facilitate imaging of taller plants, the PET system can expand the axial FOV by changing the subject height using a lab jack, enabling the imaging of taller plant species. The measured spatial resolution at the central FOV was 3.3 mm FWHM, and the sensitivity was 3.7%. The timing resolution was 6.78 ns FWHM with a lower energy threshold set to 350 keV. Phantom images simulating plants were successfully measured using the developed plant PET system. We conclude that the developed plant PET system holds promise for effective plant imaging.

Fabrication of double-cladding fluoroindate glass fibers and watt-level 2.7 μm laser application

We report a double-cladding fluoroindate fiber with three layers of step refractive index glass. The gain double-cladding fluoroindate fiber is fabricated by the traditional rob-in-tube method. It possesses a tunable refractive index from 1.46 to 1.51 via modification of glass components. The good thermal–mechanical performance of fluoroindate glass protects fiber against thermal damage. The luminescence characteristics and radiation parameters are studied, and the results indicate enhanced emission with the higher Er3+ doping concentration. The maximum output power and the slope efficiency of heavily erbium-doped fiber laser are about 1.2 W and 16.3 %, respectively. The numerical modeling of a heavily-erbium-doped fluoroindate fiber laser is also carried out and compared with our experimental results. The low phonon energy and high laser damage threshold of fluoroindate glass allow for higher output power and efficiency of fiber lasers. In brief, this is the first study on watt-level 2.7 μm fiber laser in all-solid double-cladding fluoroindate fibers to our knowledge.

The spectrum-averaged cross section investigation of 117Sn(n,n')117mSn and 67Zn(n,p)67Cu reactions

The spectrum averaged cross sections (SACS) in standard neutron field, e.g. 252Cf(s.f.), is a preferable tool for cross section evaluation and validation. A set of reaction measurements with high energy thresholds was previously performed. The presented work focuses on lower energy threshold reactions, namely on the inelastic scattering of the tin foil, more specifically the reaction 117Sn(n,n')117mSn, and the zinc foil reaction, namely 67Zn(n,p)67Cu. These reactions are of special interest due to their intermediate energy range, which is essential in classical reactor dosimetry and fast reactor dosimetry. The experiments were carried out in a standard neutron field formed by 252Cf(s.f.) source in Řež. The experimental results were compared with calculations using MCNP6.2, ENDF/B-VII.1 transport library, and ENDF/B-VIII.0 and IRDFF-II cross section data library. Additionally, the calculations using CEA code DARWIN/PEPIN2 using JEFF-3.0/A were executed. The obtained experimental SACS of previously measured reactions were in good agreement with the SACS calculated using the IRDFF-II library. Additionally, the calculational reaction rate of 67Zn(n,p)67Cu was in accordance with the experimental data in case of ENDF/B-VIII.0 nuclear data library. Moreover, the calculational results of 117Sn(n,n')117mSn obtained by DARWIN/PEPIN2 code (using JEFF-3.0/A nuclear data library) are closest to the experimental results.

Exploring the Potential of Residual Impurities in Germanium Detectors for MeV-Scale Dark Matter Detection

The direct detection of MeV-scale dark matter (DM) particles hinges on achieving an exceptionally low-energy detection threshold. Germanium (Ge) detectors, meticulously tailored with precise impurity compositions, hold the potential to enhance sensitivity to energy levels below the sub-electronvolt (sub-eV) range. This study explores the behavior of residual impurities inherent to Ge detectors at helium temperatures, unveiling a captivating freeze-out phenomenon leading to the formation of excited localized states known as dipole states. Using compelling evidence from relative capacitance measurements obtained from two detectors, we elucidate the transition of impurity atoms from free charge states to these dipole states as the temperature drops from 11 to 6.5 K. Our investigation comprehensively covers the intricate formation of these dipole states in both n-type and p-type impurities. Furthermore, we shed light on the electric field generated by these dipole states, revealing their ability to trap charges and facilitate the creation of cluster dipole states. Confirming findings from previous measurements, we establish that these excited dipole states exhibit a binding energy of less than 10 meV, offering an exceptionally low detection threshold for MeV-scale DM. Building upon this concept, we propose the development of a 1-kg Ge detector with internal charge amplification—an innovative approach poised to surpass electrical noise and enable the detection of MeV-scale DM with unprecedented sensitivity.

The Case of the Missing Very High-energy Gamma-Ray Bursts: A Retrospective Study of Swift Gamma-Ray Bursts with Imaging Atmospheric Cherenkov Telescopes

Gamma-ray bursts (GRBs) are particle acceleration sites that can emit photons in the very high-energy (VHE) domain through nonthermal processes. From 2004 until 2018, the current generation of Imaging Atmospheric Cherenkov Telescopes (IACTs) did not detect any GRB in the VHE domain. However, from 2018 to 2020, five detections were reported. In this work, we try to solve the case of the missing VHE GBRs prior to 2018. We aim to identify GRBs that might have eluded VHE detection in past years by the H.E.S.S., MAGIC, and VERITAS IACTs. To do so, we study GRBs with known redshift detected by Swift from 2004 until June 2022. We first identify all GRBs that could have been observed by these IACTs since 2004, considering observation conditions and visibility constraints. We assume a relation between the X-rays and the VHE gamma rays based on the VHE GRBs detected to date and combine this with the redshift measurements, instrument response information, and observation conditions to predict the observed VHE gamma-ray flux from the Swift X-Ray Telescope measurements. We report findings on 12 bright low-redshift GRBs that we identify as most likely to have been detected in the VHE domain by current IACTs. The rate of IACT-detectable GRBs with ideal observation conditions is <1 VHE GRB per year with the current configuration. With its lower energy threshold and higher sensitivity, this rate increases to ∼4 VHE GRBs per year with the Cherenkov Telescope Array.

Gamma/Hadron Separation Method for the HADAR Experiment

Ground-based arrays of imaging atmospheric Cherenkov telescopes (IACTs) are the most sensitive γ-ray detectors for energies of approximately 100 GeV and above. One such IACT is the High Altitude Detection of Astronomical Radiation (HADAR) experiment, which uses a large aperture refractive water lens system to capture atmospheric Cherenkov photons (i.e., the imaging atmospheric Cherenkov technique). The telescope array has a low threshold energy and large field of view, and can continuously scan the area of the sky being observed, which is conducive to monitoring and promptly responding to transient phenomena. The process of γ-hadron separation is essential in very-high-energy (>30 GeV) γ-ray astronomy and is a key factor for the successful utilization of IACTs. In this study, Monte Carlo simulations were carried out to model the response of cosmic rays within the HADAR detectors. By analyzing the Hillas parameters and the distance between the event core and the telescope, the distinction between air showers initiated by γ-rays and those initiated by cosmic rays was determined. Additionally, a Quality Factor was introduced to assess the telescope’s ability to suppress the background and to provide a more effective characterization of its performance.