Xing-Chi Han · Shuo Wang · Hong-Yi Wu · Zhi-Huan Li · Guo-Zhu Shi · Ruo-Fu Chen ·Yue-Zhao Zhang · Shou-Yu Wang

Abstract A particle detector array designed for light-charged particles, known as the CsI-bowl, was built for exit channel selection for in-beam γ-ray spectroscopy experiments.This device is composed of 64 CsI(Tl) detectors, organized in a structure reminiscent of a tea-bowl.High quantum efficiency photodiodes, characterized by their minimal mass, were employed to collect scintillation light.Its design, construction, particle identification resolution, and its effectiveness in relation to exit channel selection are described in this paper.In source tests, the optimal figure of merit for the identification of α-particles and γ-rays using the charge comparison method was found to be 3.3 and 12.1 for CsI detectors coupled to photodiodes and avalanche photodiodes, respectively.The CsI-bowl demonstrated effectiveness in identifying particles, specifically the emission of protons and α-particles in the 58Ni(19 F, xpyn) fusion–evaporation reaction, thereby enabling the selection of the desired exit channels.

Keywords Light charged particle detector array · Particle identification · Charge comparison method

1 Introduction

For many decades, high-resolutionγ-ray spectroscopy,facilitated by large HPGe detector arrays, such as GAMMASPHERE [1] and EUROBALL [2], has been instrumental in accumulating extensive knowledge about nuclear behavior at high spins and excitations.This method has served as a fundamental element in the study of nuclear structures [3].The fusion–evaporation reaction induced by heavy-ions is the most effective method to maintain a balance between high angular momentum and a large cross section.However, this reaction populates a considerable number of exit channels, which complicates the selection of a specific channel for investigation.As a typical solution,data analysis often involves the use ofγ-ray multiplicity gating techniques.

However, for nuclei of medium or low atomic number (A),the large cross sections for charged particle emissions result in numerous channels with comparableγ-ray multiplicities.These channels cannot be separated by employingγ-ray multiplicity gating alone.In these instances, conducting measurements for the simultaneous emission of protons,α-particles, andγ-rays in the fusion–evaporation reaction can enhance the selection of the channel of interest.Hence, a light-charged particle detector array is commonly used as an auxiliary detector to enhance the peak-to-background ratio inγ-ray spectroscopy, as observed in examples like Microball[4] and DIAMANT [5].These two ancillary detectors were composed of CsI(Tl) scintillators, which exhibit several advantages over silicon detectors, including low cost, ease of machining, and strong resistance to radiation.

The advent of high-speed analog-to-digital converters(ADCs) with excellent electronic resolution has catalyzed novel opportunities for digital processing of pulses from detectors, simplifying the configuration of readout electronics [6].Leveraging these benefits, digital data acquisition systems have become increasingly prevalent in recent years, surpassing traditional analog data acquisition systems in the realm of nuclear physics research [7, 8].

When comparing pulse shape discrimination (PSD) [9]with the traditional analog electronics used in microballs and DIAMANT, it was observed that the capability of particle identification (PID) with a digital data acquisition system can be further improved via pulse shape analysis (PSA) of the stored waveforms [10].

Recently, research interest has significantly increased in the study of nuclei withA≈ 80, owing to the presence of numerous intriguing nuclear structure phenomena, such as multiple chiral doublet bands with octupole correlations.For a detailed study of this mass region, a light-charged particle detector array, termed as the CsI-bowl, with a digital data acquisition system was constructed at Shandong University(Weihai) forγ-ray spectroscopy experiments.The following section outlines the technical layout of the proposed array.Section 2 describes the design and construction of the CsIbowl.In Sect.3, the PID resolutions of CsI-PD and CsIAPD detectors in the source tests are presented using charge comparison and waveform fitting methods.The results of the channel-selection performance in the58Ni(19F,xpyn) reaction with the AFRODITE array at iThemba LABS are shown.

2 Detector configuration

To select the exit channels of interest in the fusion–evaporation reaction, the light-charged particle detector array should satisfy the following criteria: (i)good PID resolution for light-charged particles, especially for protons andα-particles; (ii) large solid-angle coverage;(iii) minimum absorption and scattering forγ-rays; (iv)reasonable granularity to reduce the probability of multiple hits in one detector; (v) good resistance to radiation damage;and (vi) small size

For the detection of charged particles, scintillators,such as NaI(Tl) [11–13], CsI(Tl) [9, 14, 15], and LaBr3[16–18], are undoubtedly some of the most frequently and widely used in nuclear and particle physics experiments.The overall decay time of the emitted light pulse in certain scintillators will vary with the type of exciting radiation and such scintillators are thus capable of pulse shape discrimination.The CsI(Tl) scintillator is known to exhibit a predominantly two-component scintillation process, with a fast component having decay times ofτf= 0.4 ～ 1.0 μs ,and a slow component having decay times ofτs= 7.0 μs[4].The relative populations of the two components depend on the type of ionizing radiation [19, 20].In general, the fast decay timeτfincreases as the ionization density of the detected particle decreases, whereas the ratio between the fast and slow components decreases.Additionally, the high conversion efficiency of the CsI(Tl) scintillator (64,000 photons/MeV) results in signals with high signal-to-noise ratios, which can lead to good PID resolution [21].Thus, the CsI(Tl) scintillator has been chosen as the detector material.

Considering the geometric structure, the assembly of the AFRODITE array, absorption of theγ-rays, and size of the CsI(Tl) crystals were chosen as 10 mm ×10 mm ×5 mm,and the maximum energy deposited in those crystals was 35 MeV for protons and 130 MeV forα-particles.This is sufficient for depositing the emission-charged particles in the fusion–evaporation reaction.The crystals were produced at the Institute of Modern Physics at the Chinese Academy of Sciences (IMPCAS) [22].The scintillation light was collected using a Hamamatsu S3590 PIN photodiode (PD)or S8664-1010 avalanche photodiodes(APD) with photosensitive areas of 10×10 mm2.The crystals were coupled to the PDs or APDs using light-guided silicone grease (EJ-550), and special care was taken to ensure that the silicone grease was free of bubbles, typically by pressing and sliding to squeeze out the bubbles.Each crystal was placed in a SiO2container, which served as a reflector for the scintillation light and provided a compact CsI-PD or CsI-APD assembly.The face of the crystal was covered with aluminized Mylar foil of 2.0 μ m to improve light collection.All CsI(Tl) detectors were soldered onto a PCB substrate via two brass pins of PDs or APDs, and the signals were sent to a compact 16-channel integrated charge-sensitive preamplifier MODEL SPA02-16 [23].The preamplifiers were designed and manufactured by the China Institute of Atomic Energy (CIAE).

3 Performance of the detector array

The capability of PID for each CsI(Tl) detector was determined byγ-source60Co andα-source (241Am and239Pu mixed) a schematic view of the experimental setup is shown in Fig.1.The channel-selection performance of the CsI-bowl was examined in the58Ni(19F,xpyn) reaction with the AFRODITE array at iThemba LABS.During the experiments, the waveforms were recorded and saved on a computer using XIA Pixle-16 modules from XIA LLC [24].Offline digital signal processing was performed based on the ROOT data analysis framework [25].

Fig.1 (Color online) Schematic view of the experimental set-up using γ-source and α-source

3.1 PID resolution of the detector

Figure 2 shows four types of sample waveforms after amplitude normalization and baseline subtraction: (1) photon conversion directly in the PD, (2)α-particles impinging on the CsI(Tl)crystal, (3) protons impinging on the CsI(Tl) crystal, and 4)γrays impinging on the CsI(Tl) crystal.The rise time (from 10%to 90% peak amplitude) of these waveforms is approximately 0.08, 1.2, 2.3, and 2.9 μs , respectively.According to the different rise times of the waveforms, the PID resolution was examined using the charge comparison [26, 27] and waveform fitting methods [28, 29].The charge comparison method is based on an analysis of the relative population of fast and slow decay components to distinguish different types of particles by integrating a shorter window near the peak (Qshort) and longer window encompassing the majority of the waveform (Qlong).In the source tests, different integral intervals were employed for the analysis of the waveforms to determine the optimal parameters for the in-beam experiments.The PID resolution of a detector is usually quantified using the figure of merit (FOM)[30, 31] defined as

Fig.2 (Color online) Sample waveforms show different pulse shapes due to protons, α-particles and γ-rays impinging on the CsI(Tl) crystal, and the fast rise-time trace s(labeled 1) is from γ-ray interaction directly on the silicon photodiode.To facilitate the comparison, the waveforms were normalized to the same height at t=10 μs and their baselines were shift to zero

whereμand FWHM denote the centroids and full-width-athalf-maximum of theα- andγ-ray peaks, respectively.Give that the FOM values are not sensitive to theQlongvalue [30],only a short gate interval is adjusted to optimize the PID resolution.Figure 3(a) shows the FOM values with different short-gate integrals for the CsI-PD detector; the optimal value of FOM is obtained as 3.3 with a short gate interval of approximately 5 μ s.For the comparison of PID resolution,a PID parameter is usually set as the ratio ofQlongtoQshort,and the PID spectrum atFOM=3.3 is shown in Fig.3(b) to distinguishα-particle andγ-ray.Owing to the internal gain(50–100) of APD, the PID resolution of CsI-APD detector is significantly improved with a high signal-to-noise ratio,and the optimal value of FOM obtained is 12.1 with a shortgate interval of approximately 4.0 μ s, as shown in Fig.3(c)and (d).

In the waveform fitting method, the light output of the CsI(TI) scintillator can be modeled using a combination of two exponential functions with fast and slow decay time constants[32]:

Fig.3 For α-particles and γ-rays, the FOM values change as a function of short gate interval for a CsI-PD detector and c CsIAPD detector, and the PID spectra were obtained using the charge comparison method for b CsI-PD detector and d CsI-APD detector

whereL(t) denotes the light pulse amplitude at timet;τfandτsdenote the decay time constants for the fast and slow components, respectively, with amplitudesNfandNs.The preamplifier output signalU(t) is described by the convolution between the light output function,L(t) and preamplifier impulse response function,G(t), an exponential function with decay time constantτp.The functional form of the signalU(t) is obtained via Eq.4 for the times before and after the beginning of the signal att0

where constantCdenotes the baseline of the signal.The relative intensities of the fast and slow components depend on the stopping power dE∕dxof the particle in the detector.A higher ionization loss produces a higher density of free electrons and holes, leading to a faster component.

To accelerate the fitting procedure and improve the fitting accuracy, it is necessary to reduce the number of free parameters and provide suitable initial values.The baseline constantC, which is determined from the average value of the samples proceeding to the leading edge, is set to a fixed value for each waveform.In the fitting processes,the parameterτpis set to 15.4 μs.This value is determined by performing a large number of fittings of the pulse tails(ranging from 20 to 30 μs ) using a single exponential decay functionG(t).The fast decay constantτfwas limited from 0.4 to 1.0 μs , and the slow decay constantτswas fixed at 7.0 μs ,as mentioned above.The other four parameters (t0,τf,NfandNs) are determined by applying a least squares fit to the waveform.A linear fit to the beginning of the leading edge is implemented to obtain the initial valuet0.The intersection between the fitting function and baseline can be calculated to obtain the valuet0, where the signal begins to rise.The initial valueτf=0.7 μs was set as the central value for the limited range of 0.4–1.0 μs.Meanwhile,NfandNsare initially set to half the amplitude of the waveform.The fitting results forα-particle andγ-ray waveforms obtained using the analytical function Eq.3 are shown in Fig.4(a) and(b).The fast and slow components can be clearly obtained,and the proportions of the fast and slow components are obviously different forα-particle andγ-ray.Thus, the PID for different types of incident particles can be obtained from the ratios of the fast and slow components of the waveforms.Figure 5 shows the optimal two-dimensional PID spectra for the CsI-PD and CsI-APD detectors using charge comparison and waveform fitting methods.Theα-particles andγ-rays were clearly separated by both methods, whereas the PID capability could be slightly improved by the waveform fitting method when compared to that by the charge comparison method.However, the fitting method requires the storage of the entire waveform, which results in high memory consumption during the experiments.Therefore, the charge comparison method was employed to obtain PID information in the in-beam experiment, and the integral values ofQshortandQlongwere recorded.Although the PID capability of CsI-APD detector is better than that of CsI-PD detector, the CsI-PD detector is applicable for the in-beam experiments due to the disadvantage of the gain instability with respect to the temperature of APD [33].

Fig.4 (Color online) Experimental baseline subtraction waveforms(black dash) and the least square fit results (red line) of a α-particle and b γ-ray.The blue and magenta lines illustrate the fast and slow components of the signal, respectively

Fig.5 (Color online) Particle identification plots for α-particles and γ-rays in the source tests using the charge comparison method and waveform fit method, (a) and (b) for the CsI(Tl)-PD detector and (c) and (d) for the CsI(Tl)-APD detector.All units in the figure correspond to arbitrary units

3.2 In-beam performance of the CsI-bowl

The CsI-bowl was paired with the AFRODITE array [34] at iThemba LABS as a channel selection device during the58Ni(19F,xpyn) in-beamγ-spectroscopy experiment.This setup can cover the angular range from 35?to 85?downstream,and from 143?to 155?upstream, as depicted in Fig.6.The19F ions with an energy of 62 MeV were used to bombard a target comprising 1.0 mg/cm2of highly enriched58Ni evaporated onto a 12.8 mg/cm2Au backing.Throughout the fusion–evaporation reaction, several nuclei in excited states were produced, such as (1α3p)70Ge, (1α2p)71As, (αp)72Se,(3p)74Se, (2p1n)74Br, and (1p1n)75Kr.The energy regions of the emission protons andα-particles were approximately 3–18 MeV and 7–24 MeV, respectively.Theγ-rays from the decay of the excited nuclear states were recorded by the eight Compton-suppressed Clover detectors, while charged particles were detected by the CsI-bowl.Particle andγray coincidence events were selected by applying a 300-ns time-window.

Fig.6 (Color online) CsI-bowl inside the reaction chamber coupled with the AFRODITE array at the iThemba LABS.The inset shows the schematic drawing of the CsI-bowl, and 58 Ni target was placed in the center of the chamber

Fig.7 (Color online) PID resolution for protons and α-particles of the CsI-bowl using the charge comparison method in the 58Ni(19 F, xpyn)fusion–evaporation reaction.a Plot of Qlong with respect to Qshort.b Plot of PID parameter with respect to Qlong value.The PID parameter is obtained by the ratio of Qlong and Qshort

The two-dimensional PID histogram obtained using the charge comparison method and PID resolution plot of the CsI-bowl are shown in Fig.7(a) and (b), and approximately 93% of the detected light-charged particles can be distinguished as protons andα-particles.The exit channels of interest can be selected by gating different numbers ofα-particles and protons.Theγ-γcoincidence spectrum of the58Ni(19F,xpyn) reaction is shown in Fig.8(a).The prominent peaks in the spectra at 635, 728, 868, 714, and 1000 keV stem from74Se and71As.based on calculations using the PACE4 [35, 36] computer code, these nuclei are the frequent product of evaporation of19F with58Ni.Theγ-ray spectrum gated to two protons is shown in Fig.8(b).The yrastγray transitions of74Se under 10+were strongly enhanced,and almost all contaminations from other exit channels were eliminated.Theγ-ray spectrum in coincidence with oneα-particle andxprotons, achieved with PID gates, is shown in Fig.8(c), and theγ-rays of71As,72Se, and70Ge are clearly enhanced with respect to the channels withoutα-particle emission.Furthermore, it is an effective way to distinguish similar energies ofγ-rays from different exit channels using the light-charged particle gating technique.As shown in Fig.8(a), the 862-keVγtransition from 2+to 0+of the past states in72Se [37] is mixed with the 868-keVγtransitions from 6+to 4+of the past states in74Se [38].The exit channel of72Se is populated with oneα-particle and one proton emission from the compound nucleus, whereas the exit channel of74Se has three proton emissions.The 862-keVγtransition from72Se can be completely eliminated by gating the two protons, as shown in Fig.8(b), and the 868 keVγtransition from74Se disappeared by gating oneα-particle andxprotons, as shown in Fig.8(c).Another advantage of exit-channel selection is the enhancement of weakγtransitions against the background.As shown in Fig.8(a), the 821- and 831-keVγtransitions in71As and72Se cannot be discriminated using only theγ-γcoincidence spectrum because of the strongγ-ray background from other channels.Afterα-xp-γco-incidence gating, the twoγtransitions can be clearly distinguished, as shown in Fig.8(c).Consequently, these results provide unambiguous verification of the capability of the CsI-bowl as a channelselection device.

Fig.8 Gamma-ray spectra from the reaction of 58Ni(19 F, xpyn)with beam energy 65 MeV.a Total γ-γ projection spectrum, b γ-ray spectrum for two protons emission channels, and c γ-ray spectrum for α-xp emission channels

4 Summary and outlook

The CsI-bowl is designed, constructed, and commissioned as an ancillary detector for in-beamγ-spectroscopy experiments to select the exit channels.The particle identification resolutions of the CsI coupled with photodiode and avalanche photodiode detectors were studied usingγ-source60Co andα-source (241Am and239Pu mixed) via the charge-comparison and waveform fitting methods.The optimal value of the figure of merit for distinguishingα-particles andγ-rays was found to be 3.3 using the charge comparison method for CsI coupled with photodiode detectors, and the optimal value of the figure of merit was equal to 12.1 for the CsI coupled with avalanche photodiode detectors due to the internal gain of avalanche photodiode.The particle identification resolution can be slightly improved via the fit waveform method when compared to that obtained via the charge comparison method.In the58Ni(19F,xpyn) experiment, approximately 93% of the protons andα-particles were clearly identified, and the capability of the CsI-bowl as a channel selection device was determined.The present measurements highlight the improved sensitivity ofγ-ray spectroscopy using the CsI-bowl coupled with an AFRODITE array.

With the development of the new scintillation materials LaBr3(Ce) and CeBr3[39] in the last decade, hybrid LaBr3(Ce)-HPGe arrays [40, 41] have become mainstream detection arrays in nuclear spectroscopy experiments for the measurement of excited nuclear state lifetimes [42–44],providing essential insights into the structure of nuclei.A hybrid CeBr3/LaBr3(Ce)-HPGe array is currently under construction for nuclear structure studies at Shandong

University (Weihai).The hybrid CeBr3/LaBr3(Ce)-HPGe array allows the installation of up to 10 HPGe and 16 CeBr3/LaBr3(Ce) detectors.An updated version of the CsI-ball,which has nearly 4πsolid angle coverage, as an ancillary detector for the hybrid CeBr3/LaBr3(Ce)-HPGe array is also under construction.The scintillation light is collected by a silicon photomultiplier [45, 46] to further improve the particle identification capabilities.With the high gain (105–107) of silicon photomultiplier, the waveform of signals can be recorded directly via the digital data acquisition system without the preamplifier.The length of the waveform is approximately 3 μ s, which allows the entire waveform to be recorded at a high counting rate during the experiment.Owing to the strong function of the gain and temperature of the silicon photomultiplier, a temperature-compensated bias voltage system has already been developed, and the bias voltage is precisely controlled using a thermosensor to maintain a constant gain during the experiment.

AcknowledgementsThe exceptional skill demonstrated by the stafffrom the Crystal Group at IMPCAS during the construction of the CsI-bowl, as well as the outstanding collaboration with the Nuclear Reaction Group from the China Institute of Atomic Energy during preamplifier testing, is greatly appreciated.The authors also extend their gratitude toward the technical staffand accelerator group at iThemba LABS for their consistent support throughout the experiment.

Author ContributionsAll authors contributed to the study conception and design.Material preparation, data collection and analysis were performed by Shuo Wang, Xing-Chi Han, Hong-Yi Wu, Zhi-Huan Li,and Shou-Yu Wang.The first draft of the manuscript was written by Xing-Chi Han and all authors commented on previous versions of the manuscript.All authors read and approved the final manuscript.

Data availabilityThe data that support the findings of this study are openly available in Science Data Bank at https:/www./doi.org/10.57760/sciencedb.10201 and https:// cstr.cn/ 31253.11.scien cedb.10201.

Declarations

Conflict of interestThe authors declare that they have no competing interests.