Demonstration of nuclear gamma-ray polarimetry based on a multi-layer CdTe Compton camera

To detect and track structural changes in atomic nuclei, the systematic study of nuclear levels with firm spin-parity assignments is important. While linear polarization measurements have been applied to determine the electromagnetic character of gamma-ray transitions, the applicable range is strongly limited due to the low efficiency of the detection system. The multi-layer Cadmium-Telluride (CdTe) Compton camera can be a state-of-the-art gamma-ray polarimeter for nuclear spectroscopy with the high position sensitivity and the detection efficiency. We demonstrated the capability to operate this detector as a reliable gamma-ray polarimeter by using polarized 847-keV gamma rays produced by the \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$^{56}\textrm{Fe}$$\end{document}56Fe(\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${ p},{ p'}\gamma $$\end{document}p,p′γ) reaction. By combining the experimental data and simulated calculations, the modulation curve for the gamma ray was successfully obtained. A remarkably high polarization sensitivity was achieved, compatible with a reasonable detection efficiency. Based on the obtained results, a possible future gamma-ray polarimetery is discussed.


Introduction
In nuclear physics, various exotic phenomena have been found by changing neutron and proton numbers.A breakdown of the neutron magic number N = 20 1,2 suggested the concept of the "Island of inversion".Shape coexistence in excited states were discussed in the neutron-rich "doubly-magic" nuclei 78 Ni 3 .Reflection asymmetry driven by octupole deformation produces parity doublet bands in the odd-mass nuclei 4 .Sudden shape changes recognized in neutron-rich zirconium isotopes 5,6 were interpreted as an occurrence of the quantum phase transition 7 .To detect and track these structural changes, the systematic study of the energy levels with firm spin-parity assignments is indispensable.However, the experimental methods have been strongly limited due to the available statistics of rare isotopes.In such circumstances, tentative assignments are often carried out by taking into account the transition strengths derived from experimental data, extrapolation of the systematics from less-exotic nuclei, and theoretical predictions.The development of efficient gamma-ray polarimeter, which can be applied for the spin-parity assignments in the wide range of nuclei, is expected to open up new perspectives in nuclear physics research.
Linear-polarization measurements of gamma rays have been applied to determine the electromagnetic character of the nuclear transitions.The differential cross section for Compton scattering of linear-polarized photons is described by Klein-Nishina formula 8 dσ dΩ where r e is the classical electron radius, E 0 /E is the ratio of energy of incident and scattered photons, θ is the polar scattering angle, φ is the angle between the polarization plane of the incident photon and the Compton scattering plane, respectively.
The formula represents the anisotropic dependence of the differential cross sections.While electric dipole (E1) radiation is dominant in atomic systems, electromagnetic multipole radiation often occurs in nuclear systems due to the actual level structure and the relevant nuclear wave functions.Therefore, linear polarization measurements in nuclear spectroscopy play an essential role for spin-parity assignments of the excited states.The analyzing power is known to be described as where E γ and E ′ γ are energies of incoming and scattered gamma rays.The analyzing power as a function of polar scattering angle θ is shown in Fig. 1.For sub-MeV gamma-ray polarimetry, detection areas which cover around 70 • are important to achieve high analyzing power.To obtain the angular distribution of the Compton scattering events, the detection system typically consists of two types of detectors, the scatterer and the absorber.In the early stage of the developments, the polarization sensitivity was achieved by combining scintillation counters 9 .
The high-energy resolution spectroscopy was achieved with different configurations of Ge(Li) spectrometers, such as in the separate crystals [10][11][12][13][14][15] , planer-type crystal 16 , and single crystal with electrically divided regions 15,[17][18][19][20][21] .Among those intensive research activities, the advent of the Clover Ge detector 22 , which contains four large-volume hyper pure Ge crystals, provided the reasonable polarization sensitivities and detection efficiencies 23,24 .By selecting separately the scattered events with the different crystals, the detector provided a simple identification of the transition character.Recently, gamma-ray tracking arrays such as the Advanced GAmma Tracking Array (AGATA) 25 and the Gamma-Ray Energy Tracking In-beam Nuclear Array (GRETINA) 26 have begun to provide a detailed scattered angular distribution of polarized gamma rays [27][28][29] .These detector arrays have a few-millimeter position resolution through pulse-shape analysis and signal decomposition.The high-position granularity has significantly improved the performance as a polarimeter without losing the absolute efficiency and the polarization sensitivity even when the detector is placed close enough to the radiation source.
Under these circumstances, CdTe-based detectors started showing its capability for Compton polarimetry.A fine-pitch CdTe-based polarimeter 30 showed a high sensitivity for 70-300 keV energy range.A two-layer configuration based on CdTe pixelized detector 31,32 indicated fine polarimetric potential of multi-layer CdTe detectors.Therefore, multi-layer position sensitive CdTe Compton cameras 33 serves as a new device for the linear polarization measurements in nuclear gamma-ray polarimetry.The detector is expected to offer the high-position resolution, high detection efficiency, and reasonable energy resolution for sub-MeV gamma rays.In this paper, we demonstrated the capability to operate this detector as a reliable gammaray polarimer by using polarized 847-keV gamma rays produced by proton inelastic reactions.By combining the experimental  data and simulated calculations, a modulation curve was successfully obtained.A remarkably high polarization sensitivity with a reasonable detection efficiency is reported.Based on the obtained performance as a polarimeter, a future-possible gamma-ray polarimetry is discussed.

Multi-layer position-sensitive CdTe Compton camera
The multi-layer position-sensitive semiconductor Compton cameras have been developed originally for astronomical observations on board the Hitomi satellite 33,34 .Detectors which employ both Si and CdTe layers 35,36 have brought the technological advances to various studies such as in the X-ray polarimetry from highly-charged ions on the ground 37,38 , the localization of radioactive materials dispersed following the nuclear power plant accident 39,40 , the multi-probe tracker in nuclear medicine and small animal imaging 41 .
As a gamma-ray polarimeter, the stacked CdTe layers are suitable 33,42,43 to achieve high detection efficiency for sub-MeV energy region.To demonstrate the capability of the gamma-ray polarimetry, a detector which is composed of twenty CdTe layers was employed for the present work.The detector and the schematic view of the arrangement of the sensors are shown in Fig. 2. The CdTe pad devices were produced by ACRORAD Co. Ltd.Each sensor has 8 × 8 pixels with the effective area of 25.6 × 25.6 mm 2 .Thus, the uniform position resolution of 3.2 mm was achieved over the effective area.The CdTe sensors were arranged in a 2 × 2 array for each layer.The thickness of each layer is 0.75 mm, and accordingly the total thickness amounts to 15.0 mm.The distance of each layer is 3.6 mm from 1st to 16th layers, 5.2 mm from 16th to 17th layers, 4.4 mm from 17th to 20th layers, respectively.The sensor configuration is similar to that described with more detail in reference 35 .A high voltage around 1000 V was applied for the operation.The Schottky-barrier diode type, which employed an aluminium (Al) anode and platinum (Pt) cathode, enabled us to apply the high-bias voltage with low-leakage current.The temperature around −18 • C were kept with a refrigeration system to prevent the sensors from the time instability under bias voltage (known as polarization phenomenon 44 ) and also to reduce the leakage current for better noise performance.
The charge signals from the electrode are amplified and converted to digital signals in the ASICs developed for Hitomi SGD 45 .Each channel has two charge sensitive amplifiers.One has a short shaping time (∼ 1 µs) for generating a trigger signal, and the other has a longer shaping time (∼ 3 µs) for holding the pulse height.The conversion time of the Wilkinson-type analogue-to-digital converters (ADCs) in the ASIC is less than 100 µs.The data were collected under the self triggering in zero-suppression mode.To subtract the common-mode noise level of each event, an ADC of the 32nd pulse height value (half the number of ASIC channels) was detected in every readout timing 45 .

Production and Detection of polarized gamma rays
The inelastic reaction, 56 Fe(p, p ′ γ), was selected to produce the polarized gamma rays.The reaction has been known to produce 847-keV gamma ray with nearly 50% polarization 12,23 .The schematic drawing of the measurement and the definition of angles are shown in Fig. 3.The target was placed at the origin.Θ is the polar angle to the beam axis.θ and φ are polar and azimuthal scattering angles.The natural iron foil of 10 µm was irradiated with the 3.0-MeV proton beam provided by the RIKEN Pelletron accelerator 46 .The detector was placed perpendicular to the beam axis, where the high degree of polarization is expected.The distance between the target and the front side of the first layer was set to 18.0 cm.The polar angle corresponded to Θ = 90 • ± 8 • .The beam intensity was appropriately tuned during the measurements to keep the total count rate around 1 kcps.In order to obtain the reference data, a segmented Ge detector (CNS-GRAPE 47 ) was placed on the opposite side of the Compton camera (Θ = −90 • ).The intensity of the gamma rays and the experimental asymmetry was monitored in parallel.

Energy spectra
The calibration parameters from each ADC value to energy of the gamma ray were determined from the measurements using standard gamma-ray sources, 137 Cs (662 keV), 133 Ba (81, 276, 303, 356 keV), 57 Co (112, 136 keV), 22 Na (511 keV), 54 Mn (835 keV) and 154 Eu (39.5+40.1,779 keV).The systematic uncertainties of the calibration based on the sources were estimated to be 0.1% up to 900 keV.The energy spectra with different multiplicity are shown in Fig. 4. The transition from the first 2 + to the 0 + ground state in 56 Fe is clearly identified at 847(1) keV.The energy resolution was 1.3% (FWHM) for the transition.There was no significant degradation on the energy resolution depends on the multiplicity.Low-intensity peaks at 352 and 367 keV were both assigned as the transitions from the excited state at 367 keV in 57 Fe, which is inherent in the iron target with the natural abundance of 2.1%.The average saturation energy of each pixel was estimated to be 890 keV, and about 10% of total number of pixels could not measure the 847-keV transition in a single interaction.However, the effect was negligible in Compton scattering events, which split the total energy in several pixels.The pixels which had the readout issue were also taken into consideration both in the data analysis and the simulations to compare them in a consistent way.The relative efficiencies from a single-hit event to 2-hit and 3-hit events to produce the full energy peak at 847 keV were estimated to be 63.8(3)% and 7.6(1)%, respectively.The 2-hit events which produced the 847-keV peak were used for the analysis.

Azimuthal distributions and degree of polarization
The degree of polarization of incoming gamma rays is obtained by comparing the azimuthal angular distributions of the measurement and the simulated calculation.The optimal azimuthal angle which satisfies the kinematics of Compton scattering is selected in the 2-hit event.In the procedure, the detector responses are normalized by the aid of the energy spectra in the early stage, and then the degree of polarization is estimated by the maximum-likelihood method.The modulation curve is obtained by taking the ratio of the measured and the simulated azimuthal angular distribution.The amplitude of the modulation curve indicates the polarimetric performance of the detector.We utilized the ComptonSoft toolkit 48 , which employed Geant4 simulation framework 49 for physical calculations.More details on the method can be found in a previous work on hard X-ray scattering polarimetry 37 .
First, we defined the divisions to minimize the effect of the noisy sensor regions and the background.After defining the division, we constructed φ histograms designated as d i j , where i denotes the bin number of the φ histogram and j denotes the index of division.The normalization factors ν j for each division were obtained by fitting the measured energy histograms (e (exp) j ) with the simulated histogram obtained by e (model) (P) = ν j [(1 − P)e (0) + Pe (100) ] + e (bg) , where P is the degree of polarization, e (0) and e (100) represent the simulated energy spectra with unpolarized, 100%-linearly polarized gamma rays and e (bg) represents the background energy histogram from the measurements.The background spectra were appropriately scaled by the live time of the measurements.For the iteration purpose, the degree of polarization was set to an initial value.Second, the logarithmic likelihood M(P) was computed with the fixed ν j given by where For here, the Poisson distribution (p i j ) is assumed for d . The procedure was iterated by updating the degree of polarization P until M(P) reached its maximum.In addition to the iteration process, the optimal polar scattering angle was investigated to increase the polarimeter performance.As a result, we found out that the restriction of the polar scattering angle (40 • ≤ θ ≤ 120 • ) significantly improved the polarization sensitivity with a certain loss of the detection efficiency.
The measured azimuthal angular distribution and the best-fit simulated results are shown in Fig. 5.The figure shows the source (black dots) and the simulated (red line) distributions and the ratio between experimental and simulated values.As a result, the degree of polarization was determined to be P = 0.57 (4).The positive sign of the degree of polarization is defined as an electric transition character, and the known electric character was reproduced 12,23 .The uncertainty of the degree of polarization was dominated by systematic one associated with different event selections.The systematic uncertainty originating from the detector rotation was also estimated by our simulation, and was negligibly small compared to the precision in the present work.A strength of this method is the deduction of the degree of polarization without obtaining the angular-intensity distribution with different polar angles 14 .

Modulation curve and Polarization sensitivity
The modulation curve for the 847-keV gamma ray is shown in Fig. 6.The curve was obtained by taking the ratio of the measured-and the simulated-azimuthal angular distribution with the obtained degree of polarization, as expressed: The curve was fitted with the function A(1 + Q ′ cos 2(φ − φ 0 )) where A, Q ′ , and φ 0 are the fitting parameters.Q ′ is called modulation factor, which is proportional to the degree of polarization.The modulation factor indicates the maximum experimental asymmetry observed in the modulation curve with a given degree of polarization.In the present work, we obtained Q ′ = 0.203 (9) with P = 0.57(4) from the fitting.The polarization sensitivity for 100% polarized photon is calculated to be Q = Q ′ /P = 0.35(4) at 847 keV.
To compare the polarization sensitivity with those of other detectors, the interpolation has been conducted by scaling the ideal sensitivity with a linear function , where a and b are fitting parameters from experimental data.The theoretical limit of the sensitivity (Q 0 ) is calculated by assuming the point-like scatterer and the absorber 20 as expressed: where m e c 2 is electron rest energy.For example, the sensitivities of GRETA and Clover detectors can be calculated as Q(847 keV) = 0.200(5) (a = 3.25(7) × 10 −4 , b = 0.131(6) 29 ) and 0.229 (a = 1.470 × 10 −4 , b = 0.3404 24 ), respectively.The similar sensitivity, Q(847 keV) = 0.32, was reported by POLALI 15 which adopted a low efficiency germanium crystal as a scatterer.Regarding polarimeters adopted CdTe sensors, a modulation factor 0.78 was reported in 200-300 keV region with a single layer, fine-pitch CdTe sensor 30 , and 0.13 was reported in 278 keV with two-layer configurations 31,32 , respectively.The estimation of the polarimetric performance of the energy range will be important for quantitative comparison.While the  (8), and φ 0 = 0.06(3), respectively.polarization sensitivity is an important design parameter for polarimeters, the sensitivity does not necessarily take into account the aspect of the detection efficiency.One can achieve a higher polarization sensitivity by imposing more strict constraints on the event selection.For example, the narrower polar scattering angle (60 • ≤ θ ≤ 100 • ) for the present data set provides a higher sensitivity of Q(847 keV) = 0.40(4), but results in the loss of 65% of the detection efficiency.
In position-insensitive detectors, the polarization sensitivity and the coincidence efficiency are typically in a trade-off relationship for a given detector volume and segmentation.To consider this aspect, a figure of merit 19,50 has been used as a quality description for polarimeters F = εQ 2 , where ε is the detection efficiency of Compton scattering events.The increase of the figure of merit results in a proportional decrease of the time to achieve required statistical significance.This efficiency is not necessarily identical with the intrinsic absolute coincidence efficiency, and here needs to consider the loss of the statistics by the event selection.With the use of ComptonSoft toolkit 48 , the absolute coincidence detection efficiency for 847-keV gamma-ray in the present setup was estimated to be ε coin = 3.0 × 10 −5 by considering the distance between the target and the detector, obtained the size of the radiation spot, and passive materials around the detector.For the present work, we imposed the restriction on the polar scattering angle (40 • ≤ θ ≤ 120 • ) to maximize the figure of merit.The polarization sensitivity increased by 1.5 times in exchange for a 36% the efficiency loss due to event selection.By considering these facts, we adopted ε ∼ 1.9 × 10 −5 for the efficiency.The figure of merit was estimated to be F ∼ 2.4 × 10 −6 .Note that the efficiency took into consideration the pixels which were not used for the analysis, but the materials such as the beam duct and the refrigerator were not considered.The performance is summarized in Table 1 with those of other gamma-ray polarimeters.The highest polarization sensitivity for 847 keV was demonstrated by the multi-layer CdTe detector with the comparable figure of merit.The strength of the position sensitive detector lies also in the capability to adjust its sensitivity for the transition of interest by optimizing event selection.

Possible future gamma-ray polarimetry
In the present work, a remarkably high polarization sensitivity by the multi-layer CdTe Compton Camera was demonstrated for the 847-keV gamma ray.In Table 1, one notices the significantly high figure of merit by the tracking Ge detector 29 .The high value is attributed to the high-detection efficiency of the array (ε = 5.3(1) × 10 −2 ).Increasing the absolute detection efficiency will be important for reliable gamma-ray polarimetry.Several possibilities for possible future gamma-ray polarimetry  29 .c interpolated by (a = 1.470 × 10 −4 , b = 0.3404: setup II) 24 .
First, achieving the large solid angle by modifying the present experimental setup improves the figure of merit.In the present work, a commercial refrigeration system to reduce the leakage current was employed.The system prevented the detector from placing close enough to the target.Minimizing physical conflicts can improve the absolute detection efficiency.If the distance between the target and the detector is shortened from 18.0 cm (current distance) to a few centimeter, more than 10-times larger solid angle will be achieved.
Second, adding the new CdTe sensors on the side close to the stacked layers will increase the coverage of the scattered gamma ray and analyzing power.The technology of stacking layers on the side has already been demonstrated such as in the Si/CdTe Compton camera 35 .In the present study, we found out that the restriction on the polar scattering angle (40 • ≤ θ ≤ 120 • ) increased the polarization sensitivity and maximized the figure of merit for the transition, compatible with a reasonable detection efficiency.Adding the side detectors will provide a better coverage for the polar scattering angle.
Lastly, the development of large volume of the CdTe sensors and high-position resolution are quite important for the next-generation gamma-ray polarimetry.With recent detector developments, CdTe layers with high-position resolution became possible with double-sided strips.As an example, the spacial resolution of 250 µm was recently reported 51 .By improving the position resolution, the close placement to the target will be achieved without the loss of the polarization sensitivity.
For an another perspective, the detector with ultrahigh position resolution will simultaneously provide the detailed polarangle intensity distribution, which has been conventionally measured by changing the angle of detectors or employing the large detector arrays.The intensity distribution independently provides the information about the degree of polarization at Θ = 90 • by fitting it with Legendre polynomials 14 .The development of high-position resolution detector will serve for a highly-efficient gamma-ray polarimeter not only for rare isotopes, but also for the systematic study of the degree of polarization, which would allow us to gain insight into the population patterns of magnetic substates in excited states.

Summary and Conclusion
We demonstrated the gamma-ray polarimetry with a multi-layer CdTe detector.The polarized 847-keV gamma rays were produced by the 56 Fe(p, p ′ γ) reaction at 3.0-MeV proton beam energy.By combining the experimental data and the simulated calculations, the modulation curve for the transition was successfully obtained.The highest polarization sensitivity was achieved for the transition with a reasonable detection efficiency.Based on the obtained performance, Compton polarimeter improvements for nuclear physics was discussed.

Figure 1 .
Figure 1.Analyzing power depends on the gamma-ray energy.

Figure 2 .
Figure 2. (a) Exterior of the multi-layer CdTe Compton camera.(b) Schematic drawing of the inside of the detector.Only sensor chips are displayed, and circuit boards and ASICs are not shown.

Figure 3 .
Figure 3. Schematic drawing of the measurement of polarized γ ray and definition of angles.The target was placed at the origin.Θ is the polar angle to the beam axis.(θ , φ ) are polar and azimuthal angles between the incoming and scattered photon directions.

Figure 4 .
Figure 4. Energy spectra with different multiplicity by the multi-layer CdTe Compton camera.The typical count rate of the 847-keV transition was around 3 cps under the 2-hit condition.The data acquisition time was approximately twelve hours.

Figure 5 .
Figure 5. (a) Measured azimuthal angular distributions for the 847-keV gamma rays.The black dots and red line represent experimental data and the simulated distributions, respectively.The normalization factors were applied into the simulation data before the cut of the polar scattering angle.The polar scattering angle was restricted to be 40 • ≤ θ ≤ 120 • both for experimental and simulated distributions.(b) Ratio between the experimental data and the simulated results.

Table 1 .
Comparison of the performance of CdTe Compton camera and other gamma-ray polarimeters.interpolated by (a = 3.25 × 10 −4 , b = 0.131) b