Perspectives of a single-anode cylindrical chamber operating in ionization mode and high gas pressure

As part of the R2D2 (Rare Decays with Radial Detector) R&D, the use of a gas detector with a spherical or cylindrical cathode, equipped with a single anode and operating at high pressure, was studied for the search of rare phenomena such as neutrinoless double-beta decay. The presented measurements were obtained with a cylindrical detector, covering gas pressures ranging from 1 to 10 bar in argon and 1 to 6 bar in xenon, using both a point-like source of $^{210}$Po (5.3 MeV $\alpha$ ) and a diffuse source of $^{222}$Rn (5.5 MeV $\alpha$). Analysis and interpretation of the data were developed using the anodic current waveform. Similar detection performances were achieved with both gases, and comparable energy resolutions were measured with both sources. As long as the purity of the gas was sufficient, no significant degradation of the measured energy was observed by increasing the pressure. At the highest operating pressure, an energy resolution better than 1.5% full-width at half-maximum (FWHM) was obtained for both gaseous media, although optimal noise conditions were not reached.


Introduction
The search for neutrinoless double beta decay (β β 0ν) is currently one of the main goals in particle physics.Several ton-scale projects have been internationally identified such as LEGEND [1], CUORE [2], or nEXO [3].They aim to reach an effective Majorana neutrino mass, m β β , of ∼10 meV, ruling out the inverted hierarchy region in the absence of any signals, employing at least one ton of material enriched in double beta decay isotopes.However, reaching the desired a e-mail: pascal.lautridou@subatech.in2p3.frb e-mail: anselmo.meregaglia@cern.chc Now at IFIC Valencia d Now at CBMN Bordeaux low-background level with a complex detector counting a large number of readout channels is not a trivial task.This weakness is also highlighted in other techniques, such as in time projection chambers (TPCs) utilising pressurized xenon gas to image the two expected beta traces [4].Finally, the cost of the detector is also an important parameter to be considered.
An alternative approach could be provided based on the concept of the Spherical Proportional Counter (SPC).This kind of radial TPC, first developed at CEA Saclay [5], was initially devoted to the physics of low-energy neutrinos, such as neutrino oscillations, neutrino coherent elastic scattering and supernova neutrinos detection [6,7].In the last two decades, many important developments have been made.The SPC's simple design, consisting of a large grounded sphere with a small central spherical anode as a unique readout channel, minimizes the detector material budget.Furthermore, the SPC exhibits both a low energy threshold and a high gain, allowing the observation of the signal created by the production of a single primary electron.This technology was chosen by the NEWS-G collaboration [8][9][10] for the search of light dark matter candidates in the mass range between 0.1 and 10 GeV.Combined with a percent-level energy resolution, this detector could be an appealing device in the search for β β 0ν decay.Monte-Carlo studies of a 136 Xe-filled detector at 40 bar suggest a sensitivity covering the inverted mass hierarchy region [11] with a ton scale detector.R2D2 R&D began in 2018, with the hope of proving that an energy resolution of 1% full width at half maximum (FWHM) was achievable at 2.458 MeV, corresponding to the transition energy (Q β β ) of the 136 Xe double beta decay.Preliminary studies reported the operation of the SPC in proportional mode with argon [12,13] and xenon [14].As of 2021, a cylindrical proportional counter (CPC), using the same detec-Fig.1 CPC working principle and electric field computed with COM-SOL [15].
tor principle as the SPC, but with a wire as a central anode, was developed [14].The results showed a preference for a cylindrical cathode, therefore making the CPC the baseline design for the R2D2 detector.The main advances arise from three factors: the grounded anode, which results in noise almost independent of the high voltage applied to the cathode; the 1/r dependence of the electric field, inducing a stronger field away from the anode and minimizing the impact of electronegative impurities; and, finally, greater homogeneity of the field compared to that caused by the rod supporting the anode.However, as discussed in Sec. 7, the main difference lies in choosing to operate the CPC in ionization mode rather than in proportional mode.This latter mode offers the significant advantage of reducing the bias demand (which was only available up to 6 kV), and eliminating the gain fluctuations inherent to the avalanche mechanism for energy resolution.
In this report we present the current state of this development.The experimental setup is discussed in detail in Sec. 2. Since gas purity is paramount to a correct detector operation, the gas recirculation system and the obtained purity are discussed in Sec. 3 and 4, respectively.The effect of temperature variation is addressed in Sec. 5. Signal modeling and processing are the subject of Sec. 6, and data analysis is presented in Sec. 7.

Experimental setup
The detector is a CPC, consisting of three basic elements: a supporting structure, a cylindrical cathode, and a central wire.Moreover, in order to operate the CPC at pressures higher than atmospheric, a certified vessel withstanding the pressure is required.The CPC working principle is identical to any gas proportional counter: ionisation tracks traversing the gaseous medium deposit energy and produce ionisation electrons.An electric field, produced by a negative high-voltage applied on the cylinder (cathode), drifts the electrons towards the grounded central wire (anode), which is used to read out the induced signal.Depending on the value of the electric field, the electrons could undergo an avalanche process near the wire (proportional mode) or simply be collected with no amplification (ionisation mode).The electric field, computed with COMSOL software [15], as well as the detector working principle can be seen graphically in Fig. 1.

High pressure vessel
To avoid additional costs, the available spherical vessel (of 20 cm radius and volume of approximately 33 l), built in the first R2D2 phase and tested as SPC, was used as the high pressure vessel to host the CPC.The cathode dimensions (with height and radius of 27 cm and 11 cm, respectively) were therefore chosen in order to fit inside the existing vessel.To hold the counter vertically inside the sphere, two specifically designed pieces were attached to both end-caps of the supporting structure.A picture of the system (with an epoxy-copper cathode), set in the bottom hemisphere of the spherical vessel, is shown in Fig. 2.

CPC supporting structure
The skeleton consists of three cylindrical pillars holding two end-caps together.These end-cap pieces have the shape of a spoked wheel (to provide stability to the detector) surrounded by a groove (to hold the cathode in place).In order to ensure an optimal gas purity, as discussed in the following sections, all the materials were selected to minimize the outgassing.Several plastics polymers were considered, such as PVC (polyvinyl chloride), Teflon (PTFE polytetrafluoroethylene), PEEK (polyether ether ketone) and POM (polyoxymethylene), and their outgassing was measured.The procedure involved inserting the material in a dedicated vacuum chamber and pumping for 30 minutes.Afterward, the pump was turned off, and the subsequent pressure increase in the vacuum chamber was monitored over a few minutes to assess outgassing.Teflon was chosen as the optimal material for the CPC structure, as indicated in the results presented in Fig. 3.

CPC cathode
The cathode choice was also motivated by the lowest possible outgassing.As part of the CPC skelton, it must to be thin enough to be bent and inserted in the end-cap grooves, whilst assuring rigidity to the detector.A sheet of aluminium 200 µm thick, immediately available, was used as a cathode.The distance between the CPC wall and the grounded spherical vessel is large enough to safely apply high-voltage at the level of 5 kV without sparking.An insulating layer on the outer side of the cathode would be desirable, but this greatly increases the risk of outgassing.Two copper cathodes covered with kapton and epoxy as insulating materials were tested, revealing important outgassing: this is likely due to the material itself (epoxy) or to the process used to attach the two materials (glue between kapton and copper).The outgassing results of different components of the CPC and different cathodes are shown in Fig. 4: the full CPC with an aluminium cathode is almost identical to that obtained with the CPC structure only, proving the aluminium cathode is not a source of oxygen or other electronegative impurities.This cathode was therefore chosen for the presented measurements.The final setup was baked at 100 degrees Celsius for one week in a vacuum chamber to reduce the outgassing of the materials during operation.

CPC central wire
The diameter of the central wire is a principal parameter and is selected according to the chosen detector operation mode (i.e.ionization or proportional).In the case of ionization, a thicker wire ensures a larger electric field far from the wire, and the possibility to increase the electron's drift velocity, without starting an avalanche near the anode.In proportional mode, a thinner wire allows reaching the multiplication region with a relatively small high voltage (i.e.below 10 kV).The measurements presented were conducted with two different anodes: a tungsten wire with a 50 µm diameter, and a metallic rod with a 1.2 mm diameter.The wire/rod is connected to a feed-through on the top of the vessel and read through the electronics chain and DAQ to record the signal.The high voltage is applied to the CPC wall via a wire connected to a feed-through at the bottom of the vessel.The coupling between the signal and high-voltage is minimized (at the level of the irreducible anode-cathode capacitance) ensuring the weakest dependence of the high-voltage on the noise.This choice was a major improvement over previous prototypes [12,14], where the increase in high-voltage at high pressure subsequently increased the electronic noise, limiting the energy resolution.

Read-out electronics
The electronics chain consists of a low-noise charge preamplifier sending its waveform to a CALI card (ADC) which was controlled by the SAMBA acquisition software [16].The integrator output signal is continuously sampled at 2 MHz.An internal digital amplitude threshold discriminator commands the capture of the waveform of interest for a duration of 3 ms.The first 1.5 ms is devoted to estimating the baseline of the event, while the second half of the time window, containing the triggering physical transient, undergoes further digital processing (see Sec. 6).The event trigger was carried out by positioning a fixed digital threshold on the amplitude of the wave train.For this program of operation in ionization mode, the charge-sensitive low noise and high gain amplifier ORTEC 142PC was adopted to read the anode sensor.

Calibration source
A 210 Po source, emitting α's at 5.3 MeV, is used to assess the detector energy resolution.The polonium is deposited on a silver plate of 0.6 × 0.6 cm 2 , and positioned on the outside of the cathode behind a hole of 1 mm radius, through which the α particles enter the CPC active volume.The activity of the source was roughly 10 Bq, corresponding to a rate of 0.8 events per second in the detector.

Gas recirculation system
As with all detectors utilising significant drift distances of the primaries ionisation electrons, a clean gas is crucial to limit the effects of electronegative impurities, such as oxygen or water.These molecules can capture the ionization electrons, which then migrate towards the anode, thus reducing the charge collected.In the case of noble liquid detectors, several solutions exist to guarantee a purity at the level of tens of ppt (part per trillion) or, as more commonly measured, an "electron lifetime" of several milliseconds.Currently, the most effective purification method, used by the XENON collaboration, is based on cryogenic filters built from copper [17]: oxygen impurities react with copper (2Cu+O 2 → 2CuO) and become permanently bound to its surface.Other purification processes could be applied such as the spark discharge technique [18], or cold and hot getters [19][20][21].The use of commercial getters does not require cryogenic and distillation infrastructures, and are much cheaper compared to the cryogenic filter.However, the best achievable gas purity is then reduced to the ppb (part per billion) level.
An additional aspect which is normally taken into account in gas/liquid purification is the radon contamination [22].Radon is the dominant source of background for experiments searching for rare events and it is typically removed through distillation columns [23].However, in the present phase of the R2D2 R&D, there is no low background requirement and the presence of radon does not represent an issue.On the contrary, as will be discussed in the following, radon will be exploited as a diffuse source of calibration.
In order to operate the R2D2 prototype with sufficient gas purity conditions, a combination of cold and hot getter was chosen.A PS3-MT3-R-2 MonoTorr Entergris gas purifier was used as hot getter, granting outlet impurity levels for oxygen reduced to ppb levels.The hot getter was coupled to a Mott micro-bulk gas purifier (MGP-30-125-IG-106-1.5NM-V3-V) cold getter, located upstream, to provide a first purification stage and increase the lifetime of the hot Fig. 5 Recirculation system including cold and hot getter, as well as the flow meter and the recirculation pump chamber.getter filtering components.This decision was mainly motivated by the allocated R&D budget, which was not sufficient to have purification based on cryogenic filters.Nonetheless, in terms of electronegative impurities, the contamination at the level of few ppb is sufficient to operate the detector up to 10 bar.The 10 bar limit is determined by the maximal operating pressure of the getter.The gas recirculation system is the second component ensuring the detector operates under the required conditions.The gas was recirculated using a KNF pump (NPK03KVDC-B4) installed in a dedicated chamber to prevent leaks and withstand the gas pressure.The gas flow, measured by a Bronkhorst El-FLOW flowmeter, can be adjusted by varying the voltage applied to the pump.The two limiting factors are a maximal allowed flow of 10 l/m in the hot getter, and a limit on the rotation speed of the pump which induces vibrations on the setup increasing the noise.The latter limit depends on the gas pressure and was set, in practice, by requiring the same noise level when taking data with and without the recirculation system.A picture of the recirculation setup can be seen in Fig. 5.
To establish the need of a persistent recirculation while taking data, the gas purity degradation was assessed as a function of time.For this, a dedicated run was taken with the detector filled with xenon at 3 bar, exploiting the 210 Po calibration source on the cathode.The mean reconstructed charge corresponding to α events was monitored over the whole run.At the beginning of the run the detector was in a stable condition, with the gas recirculating at the speed of 2.2 l/m for a few days. 1 h after the run began, the recirculation was stopped, and a constant reduction in the collected charge of about 1.5% per hour was observed.Such a reduction demonstrates that a permanent recirculation (or at least regular) is needed during data taking to guarantee an optimal gas quality.16 h after the run start, the recirculation was restarted, and the same reconstructed charge achieved at the beginning of the run was retrieved after only ∼ 4 h.These observations are summarized in Fig. 6.

Gas purity
Gas purity effect is usually summarized in terms of "electron lifetime".This is a typical estimator of the contamination in terms of electronegative impurities in noble liquid/gas experiments.Indeed, in such TPCs, the number of charges present in the drift space (and finally the collected charge) for a given energy deposit, depends on the electron drift time (and therefore on the position of the primary ionization along the drift direction).For a deposited charge producing a total number of N(0) electrons at time t = 0, the collected number of electrons N(t) after a drift time t should follow the relation [17]: where τ is the so-called "electron lifetime".A long electron lifetime implies a smaller recombination during the drift, and thus a medium with a small contamination of electronegative impurities.It must be noted that the impact of the electron lifetime on the reconstructed signal is different for detectors operated in proportional and ionisation mode and, from the point of view of the gas purity, operating the detector in the ionization regime should provide a slight advantage.The ionization signal begins as soon as the primary electrons are in motion, whereas the proportional signal, formed by the Fig. 7 Ratio between attenuated and total signal for the induced current (orange line) and for the total collected charge (blue line).The green line shows the difference between the two signal ratios.
backwards drift of the produced ions, depends instead on the number of primaries which enter in the avalanche region, i.e. those that reach the anode.The effect of a loss of charges during migration is presented for the two regimes in Fig. 7.The calculation is carried out in the extreme case of a point-like interaction at the cathode position.The ionization signal ratio with and without charge losses (IC Att /IC) is compared to the ratio of charges (Q Att /Q) for different values of the electron lifetime.The curves show that ionization mode presents better immunity to this effect and that the constraints on the gas purity could be less stringent.Furthermore, in ionization, this effect could likely be tuned by adjusting the induced signal shape, thereby accounting for the loss of primaries during transport (see Sec. 6).
In R2D2, the 210 Po source is located on the cathode and it is not possible to vary its position radially along the drift direction.Another way to vary the drift time is to change the high-voltage applied on the cathode, and therefore the electric field, leaving the α source position unchanged.The two methods are not exactly equivalent, since the modification of the electric field to change the electron drift time could also modify the interactions between drifting electrons and electronegative impurities, and therefore the attachment.Considering such an effect, and that R2D2 is operated in ionisation mode, the exponential fit of Eq. 1 is only an approximation for the electron lifetime evaluation.Considering the method's limitations, measurements of the signal integral were taken at various high-voltages, i.e. different electron drift times.However, in the present experimental setup, there is no independent measurement of the drift time since the α emission time is unknown, and no trigger based on the scintillation light is present.As is well established for chambers working in ionization regime, an experimental measurement of the drift time (and therefore the position of the interaction) can instead be deduced from the estimate of the signal duration (see Sec. 6).For validation, the GARFIELD++ [24,25] simulation toolkit was used to evaluate the drift time of electrons from the cathode to the central anode at different pressures and high-voltages.GARFIELD++, developed for the detailed simulation of particle detectors based on ionisation measurement in gases, had already been benchmarked against the R2D2 data in a spherical geometry [12], and showed very good agreement.Data were taken at different pressures between 1 and 10 bar with small high-voltages, ensuring a high drift time and no avalanche near the anode.The values of the high-voltage used in this study, which also enhance the electron attachment phenomena, are typically much lower than those used in standard detector operation conditions described in Sec.Nonetheless, comparing different GARFIELD++ distributions, large discrepancies on the drift time were observed at small electric field.In particular, the distribution available in 2021 was benchmarked against the one released in 2023.A summary of the obtained results can be seen in Fig. 9.
The discrepancies between the different GARFIELD++ releases are important at pressures below ∼ 3 bar and almost disappear at pressures above 5 bar.
An electron lifetime below 1 ms is found in argon at 10 bar, and below 2 ms in xenon at 6 bar, showing that the gas purity does not match the quality to that already achieved by liquid noble gases experiments with electrons lifetimes greater than 20 ms in argon [26] and 10 ms in xenon [17].Present performances of the R2D2 setup could certainly be improved by means of a qualitative leap in gas purity.Two identified possibilities are: baking the full device at a higher temperature, and the replacement of the recirculation pump (an important source of electronegative impurities) with a cleaner one e.g. a magnetically driven piston pump [27].

Temperature effects
To precisely evaluate the energy resolution of the detector, the impact of some external factors on the detector stability were also studied.In particular, measurements conducted during the R&D phase revealed that the electronic response is sensitive to temperature changes.Furthermore, it is plausible that temperature fluctuations in the detector's surrounding environment directly influence the internal gas pressure, thereby contributing to the observed behavior.These effects are negligible in short runs exploiting the 210 Po calibration source, since the α rate of few Hz allows for a large statistics in less than one hour of data taking.However, temperature changes become relevant in long runs, needed to exploit 222 Rn in the detector, which last more than one day.Dedicated measurements taken while monitoring the detector temperature with a precision of 0.05 degrees Celsius, showed that the α signal integral variations are in phase with temperature changes.Temperature variations of about 1 degree could re-  sult in energy variations at the per-cent level, as can be seen in Fig. 10.
The operation of the detector in a temperature controlled environment is desirable, however not possible for the current experimental setup.Therefore it was decided to constantly monitor the detector temperature and, as explained in Sec. 7, apply corrections for runs longer than few hours if needed (i.e. if temperature changes larger than 0.5 degrees were registered).
6 Signal modelling and processing

Signal modelling
The principle of calculating the current induced by the movement of charges between the electrodes has been established for decades [28,29].This method is widely used in the field of particle detection simulation [30].In a parallel plate ionization chamber (PPC), examining the current signal enables the extraction of various physical information related to the interaction.This includes details such as the deposited energy and the distance of the interaction from the anode [31,32].These features led to the selected approach based on the exploitation of the current signal.Such a signal is extracted by digital processing of the original output waveform of the integrator [12,14,33].To interpret the results, a dedicated modeling of the expected current waveform was jointly developed.Compared to GARFIELD++, it is a simplified and fast simulation based on the average macroscopic behavior of the primaries during their transport in the gas.This choice was motivated by two observations: first the main parameters which shape the signal are more easily identifiable; and secondly, in the ionization regime (where the signals are very weak) the fluctuations observed are essentially dominated by Radio Frequency Interferences (RFIs) and electronic noises.Starting from the GARFIELD++ results giving the drift time of a primary as a function of the distance, the average veloc- ity of the charge v e (r) is first deduced.The expression of the electric field of a cylindrical counter is given by: where r is the radial distance from the anode, V 0 the potential difference between anode and cathode, E w the so-called weighted electric field, and r cathode /r anode the ratio of cathode and anode radii, respectively.By taking into account Eq. 2 and referring to the Shockley-Ramo theorem, for this simple geometry the current I(r) induced on the anode by an elementary charge drifting all along the radius of the CPC can be calculated via: where e is the electron charge and v e (r) the drift radial velocity at the r position.An example of the shape of this current is presented in Fig. 11 for the CPC filled with a Xe gas at 3 bar and a high-voltage of 3000 V.Note that in the developed simulation, the contribution of ions is omitted because their speed is several orders of magnitude lower (while in proportional mode, most of the signal is induced by the ions produced in the avalanche).Furthermore, the very long ion component would be well below the detector sensitivity, particularly in xenon.The introduction of the distribution of primary charges and the conversion of the distance into sampled time then allow the construction of the signal induced by a full track.The penultimate stage of the simulation chain introduces the attachment effect as well as the instrumental noise.Convolution of the resulting waveform with the electronics response completes the construction of the induced current signal, called S(t).While constant for a PPC, in case of the CPC, the S(t) signal varies considerably in time (see Fig. 12), usually increasing slowly at first.Consequently, the conventional observables used for PPCs are no longer suitable, and new physical observables must be constructed to correctly classify the registered events.This necessitates first fixing the analyzed relative fractions of the signal.Since S(t) is typically small at its onset, the ideally reintegrated waveform Sr(t) is employed for this purpose.Indeed, the latter makes it possible to immediately deduce the total induced charge Qt, which is given by the maximum of Sr (the integral of S).This allowed us, for temporal estimation, to precisely set the crossing points above the thresholds, defined in terms of a fraction of Qt.With both temporal representations, it becomes possible to build the three observables used in the analysis, namely Qt, Dt and Dh, which contain different information about the interaction in the gas.An example of a signal with the associated observables is presented in Fig. 12.
The left part of the signal is due to the induction in the weakly increasing electric field, far from the anode.In this region of electric field, the detector functions almost as a PPC.The strong growth on the right part of the signal is caused by the rapid increasingly variation of the electric field (as 1/r) in the region ∼ 1 cm from the anode [34].This is typical of a signal observed with a wire proportional chamber (WPC) [35].
After reaching a maximum, as the primaries are collected by the anode, the induction vanishes.The descending part of the transient becomes mostly governed by the impulse response of the amplifier (undershoot), rather than by the arrival tail of the electrons.Each observable describes a specific feature of the registered signal: -Qt is the maximum value of the integrated signal Sr, noting that because of the non constant electric field, it is not directly the true initial ionization charge Q.To reconstruct the true energy deposited (at first order proportional to the number of primary charges created), Qt must be corrected for the radial position of the interaction, especially when the interaction is close to the anode.However, by selecting a narrow range in energy deposit positions, Qt can provide a rough estimate of the deposited energy.
To complete the specificity of the ionization signal as a consequence of the Shockley-Ramo theorem, it should be also noted that Qt is not affected by diffusion during the transport of primaries.-Dt is the total duration of the signal S from the moment it rises above threshold (constant fraction of the total integral) up to the maximum of the signal Sr (i.e. when S goes to zero).As is the case in a PPC, Dt gives a direct measurement of the maximum radial distance from the track to the anode (see Figs. 13(a) and 13(b)).-Dh is the width of the signal S at half height.This observable has no specific meaning in a PPC, and is linked to the radial extent of the track demonstrated by the developed signal simulation for the CPC.This feature makes Dh particularly interesting in the event selection applied for the energy resolution analysis.From the Shockley-Ramo theorem it follows that for a point-like energy deposit, Dh is independent on the radial position inside the CPC.However, a track can be considered as the sum of energy deposits at the same time and different radial distances.Each one would be the seat of the same radial evolution, but shifted by the difference in drift time at the considered positions.Such a superposition results in a signal with a larger Dh for tracks extending in the radial direction (see Fig. 13(c)).
As an example of the behavior of these quantities, the current waveforms, including noise and electronic response, were modeled in 1 mm steps from the cathode to 1 cm from the anode.Predictions of the observables Qt and Dt as a function of the distance r to the anode for a point-like deposit are presented in Fig. 13(a).The dependencies demonstrate that an effective recognition of the endogenous decays (betas and alphas) emitted by the anode or the cathode can be obtained, making it possible to better reject this contribution from the radioactive background.In Fig. 13(b) the dependence of Dt as a function of Qt is shown, whereas a simulation of the evolution of Dh as a function of the track length is presented in Fig. 13(c).Ionization lines starting from the cathode and pointing towards the anode were selected.For inclined tracks, projecting the tracks onto the radial axis ultimately results in a similar distribution being followed.The computed response of the CPC in ionization mode combines advantages specific to PPC with some distinctive of WPC.The CPC signal is produced throughout the drift of the electrons towards the anode, offering original and attractive detection functionalities, in particular with a single read-out anode.

Signal processing
Regarding the experimental data, which relies on the treatment of the registered integrator waveform, several numerical operation are needed to retrieve the current signal.The first step consists of the application of a Fourier transform base fil- ter with the aim of mitigating the noise caused by the electronics and the detector vibrations.A deconvolution algorithm is then applied in order to retrieve the current signal.The registered waveform is indeed a convolution of the physical signal with the typical RC shaping of the preamplifier: such a component has to be removed to fully exploit the temporal shape and correctly infer physical properties of the event.A last filtering step in temporal space serves to mitigate high frequency noise brought by the deconvolution algorithm.An ultimate adjustment of the baseline completes the processing, leading to the final deconvolved time series S. Like for the simulation, a running sum of S provides the signal Sr.The comparison of a registered waveform corresponding to an α particle issued from 210 Po, with the simulation of the signal due to ionisation electrons created at the cathode surface, can be seen in Fig. 14.The small difference in shape was attributed to the slight deviation of the average distribution of migration times from GARFIELD++ with respect to the real behavior of the charges in the detector.
To graphically show the different event populations, an example corresponding to a run taken with the detector filled with xenon at 3 bar is shown in Fig. 15.The different sources ( 210 Po calibration source and 222 Rn contamination) are shown as well as their expected position.In the Dh versus Qt plot, the 210 Po structure is explained considering that 5.3 MeV α tracks are short (i.e.3-4 mm) and they traverse a small hole in the cathode to reach the active volume.For such short tracks the hole thickness of 200 µm is not negligible and in such a region the field is distorted, therefore not all ionization electrons will reach the central anode.This explains why longer tracks inside the volume (i.e.larger values of Dh) have a better energy resolution, whereas shorter tracks, partly absorbed in the cathode, exhibit a tail towards lower values of the reconstructed energy Qt.  7 Data analysis

Data-taking conditions
The primary goal of the R2D2 R&D is to establish the detector capability, when filled with xenon gas, of reaching percentlevel energy resolution.A second objective is the demonstration that the same energy resolution can be achieved with the same detector filled either with xenon or argon.This would be an important asset for the future of the project allowing the validation of larger prototypes with argon, in a much more economic way, and with no need of a gas recuperation system.As mentioned in Sec. 1, to reach the project goals it was decided to run the detector in ionisation mode.Such a choice is driven by several physical and practical reasons argued below.
-Provided that the signal-to-noise ratio is sufficient, operation in ionisation mode inherently yields the best possible energy resolution avoiding any stochastic fluctuation on the collected charge due to the avalanche process.Furthermore, tiny inhomogeneities of the wire diameter have no impact in ionisation mode, whereas in proportional mode they result in different gains as a function of the position along the wire impacting the energy resolution.-In the ionisation regime the signal is mostly given by the electrons drift which is similar between xenon and argon: the similar electron mobility allows for drift times of the order of 120 µs in xenon and 80 µs in argon.On the other hand, the proportional signal is mainly due to the drift of ions produced in the avalanche.With about three orders of magnitude smaller ion mobility in xenon, the time scales become very different for the two gases, and electronegative impurities are much more critical in xenon.-To work in proportional at high pressure a very thin wire or high voltages are needed.In the current setup the highvoltage power supply as well as the connections are validated up to 6 kV only.In addition when such a high bias is applied to the cathode a small leak current is observed, as well as an increase of noise.To overcome this issue several modifications would have been necessary, such as increasing the distance between the cathode and the grounded sphere, and replacing most of the connections.
A fine wire could nevertheless be used but at the cost of reducing the electric field at the cathode, making the measurements more sensitive to electronegative impurities of the gas.As explained in Sec. 4, the gas purity of the R2D2 setup is an important parameter and a weaker electric field would be unsustainable to operate without problems.
Data taken with two different anode thicknesses (50 µm and 1.2 mm diameter) were compared.The advantage of a thick anode lies in the ability to remain in the ionization regime and to increase the high-voltage (i.e., the electric field), thereby reducing the drift time and mitigating the impact of electronegative impurities.This is not the case for a thin wire since an increase of the high voltage could initiate an avalanche, bringing the detector into the proportional regime.For this reason the detector operation in ionisation mode at 6 bar with xenon was possible only with the 1.2 mm thick central anode.
Concerning the gas pressure, a scan between 1 and 10 bar was performed in argon, where the upper limit is given by the maximal operating pressure of the hot getter.The runs with xenon inside the detector were limited to 6 bar by the maximal volume allowed by the gas recovery system.

Event selection
To evaluate the energy resolution, two samples were selected, corresponding to α's of 5.3 MeV issued by 210 Po decays and α's of 5.5 MeV issued by 222 Rn decays.The first sample comes from the calibration source located behind the cathode as explained in Sec. 2. Such a sample has the advantage of being well localized with a known starting position of the α tracks.The second sample comes from radon which is introduced into the detector because of the recirculation through the cold getter.This is a well known issue of the cold getter for radio-pure experiments, but one that has turned into an advantage in providing a diffuse source for current testing.This decay allows the verification that the energy resolution is homogeneous throughout the detector volume, although is impacted by edge effects not corrected in the current setup (e.g.distortion of the field lines towards the end-caps).
The 210 Po events are selected in a Dt interval (dependent on the gas pressure and on the high-voltage applied) corresponding to events originating at the cathode radial position.Events rejected with a Dt below the lower accepted bound correspond to events originating in the gas volume closer to the anode, whereas events rejected (with Dt above the upper bound) correspond to pile-up events: those in coincidence with cosmics or low energy gammas.The selection is tuned for each run according to the gas and the high-voltage: for the run with the detector filled with xenon at 3 bar, the anode of 1.2 mm diameter and a high-voltage of -3000V (see Fig. 15(a)), events with a Dt between 100 and 110 µs are selected.A second selection for 210 Po events is based on Dh, which, as shown in Fig. 13(c), is an observable related to the radial length of the track.To minimize the impact of the energy loss in the hole through which α particles traverse the cathode, feature specific to the source location in the current setup, longer α tracks inside the detector are selected.This is done applying a lower cut on Dh in the data selection for the resolution analysis.Even in this case the cut depends on the run condition.For 3 bar of Xe (see Fig. 15(b)) events with Dh above 15 µs were selected.For the latter, the reconstructed energy is fitted with a Gaussian function that starts approximately halfway up the low energy side and continues until the end of the peak on the high energy side.All α's tracks pass through the cathode hole, but despite the selection cutoff Dh, a low energy tail due to partial energy loss persists in the selected data.An example is shown in Fig. 16.
The 222 Rn event selection follows the same logic used for 210 Po, with some additional treatments related to the temperature variation and the energy dependence according to the initial radial event position.Considering the radon event rate is approximately an order of magnitude lower than that  of the polonium source, to collect sufficient statistics the analysis duration should be at least 12 hours.As discussed in Sec. 5, variations of temperature at the level of 1 degree have an impact on the reconstructed charge at the percent level, thus degrading the energy resolution.The R2D2 setup is operated in an environment with variation of temperature between night and day which can reach up to few degrees.It follows that a correction based on the temperature variation becomes mandatory to obtain the desired energy resolution with the 222 Rn.Such a correction was systematically applied in each radon analysis.The second adjustment needed for radon events study aims to compensate for the position dependence of the reconstructed energy.As is seen in Fig. 15(a), events at small Dt (i.e.close to the anode) exhibit a small reconstructed charge, as perfectly predicted by the simulation presented in Sec. 6.Also, previous work showed that the dependence of proportional signals on the initial radial position could be well described by a power law [33].Although not optimal here, it was assumed that this dependency might still apply because both induction principles were analogous.Thus, part of the data was used to estimate such a function, which was then applied on the remaining data to correct for this effect.Once temperature and position dependence cor-  rections are accounted for, data are selected using Dt and Dh observables.An upper bound on Dt is needed to reject all events coming from the cathode and clean the sample from tails of the 210 Po source events.A lower cut on Dt is also applied to reject the region close to the anode, where the position dependence on the reconstructed charge is important and uncertainties on the model used for the correction could result in a resolution degradation.For the run with the detector filled with xenon at 3 bar, the anode of 1.2 mm radius and a high-voltage of -3000V (see Fig. 15(a)), events with a Dt between 70 and 100 µs are selected.An upper bound on Dh also helps to reject events from the cathode and clean the sample from polonium events.For the xenon run at 3 bar (see Fig. 15(b)) events with Dh below 12.5 µs are selected.An example of the fitted Gaussian function is shown in Fig. 17, where the low energy bound is almost symmetric to the high energy one.Indeed no source structure effect is present for radon, although a small tail of lower energy events can still be present due to detector edges effects.

Results
The obtained results are summarized in Fig. 18 and Tab.7.2.
The error at the level of 0.1% was estimated varying the selection cuts and the fit range.Nonetheless the detector running conditions in terms of vibrational and electronic noise were not constant between different runs.In particular, part of these disturbances was attributed to the electrical environment of the experimental room, which accommodates other measuring devices.This uncontrolled effect is probably the largest source of systematic error in the performed measurements and it could account for variations at the level of 0.2% (variations observed in runs taken in the same condition of gas pressure and high-voltage).The analysis performed on 5.3 MeV α's emitted by 210 Po (Fig. 18(a)) shows that a resolution at the level of 1.5% can be obtained independently on the gas pressure.Furthermore, similar results have been obtained for xenon and argon gases.This clearly suggests that the performance of future prototypes could be first established using argon, and extrapolated to xenon, at least if similar gas purities can be obtained.The thick anode seems to yield a better energy resolution, in particular at high pressure.The fact that the thicker anode produces a higher electric field at the cathode, thereby reducing the attenuation due to electron attachment during drift, could explain partly this observation.
The analysis performed on radon events confirms the independence of the energy resolution from the gas pressure and nature (i.e.argon or xenon).Switching from a well localized source on the cathode to a diffuse one does not degrade the energy resolution.These observations demonstrate that the energy resolution remains independent of the radial position of the energy deposit, despite a strongly varying electric field.However, a slight deviation between radon and polonium results can be observed.Unlike polonium events, the thinner wire yields a slightly better resolution for radon events at high pressures.A possible explanation is that for radon decays, happening throughout the volume and closer to the cathode, the attachment effect is weaker and the advantage of a thicker wire mitigated.Conversely, a thin wire results in a higher electron velocity at the end of their drift, and therefore a higher induced current ( I(r) ≃ v e (r)×E w (r)).Since the experimental noise is relatively constant, there is therefore an increase in the signal-to-noise ratio which could explain the improvement in the energy resolution.
The obtained results confirm that a CPC equipped with a non-polarized anode, combined with the ionization mode, offers very attractive optimization perspectives, whether in terms of polarization voltage, signal-to-noise ratio, or even anode thickness.Compared to the use of a detector equipped with a polarized anode where the noise strongly depends on the voltage (or a PPC which requires very strong polarization), the proposed CPC already offers more robust operating conditions.There is a balance to be found in the choice of detector parameters, and a reduction of the electronic noise could result in a further improvement in the energy resolution.

Conclusions
A cylindrical time projection chamber, equipped with a single anode and operating at high pressure, appears to combine many of the prerequisites necessary for a detector dedicated to the search for rare phenomena.In particular, one of the key points was to demonstrate that an energy resolution of the order of a percent could be achieved.In the present work the energy resolution was measured with α's of about 5 MeV in argon and xenon at pressures spanning between 1 and 10 bar.The presented observations validated a percent level energy resolution independently on the gas pressure and type, and on the initial deposited energy radial position.These achievements are also exciting to the extent that, with such a detector, detection adjustments could be made first under argon, before using xenon for the physics program.It is also worth mentioning that all the developments presented could be quite easily transposed into a CPC filled with liquid argon or liquid xenon.Finally, the obtained results pave the way to our final step in detector R&D for the search of the neutrinoless double beta decay program: designing a detector with the lowest possible material budget [36,37].

Fig. 2
Fig. 2 CPC set in the spherical vessel before final assembly.The wheel structure (white Teflon) supports the cathode (epoxy-copper type).The anode is attached to the central insert.

Fig. 3
Fig.3Evaluation of outgassing for some plastics polymers, potentially used for the CPC structure.

Fig. 4
Fig.4Evaluation of outgassing for different parts of the CPC and two cathodes: aluminium and copper covered with epoxy.The "Empty chamber" curve (blue line) labels the measurement of the vacuum chamber only, with no material inside.

Fig. 6
Fig. 6 Time evolution of the reconstructed average charge, corresponding to 5.3 MeV α's, for data taken at 3 bar with Xe.The stop and restart of the recirculation system are marked by the red and green arrow, respectively.

Fig. 8
Fig. 8 Reconstructed charge as a function of the electron drift time (obtained from 2023 GARFIELD++ release) in argon (8(a)) and xenon (8(b)).The different colours represents different gas pressures: 1 bar (blue), 3 bar (orange), 6 bar (green) and 10 bar (red).The dashed line represents the fit done using Eq. 1 and the obtained electron lifetime is shown on the plot.

7 .
This explains the large difference in the electron drift time (from a starting radial position at the cathode level), varying from a few ms in this study to 80/120 µs in argon/xenon, respectively, in standard detector operation.For each pressure the signal integral (i.e. the reconstructed charge) as a function of drift time was fitted with Eq. 1 to extract the electron lifetime.The results obtained at 1, 3, 6 and 10 bar in argon, with a drift time estimated using the 2023 release of GARFIELD++, are shown in Fig. 8(a).Similarly the results for xenon at 1, 3 and 6 bar are shown in Fig. 8(b).

Fig. 9
Fig. 9 Electron lifetime at different pressures for argon, based on the drift time computed with 2023 GARFIELD++ release (blue) and 2021 release (orange), and xenon based on the drift time computed with 2023 GARFIELD++ release (green).

Fig. 10 5 . 3
Fig. 10 5.3 MeV α signa integral (top) and temperature on the detector surface (bottom) as a function of time.

Fig. 11
Fig.11Evolution of the velocity of the primaries v (blue line) and the induced current I (green line) as a function of the distance r from the anode.For information purpose, the deduced mobility µ is also shown (orange line).The shown curves, built starting from the GARFIELD++ drift time simulation (t shown by the grey line), were computed for the CPC filled with a Xe gas at 3 bar and a high-voltage of 3000 V.

Fig. 12
Fig.12Example of a raw waveform (blue line), its deconvolved signal S (orange line) and the reintegrated signal Sr (green line).Note the scale factor of 30 for graphical purpose applied to the raw waveform and on the Sr signal.The horizontal solid red line represents the threshold whereas the vertical lines indicates the bounds for the signal total width Dt (dashed lines) and for the width at half maximum Dh (dotted line).

Fig. 13
Fig. 13 Signal modeling results.Qt (blue) and Dt (red) as a function of the distance r to the anode for a point-like energy deposit (13(a)).Dt versus Qt (13(b)).Dh versus the radial extension of the track, for tracks emerging from the cathode and heading towards the anode (13(c)).

Fig. 14
Fig.14Experimental waveform after processing for an α particle issued from 210 Po (blue) compared with a simulated signal due to electrons created at the cathode radial position (orange).

Fig. 15
Fig. 15 Example of a bi-plot Dt versus Qt (15(a)), and Dh versus Qt (15(b)) for a run with the detector filled with xenon at 3 bar, a central anode of 1.2 mm, and a high-voltage of -3000V applied to the cathode.

Fig. 16
Fig. 16 Resolution fit for 5.3 MeV αs from 210 Po source for a run with the detector filled with xenon at 3 bar, an anode of 1.2 mm, and a high-voltage of -3000V applied on the cathode.A selection of the events with 100µs < Dt < 110µs and Dh > 15µs has been applied.

Fig. 17
Fig. 17Resolution fit for 5.5 MeV αs from 222 Rn source for a run with the detector filled with xenon at 3 bar, an anode of 1.2 mm, and a high-voltage of -3000V applied on the cathode.A selection of the events with 70µs < Dt < 100µs and Dh < 12.5µs has been applied, in addition to the temperature and position corrections

Fig. 18
Fig. 18 Energy resolution obtained with 5.3 MeV α"s from 210 Po (18(a)), and 5.5 MeV α's from 222 Rn (18(b)).The blue lines represent xenon measurements whereas the red lines stands for argon measurements.Dotted lines are used for 50 µm anode whereas solid lines are used for 1.2 mm anode.Gas Pressure (bar) HV (V) Gain Anode diameter Resolution 210 Po Resolution 222 Rn

Table 1
Po Resolution 222 Rn Summary of all measurements taken, corresponding to Fig 18.