Prospects for measuring the 229Th isomer energy using a metallic magnetic microcalorimeter

The Thorium-229 isotope features a nuclear isomer state with an extremely low energy. The currently most accepted energy value, 7.8 +- 0.5 eV, was obtained from an indirect measurement using a NASA x-ray microcalorimeter with an instrumental resolution 26 eV. We study, how state-of-the-art magnetic metallic microcalorimeters with an energy resolution down to a few eV can be used to measure the isomer energy. In particular, resolving the 29.18 keV doublet in the \gamma-spectrum following the \alpha-decay of Uranium-233, corresponding to the decay into the ground and isomer state, allows to measure the isomer transition energy without additional theoretical input parameters, and increase the energy accuracy. We study the possibility of resolving the 29.18 keV line as a doublet and the dependence of the attainable precision of the energy measurement on the signal and background count rates and the instrumental resolution.


Introduction
The nuclear level scheme of the Thorium-229 isotope is expected to feature a long-lived isomer state, 229m Th, extremely close to the nuclear ground state.
The most recent value for the isomer energy E is , 7.8 ± 0.5 eV, obtained from indirect measurements with a NASA x-ray microcalorimeter1 [1,2], is within the reach of modern optical laser spectroscopy and could serve as a "nuclear frequency standard" [3].This standard could reach an uncertainty level of 10 −19 [4], and provide a new powerful instrument for testing the stability of fundamental constants [5,6].It has been shown that an ensemble of Thorium nuclei doped into a transparent crystal may demonstrate superradiance with a non-trivial emission dynamics [7], and may be used for building an ultraviolet (UV) laser [8].Finally, the frequency shifts and broadenings produced by such a crystal lattice environment might be used in studies of material properties, as is commonly done in Mössbauer spectroscopy [3].The necessary step towards all of these exciting applications is a direct observation and precise determination of the isomer state energy.
The existence of the low-energy state in the 229 Th nucleus was first conjectured by Kroger and Reich based on studies of the γ-ray spectrum following the α-decay of Uranium-233 [9].They concluded that this nucleus has a J π = 3/2 + isomer level lying within 100 eV above the J π = 5/2 + ground state level.The development of high quality germanium detectors (resolution from 300 to 900 eV) allowed Helmer and Reich to measure more precise γ-energies in 1989 -1993 and to predict the energy of the nuclear transition to be E is = 3.5 ± 1.0 eV, placing it into the range of optical frequencies [10].The decay pattern and combinations of transitions used by Helmer and Reich are presented in Figure 1 (a).
This unnaturally low value of E is triggered a multitude of investigations, approach discussed here using a high-resolution (∆ inst ≃ 3 -9 eV) microcalorimeter to resolve the 29.1 keV doublet (first proposed in [31]); (d): direct detection of "nuclear light" (many unsuccessful attempts [11][12][13][14] and new proposal [17]).Schemes (a), (b), and (c) are indirect measurements, involving keV energy transitions whereas scheme (d) is direct, only measuring the isomer transition of a few eV energy.
both theoretical and experimental, trying to determine the transition energy precisely, and to specify other properties of the J π = 3/2 + excited state (such as lifetime and magnetic moment).However, searches for direct photon emission from the low-lying excited state performed in the late 90's [11,12] have failed to observe a signal [13,14].In 2005, Guimãraes-Filho and Helene re-analysed the data of Helmer and Reich, taking into account new information about the nuclear decay pattern [15].They derived E is = 5.5 ± 1.0 eV.
In 2007, a cryogenic NASA x-ray microcalorimeter with instrumental resolution ∆ inst from 26 to 30 eV (FWHM) allowed Beck et al. [1] to perform a new indirect measurement of E is , involving lower energy nuclear states, as depicted in Figure 1 (b).In this measurement, the obtained transition energy (7 eV) was corrected by accounting for the theoretical branching ratios 29.19 keV → 229g Th estimated as 1/13, and 42.43 keV → 229m Th estimated as 2% in [2].This correction yields the currently most accepted value E is = 7.8 ± 0.5 eV, now placing the transition into the vacuum UV range (≈ 160 nm).
In the experiments decribed above ( [1,9,10]) the isomer transition energy E is is not measured directly but is derived from the spectrum of higher-energy (keV) γ-radiation of a spontaneously decaying 233 U source.We will refer to these measurements as indirect passive.Possible alternatives are direct passive and active approaches.
In the direct passive schemes (Figure 1 (d)), the aim is to perform spectroscopy of the ultraviolet radiation emitted from the isomer appearing in the α-decay of 233 U (2% of the nuclei decay is expected to lead into the isomer state).This method has two main difficulties: a relatively high false count rate caused by the Uranium sample radioactivity, and a high probability of nonradiative decay (quenching) of the isomer state in neutral Thorium atoms (up to 10 9 times higher than the radiative decay rate [16]).To overcome these problems, it was proposed in [17] to extract α-recoil Thorium ions ejected from an Uranium sample, and collect them in a small spot on a MgF 2 coated surface to minimize the quenching rate.Vacuum ultraviolet spectroscopy of the emitted fluorescence radiation may then allow to measure the isomer transition energy.
On the contrary, in active approaches, Thorium nuclei (in the ground state) will be illuminated by tunable radiation to excite them to the isomer state.In the solid-state approach a macroscopic (10 12 − 10 18 ) number of Thorium ions doped into UV transparent crystals can be excited, for example, by synchrotron radiation, and the emerging fluorescence signal can be studied [18][19][20].Apparent advantage of this approach is the huge number of simultaneously excited nuclei.At the same time, crystal fluorescence can cause difficulties in identifying the Thorium isomer transition, and various crystal effects can hamper the precise determination of E is .Another approach is the spectroscopy of trapped Thorium ions.At PTB, Germany, work is under way to excite nuclei of Th + ions into the isomer state using a two-photon scheme, exploiting the electronic bridge mechanism [21,22].In Georgia University of Technology, USA, the laser manipulation of Th 3+ ions is under investigation [4,23].Detection of the excitation of the Thorium into the isomer state may be based on a change of the electronic hyperfine structure [3].Studies of the hyperfine structure of Thorium are also performed at the IGISOL facility in Jyväskylä, Finland, in collaboration with a group of the University of Mainz, Germany [24].
We should also mention a number of studies aimed to measure the lifetime of the isomer state without a determination of E is .In [25], the half life of the isomer state for a bare nucleus was derived theoretically based on the calculations of the matrix element of the nuclear magnetic moment and on the experimental data concerning transitions at higher energies.They predict a half-life of T 1/2 = (10.95h)/(0.025E 3 ) for the isomer transition, where E is given in eV, which yields T 1/2 = 55 min for E = 7.8 eV.Direct measurements of this lifetime were performed in several groups [26][27][28].The obtained results vary from 2 min [26] to 6 hours [28].This discrepancy may be explained either by an incorrect interpretation of the observations [29] or by a difference in chemical composition of the Thorium resulting in different internal conversion rates.
All active approaches and eventually all nuclear spectroscopy applications require irradiation of the sample with some external narrow-band tunable radiation, and study of the emerging fluorescence.The estimation of the error σ = 0.5 eV on the isomer energy presented in [1] corresponds to one standard deviation, therefore it is necessary to scan the excitation source over 2 eV (±2σ) to find the transition with 95 % probability, or over 3 eV (±3σ) to find the transition with 99.7 % probability.Sakharov re-estimated the influence of the uncertainty of the 29.39 keV peak on the isomer energy derivation in [1] and obtained an error of 1.3-1.5 eV [30].Moreover, an analysis of more recent experimental data led him to claim that the energy of the isomer state can be anywhere in the range 0-15 eV, if the isomer state exists at all.
We believe it will be technically difficult, if not impossible, to cover such a broad energy range with a tunable narrow-band source of ultraviolet radiation in a reasonable time.We therefore propose to first increase the energy resolution on E is by an improved indirect measurement compared to [1].As we show below, it appears possible to resolve the 29.18 keV doublet [31] presented in Figure 1 (c) with todays state-of-the-art x-ray spectrometers.Resolving this doublet would significantly increase confidence in the existence of the isomer state.Moreover, the isomer energy would be measured without additional theoretical input parameters like branching ratios etc.The aim of the present study is to investigate the possibility of resolving the 29.18 keV line clearly as a doublet over a broad range of values for the isomer energy splitting and the branching ratio, and to analyze the precision that can be obtained on E is depending on the relevant experimental parameters.

Statistical aspects of the envisioned experiment
The operation principle of high-resolution x-ray microcalorimeters is to detect the heat deposited by an x-or γ-ray interacting with an absorber, using a very sensitive thermometer.Interaction with the absorber material mainly proceeds through the photoeffect.The energy of the produced photoelectron as well as Auger electrons together with their thermalization cascade should be effectively deposited within the volume of the absorber [32,33].On the other hand, the absorber should have a small heat capacity C a for good instrumental energy resolution ∆ inst .Various microcalorimeters differ in geometry, absorber material, sensor, etc., which leads to different energy resolutions, stopping powers, total detector surfaces etc.Many of these parameters are connected and can not be optimized independently.For example, increasing the size of absorber increases the solid angle and/or stopping power but degrades the instrumental energy resolution.Finally, we note that after a detection event, dissipation of the deposited heat leads to a detector-specific dead time, during which the energy of a successive photon can not be measured correctly.Therefore it is impossible to improve the precision of the measurements infinitely simply by using a more active sample, or by placing the sample very close to the detector.
The total count rate can be reduced using a designed filter which will primarily absorb photons outside the 29.18 keV region of interest.
The present study aims to answer two questions: how does the possibility to resolve the 29.18 keV peak as a doublet depend on the experimental parameters, and how to attain the most precise determination of the isomer transition energy E is .As outlined above, parameters of the experimental setup can be controlled to some extent only.For a proper design of suitable detectors and experimental configurations, we analyze how the key parameters affect these two points.

Specification of the problem and statistical model
The model employed for the statistical study should not contain too many parameters to make it accessible to a multi-factor analysis.On the other hand, it should be sufficiently comprehensive for a realistic feasibility study.For the sake of convenience, we assume a fixed total measurement time of t = 10 6 s, approximately 11 days, which corresponds to the total time of the successfull measurement in [1].Also we suppose that the background count rate near the 29.18 keV doublet is flat and symmetric, and that the monoenergetic line has Gaussian shape with full width at half maximum equal to ∆ inst (see section 4 for a discussion of these approximations).
The considered total energy interval (0-70 eV) is subdivided to a set of 0.4 eV bins (approximately a factor 10 below the expected instrumental energy reso-lution).The number of counts in the ith energy bin is a Poissonian random number n i with a mean value λ i equal to where R 29 is the signal count rate, r bg is a specific rate of background counts per 1 eV energy interval, E 1 , E 2 are the centers of lines of the components of the 29.18 keV doublet, , d = 0.4 eV is one energy bin, E i is the mean energy of the ith bin, and b is the branching ratio.The set {n 1 , ..., n N } of experimental data can be represented as a position vector n in an N -dimensional "sample space" (N = 175).
For a given sample n, we perform a nonlinear regression fit by a vector function f with N components This fit has 6 free parameters {J 1 , J 2 , Ẽ1 , Ẽis , rbg , σ} = {θ 1 , ..., θ 6 } ≡ θ.For the estimation of these parameters, we use the maximum likelihood method.
The likelihood function L(n| θ) is the probability for realizing the set n, if true mean values λ i are equal to f i (θ).We also introduce the logarithmic likelihood function

Resolving the 29.18 keV line as a doublet
As it was outlined above, resolving the 29.18 keV line as a doublet would significantly reduce the doubts [30] in the existence of the isomer state in the   The red spot corresponds to the area of the branching ratio b and isomer transition energy To check whether a given set of experimental data corresponds to a single line or to a doublet, one can apply the likelihood ratio test [34][35][36].The essence of this test and the method to estimate the significance level is described in Appendix A.Here we define the significance level as the probability to identify incorrectly the single peak as a doublet or the doublet as a single peak using the likelihood ratio test in the situation when we have either a single peak or a doublet whose parameters are specified.Figures 2, 3 show level curves for the signal count rates R 29 which are necessary to attain a significance level of 1 % for various values of the instrumental resolution ∆ inst and the R 29 /r bg ratio, in the (b, E is ) plane.It is interesting to note that the optimal energy resolution is attained when b ≃ 0.25 and not for equally strong components of the doublet.This is explained by the fact that such a branching ratio leads to a noticeable asymmetry of the peak which facilitates the identification of a second component, increasing b leads to a reduction of the main peak at fixed total signal count rate R 29 .

Precision in the determination of E is : Monte-Carlo simulations
The aim of the proposed spectroscopy study is not only to resolve the 29.18 keV line in the γ-spectrum of 233 U as a doublet but to determine the energy splitting with maximum precision.
We study the standard deviation of the isomer energy as a characteristic measure of precision (in the following we call δE is the uncertainty of E is ).Here and below, angular brackets denote expectation values, E is and Êis denote "true" and measured values of the isomer transition energy respectively.For the sake of brevity, we suppose that the true values of the energy splitting E is and of the branching ratio b are equal to 7.8 eV [2] and 1/14 [10] respectively2 .As before, we assume a total measurement time t = 10 6 s.Scaling the results to other values of E is , b, or t is straightforward.
To investigate the dependence of δE is on the instrumental energy resolution ∆ inst , the signal count rate R 29 , and the specific background count rate r bg , we perform a Monte-Carlo study of δE is .For any set of parameters, we simulate the sample n as shown in Figure 4, and estimate the parameters θ = {J 1 , J 2 , Ẽ1 , E is , rbg , s} maximizing the sum (3).We repeat this procedure 10 4 times and calculate δE is according to (4).
In Figure 5 we present curves of constant level of δE is in the plane (R 29 , r bg ) obtained from the Monte-Carlo simulation.Finally, in Figure 6 we present the curves of constant δE is in the plane (R 29 , ∆ inst ) for fixed ratios R 29 /r bg .
One can see that improving the instrumental resolution by 1 eV increases the precision on the determination of δE is by the same amount as doubling the 29.18 keV signal count rate, or doubling the measurement time t.

Experimentally attainable count rates and expected precision
In this section we estimate the attainable count rates and the resulting precision in a measurement of the isomer energy E is that can be achieved with a state-of-the-art high-resolution microcalorimeter.In particular we consider the metallic magnetic microcalorimeter maXs-20 as described in ref. [33].The aim is to demonstrate, that valuable results can be obtained with currently available technology.In section 5 we describe ongoing work towards a more refined, dedicated detector setup.
Each sensor is connected to a thermal bath, the system is installed in a cryostat operating at a temperature of about 30 mK.Energy deposited into the absorber plate heats the paramagnetic sensor and causes a change of its magnetization in an external magnetic field.Measuring this change in magnetization using SQUIDS, it is possible to determine the amount of absorbed energy and hence the energy of the incoming x-or γ-rays.Note that after the detection of an xor γ-ray, the individual pixel can only detect again after a certain relaxation time of the order of 100 ms.Therefore the total count rate R T should not be too high.
As a sample, we consider 1 mCi of 233 U electrodeposited as a film onto a metal planchet with a radius of R = 10 mm.We assume the sample to be situated 40 mm from the detector (outside the cryostat) with a total surface s = 0.5 mm 2 .Also we suppose the presence of additional material related to the cryostat vacuum system (sealing, input window of the cryostat, some other intentionally positioned shielding etc.), which we refer to as filters.
According to the NuDat 2.6 database [37], each single decay of 233 U is accompanied (on average) by one α particle with an energy from 4.309 to 4.824 MeV, 0.213 conversion electrons with energies from 2.3 to 600 keV (97.7 % of the elec-trons have energies below 50 keV), and 0.0544 photons most of which (0.052 per decay) are L-shell x-rays with a mean energy of 13 keV.As no individual L x-rays are listed in NuDat 2.6 (only average energy and total intensity), we have taken the lacking data from the X-Ray Data Booklet [38].
To estimate the detector count rates, we suppose that all α particles and electrons emitted from the Uranium sample are stopped by the sample itself or by the filter materials, therefore, we consider only x-and γ-rays.Also it is supposed that all secondary electrons and photons generated in the filters are absorbed in the material locally.This assumption is correct for relatively thick filters made from light materials like Aluminium.
We take into account absorption of the photons within the sample itself, the filters, and the detector.The total count rate R T is: Here the sum is taken over all photon energies E i , A is the activity of the Uranium sample, I i is the relative intensity (quantum output) of photons with the energy E i per single decay event, ℓ Au and ℓ f are the thicknessess of the gold absorber and filters respectively.Linear absorption coefficients a κ (E i ) (κ = U, Au, f) were taken from the XCOM Photon Cross Sections Database [39].
The count rate R 29 of signal photons is: where i = 1, 2 corresponds to the two components of the doublet, ̟ = 0.331 is the Au L shell fluorescence yield [40], I X is the probability that an energy of a fluorescence photon emitted by a Au atom is equal to E X , and P es (E γ , E X , ℓ Au ) is a probability that an incoming γ-quant with energy E i will be absorbed, and an x-ray photon following this absorption leaves the absorber (escape line).Sup-posing an isotropic spatial distribution of these secondary photons, we obtain where ℓ = ℓ Au , a γ = a Au (E γ ), a X = a Au (E X ).Also we suppose that the x-ray photon is emitted from the L shell, i.e. the deepest shell that is accessible by energy conservation, and the probability I X for emission of the photon is the relative intensity tabulated in [39] normalized to the sum of relative intensities from the L shell.
Evaluating the expressions ( 5) -( 7) yields a total detector count rate R T = 1 Hz and a signal count rate of R 29 = 7.74 mHz for a 1.3 mm thick Aluminium filter.Without any filter, the count rates for the same parameters are: R T = 13.6 Hz, and R 29 = 11.78 mHz.We see that the Aluminium filter absorbs approximately 92 % of all photons emitted from the sample, but only about 34 % of the signal photons.We conclude that filtering is an effective method to decrease the total count rate R T , caused mainly by low-energy Thorium L shell x-ray.
The background count rate r bg is caused by the escape of some fraction of the dissipated energy of γ-quants absorbed in the detector.In [1], the number of background counts close to the 29.18 keV doublet was about 30 -40 events per 3 eV bin whereas the total number of counts in the 29.18 keV peak was about 2.7•10 4 events.This yields the ratio R 29 /r bg = 2000 eV.Assuming that a similar ratio can be realized with the maXs-20 detector, we find that the uncertainty δE is on the measured isomer transition energy E is will be equal to 0.06 eV for an instrumental resolution of ∆ inst = 3 eV, signal count rate R 29 = 7.74 mHz and total measurement time t = 10 6 s.Therefore the proposed experiment to determine the isomer energy E is is expected to be almost one order of magnitude more precise than the results obtained in the previous experiment [1].Reducing the experimental resolution to 6 eV and 9 eV yields δE is = 0.19 eV and δE is = 0.56 eV respectively.Increasing the total measurement time to t = 2.6 • 10 6 s, 1 month, we can measure the isomer transition energy with an uncertainty δE is = 0.037 eV for ∆ inst = 3 eV, δE is = 0.12 eV for ∆ inst = 6 eV, or δE is = 0.33 eV for ∆ inst = 9 eV.

Further statistical aspects
We are aware of certain simplifications and assumptions in the above analysis.Here we briefly resume some additional issues that could arrise, a detailed discussion of these points is beyond the scope of this work.
First, the shape of signal peaks can deviate from Gaussian.For example, a long low energy tail on the spectral lines may lead to the appearance of a noticeable step in the backgound count rate (see, for example, Figure 2 (a) in [1]).We believe that in the work of ref. [1] this effect is caused mainly by the escape of energy from the absorber material, for example in the form of athermal phonons [33].The yield of Compton scattering is not sufficient to explain this step, see Appendix B. To take this effect into account correctly, we will have to modify our model of the background.A more difficult situation arises when the escaping energy is relatively small, which would lead to an asymmetry of the line rather than the appearance of a tail.In this case, it would be useful to study an isolated single peak separated from the doublet of interest but intense enough to give good statistics, and/or to perform an independent study with another γ source, e.g. 241Am.
Another issue that may appear is a slow time-dependent fluctuation of the response function caused by an uncontrollable drift of ambient magnetic fields and/or cryostat temperature over the duration of the measurement.We believe we can suppress such drifts below 10 eV by temperature stabilization and mumetal shielding of the setup.Additionally, we will monitor the position of a series of reference x-and γ-ray lines for a correct tracing of this drift, realizing a time-dependent calibration of the detector.Note that x-ray lines generally have a much broader linewidth than γ lines [42] which simplifies the identification.
An auxiliary calibration source, for example 241 Am, can help to enhance the quality of this callibration.
Also we should mention possible interference of the 29.18 keV doublet with coincidence and escape lines of x-and γ-rays of 223 U and other elements present in the sample.We plan to study the composition of the sample using "ordinary" low-precision γ-spectrometry.

Detector development
We are currently developing a dedicated new detector for the measurement described above, to some extend interpolating between the maXs-20 (0-20 keV) and the maXs-200 (0-200 keV) series [33].It will feature a linear array of magnetic calorimeters, each with an active area of 250 × 250 µm 2 .We will increase the absorber thickness by a factor of 2 to 3 in comparison with 5 µm in the maXs-20, leading to a stopping power of about 50 % at 30 keV.We will operate the detectors in a dry 3He/4He-dilution refrigerator at about 20 mK.In this situation, the intrinsic energy resolution of the detector caused by thermal noises of all kinds is expected to be below 2 eV (FWHM).According to the calorimetric detection principle of metallic magnetic calorimeters [41], this resolution is independent of energy as long as the total gain (including operational temperature and external magnetic fields) is stable, the dependence of the detector response on the event position in the absorber is negligible, and the photon energy E is still small enough to be within the range of linear detector response Φ(E).Also we expect that the minimal time between two correctly measurable counts in a single detector element (pixel) will be about 100 ms.
So far, we have achieved resolving powers up to about E/∆ inst = 3700 (corresponding to 1.6 eV (FWHM) at 5.9 keV) with our maXs-20 devices (unpublished), being limited by a combination of both, instabilities of the operating temperature and a position dependence.We believe that we can improve the short-term stability of the total gain and keep the position dependence of absorbtion events small enough to allow resolving powers beyond 10 4 in the planned experiment.Also, the response of the present maXs-20 detector to photon energies below 60 keV has a small quadratic deviation from a linear behavior, see Figure 7.At an energy of 30 keV this deviation is only about 3 %, which yields a 6 % degradation of the intrinsic energy resolution compared to the low-energy signals, i.e. below 2.12 eV on an absolute scale.

Sample preparation and characterization
The sample should ideally consist of isotopically pure 233 U to avoid a too high count rate not carrying relevant information and possible interference with the 29.18 keV Thorium doublet signal.For this project, we have 560 mg (about  To further purify the sample, we will perform a PUREX Uranium extraction procedure.We have also observed an efficient additional element separation in the electrodeposition process.In-house analysis using γ-and α-spectroscopy, ICP-MS, and neutron activation analysis will allow us to quantify the success of this procedures and finally know the exact composition of the final measurement sample. The sample will be produced by electroplating 233 U from a liquid solution onto a stainless steel or aluminium planchets.The target activity of 1 mCi cor-responds to 104 mg of pure 233 U or 118 mg of UO 2 .Producing correspondingly thick films (15-20 µm) turned out to be difficult in electroplating [43,44].We have therefore developed a process to deposit up to 20 mg Uranium onto stainless steel or aluminium foils of only 10-50 µm thickness.These samples can easily be stacked to realize the target activity without the carrier foils significantly reducing the count rate in the 29.18 keV peak.

Conclusion
We have analyzed the feasability of an indirect measurement of the lowenergy isomer state in 229 Th using a high-resolution magnetic microcalorimeter.
We propose to resolve the 29.18 keV doublet in the γ-radiation spectrum following the α-decay of 233 U.Such a measurement would provide a strong indication for the existence of the isomer state and improve the accuracy on the energy measurement significantly.The measurement appears feasable with currently available detector technology and samples.For the sake of brevity, we denote this random value as χ 2 r (a 2 ).It is easy to see that a is just a distance between the short hypothesis surface SS and the point y * = f (θ * ) corresponing to the true value θ * of parameters θ, see Figure 9.If the short hypothesis is true (it corresponds the situation when the surface SS coincided with SS ′ in Figure 9), then a = 0 and |L − S| 2 is just a "usual" χ 2 r random value with r degrees of freedom.
To test whether the short hypothesis is true or not, one should choose some desirable significance level α, and compare the value with some critical value λ α (χ 2 r ) such, that the probability If D > λ α (χ 2 r ), the short hypothesis is rejected on significance level α, otherwise it is accepted.The probability to reject the short hypothesis incorrectly is equal to α.On the other hand, the probability to accept the short hypothesis incorrectly is P (χ 2 r (a 2 ) < λ α (χ 2 r )).Therefore if the non-centrality parameter a 2 will be larger than a 2 r (α) such that the probability to accept the short hypothesis falsely will be less than α.
In our case, the observables n i are not normal random values with known dispersion but Poissonian random values.To estimate the possibility to identify the 29.18 keV line as a doublet, we consider √ n i as observables and approximate their distribution function by a normal distribution with mean √ λ i and dispersion 1/4.This approximation is not very precise (the bias is about 20 % for λ = 1), but it seems to be applicable for a coarse estimation.Then, for specific values of b, E is t, ∆ inst , R 29 , and r bg we calculate the set √ λ i according (1), and fit it by a vector function whose components are given by the square root of (2) with additional restriction J 2 = 0. Then we calculate the non-centrality parameter Our short hypothesis (single peak fit) has 4 degrees of freedom whereas the long hypothesis (double peak fit) has 6 degrees of freedom which yields r = 2.For significance level α = 0.01, conditions (10) and ( 11) are fulfilled at In ref. [1] the 29.18 keV line shows a low-energy tail, causing a noticeable step of about 10 counts per 3 eV bin (see Figure 2 (a) in [1]).This step can be caused either by the absorbtion of photons which had 29.18 keV energy initially, but lost some fraction of it in the source or filter material, or by the escape of some energy fraction from the absorber.The first scenario could be explained by the Compton "almost forward" scattering, when a photon loses tiny parts of its energy, about a few eV.The Klein-Nishina differential cross-section of the Compton scattering into the elementary solid angle is where E γ and E ′ γ are the energies of the photon before and after the Compton scattering respectively, r e = e 2 /(m e c 2 ) is the classical electronic radius, and θ is the scattering angle.Using a well-known relation between E ′ γ and θ, it is easy to express the differential cross-section in units of E ′ γ : Note that E ′ γ = E γ − 20 eV corresponds to θ ≃ 8.8 The differential probability for the 29.18 keV photon to be scattered into the 1 eV range can be coarsely estimated as where z i , n i , and ℓ i are the nuclear charge number, the atomic density, and the thickness of the Uranium (i = U) or filter (i = f) material respectively.
In the experiment [1], the Uranium activity of one 19 mm diameter planchet was about 0.02 mCi which corresponds to ℓ U ≃ 0.3 µm.It was covered by a Titanium foil ("filter") with ℓ f = 50.8µm.An estimation following (16) yields dP/dE ′ γ ≃ 1.9 • 10 −6 eV −1 .The peak corresponding to the absorption of 29.18 keV photons has Gaussian shape with about 3 • 10 3 events per 3 eV bin height in the maximum, and 26 eV FWHM, see Figure 2(a) in [1].This corresponds to a total number of counts forming this peak of about 3•10 4 .Therefore, Compton scattering of 29.18 keV photons in the source or filter material produces only 0.053 events per 1 eV interval, or 0.16 events per 3 eV bin.This value is significantly lower than the observed step of 10 counts per 3 eV bin.A similar estimation for a ℓ U = 15 µm layer of Uranium and ℓ f = 1.3 mm of Aluminium filter (parameters of the planned experiment) leads to dP/dE ′ γ ≃ 3.16 • 10 −5 eV −1 .Note that in (16) we treat all electrons as free.A more accurate estimation requires the substitution of the incoherent scattering function instead of z but this function does not exceed z, see [45] for details.Therefore, a more accurate calculations can only decrease the contribution of Compton scattering to the step of the background count rate.

Figure 1 :
Figure 1: Partial level schemes of the 229 Th nucleus with decay paths and energies (all in keV).Boxes in each panel denote the energy combinations used to derive E is in the "indirect" methods discussed in the main text.(a): according to Helmer and Reich [10]; (b): according to Beck et al. [1, 2], the interband transitions (dashed arrows) are taken into account; (c):

229
Th nucleus.Our first aim is to discuss the feasability of such an identification depending on experimental parametes and the (yet unknown) values of the isomer energy splitting E is = E 2 − E 1 and the branching ratio b.

Figure 2 :
Figure 2: Curves of constant levels of signal count rate R 29 (in mHz) required to resolve the 29.18 keV line as a doublet at 1 % significance level for different values of the detector resolution ∆ inst and the signal-to-noise ratio R 29 /r bg at 10 6 seconds of measurement time.The red spot corresponds to the area of the branching ratio b and isomer transition energy E is according to[2].

Figure 3 :
Figure 3: Curves of constant levels of signal count rate R 29 (in mHz) required to resolve the 29.18 keV line as a doublet at 1 % significance level for different values of the detector resolution ∆ inst and the signal-to-noise ratio R 29 /r bg at 10 6 seconds of measurement time.

Figure 4 :
Figure 4: Examples of Monte-Carlo simulated "experimental data" (black dots) and fit (red curves) for ∆ inst = 3 eV (a), 6 eV (b) and 9 eV (c).Other parameters are: R 29 = 7.74 mHz, r bg = 3.9 µHz/eV, t = 10 6 s.Plots are depicted in a "square-root scale" where the Poissonian noise is mapped onto signal-independent deviations.The origin of the energy axis is chosen arbitrarily.

Figure 5 :
Figure 5: Curves of constant δE is labeled in eV in the (R 29 , r bg ) plane for different values of ∆ inst and 10 6 s of total measurement time.

Figure 6 :
Figure 6: Curves of constant δE is labeled in eV in the (R 29 , ∆ inst ) plane for different values of R 29 /r bg and 10 6 s of total measurement time.

Figure 7 :
Figure 7: Detector response versus the energy of absorbed photons (a) and deviation from the linear behaviour (b) based on the measurement of 3 characteristic lines of an 241 Am γ-spectrum (c).

Figure 8 :
Figure 8: Left: Inductively coupled plasma mass spectrometer (ICP-MS) data of the raw 233 U material composition.Mass signals above 240 amu are molecular fragments and can be ignored.Right: Photo of the electroplated UO 2 test sample (with 238 U).

5
mCi activity) raw sample material available (in oxide powder form).The origin and preparation procedure of this material is unknown, from the γ-spectrum we suspect that originally, 232 Th has been activated in a high-flux neutron reactor and the 233 U has been separated chemically.A mass spectrum of the raw material, produced by an in-house ICP-MS can be seen in Figure8.The raw material contains > 90 % 233 U, together with traces of 232 U, 234 U, 235 U, 238 U, and the decay product 229 Th.Further daughter product of the 233 U chain have not been detected.

Figure 9 :
Figure 9: Sketch of the sample space: C and SS are the long-and short hypotheses surfaces respectively.Point y * = f (θ * ) corresponds to the true value of parameters, SS ′ is an s-dimensional surface passing through y * whose points are equidistant to SS. y is an experimental point, L = f (θ L ), and S = f (θ S ) corresponds the best fits within the long and short hypotheses respectively.The non-centrality parameter is the square of the distance a between SS and SS ′ .

λ 0.01 (χ 2 2 )
= 2 log(100) ≃ 9.21, and a 2 2 (0.01) ≃ 27.4.The values of R 29 necessary to attain a 2 = 27.4 for given values of the other parameters were calculated numerically and presented in Figures 5 and 6.Appendix B. The role of Compton scattering in the formation of an asymmetry or low-energy tail in absorption lines