Received 7 June 2016; accepted 18 July 2016; published 21 July 2016
Potassium is one of the elements responsible for natural radioactivity. It is also a very important mineral for human life. It is found in soil and water, and stored in food like beans, fruit and cacao. Samples of the same kind of food with different provenance allow studying the effect of the soil on the potassium content. The radioactive isotope 40K has a natural abundance of 0.012%  . In 10.6% of the cases, it decays via a β+ decay to a metastable state of 40Ar. In the isomeric transition to the ground state, 40Ar emits a photon of 1460 keV, which can be identified by different kinds of detectors  . Measurements of the potassium content of food samples have been performed using Atomic Absorption Spectrometers  , Fluorimeters  , Sodium Iodide Based Detectors   and High Purity Germanium (HPGe) Detectors  -  .
HPGe detectors have an excellent energy resolution. Individual gamma lines up to several MeV can be measured with good efficiency, if the detectors are large enough. They are ideally suited to investigate natural radioactivity in food. The purpose of this study is to measure the potassium content of different food samples and to evaluate the possible effects of soil and water. The method used is gamma ray spectroscopy with HPGe detectors.
The energy spectrum shown in Figure 1 was obtained with a potassium chloride (KCl)  calibration sample. The peak at 1460 keV, has a FWHM of 2.03 keV, demonstrating the excellent energy resolution of HPGe detectors. All peaks expected from natural radioactivity are suppressed due to a lead shield and the strength of the 40K peak. Clearly visible is the 511 keV peak from pair production in the vicinity of the detector.
In Section 2, the experimental setup is described. A complete description of the preparation of the samples of food is found in Section 3. In Section 4, the method to extract the potassium content of the food samples from the spectra is described. The results of the analysis are shown and discussed in Section 5. The paper closes with some conclusions.
2. Experimental Setup
The experimental setup is shown in Figure 2. An HPGe Detector produced by Canberra is surrounded by a multilayer shield with a cavity to place the samples.
Figure 1. Energy spectrum obtained with a potassium chloride calibration sample in a shielded environment.
Figure 2. Experimental setup: the high purity germanium detector is surrounded by a multilayer shield leaving a cavity for the samples; (a) shows the open shield and (b) shows the closed shield.
2.1. High-Purity Germanium Detectors
The measurements described in the paper were performed with two different HPGe detectors  : a Reverse Electrode coaxial Germanium Detector (REGe) and an eXtended Range Germanium detector (XtRa). Both are commercial detectors produced by Canberra. The REGe is a closed-end n-type detector while the XtRa is a closed-end p-type detector. Both have a very thin contact on the front surface and a carbon composite cryostat window. Both detectors were cooled with liquid nitrogen. The technical specifications for REGe and XtRa are listed in Table 1.
2.2. Lead Shield
In order to be sensitive to the low level of radioactivity in food, it is necessary to reduce the environmental background. This was done with a multilayer shield composed of an outer lead shield with a thickness of up to 10.5 cm and an inner copper shield with a thickness of 0.4 cm (see Figure 2).
Figure 3 shows a comparison between energy spectra obtained with three different shield configurations. The side walls alone already substantially reduced the high energy gamma lines, especially the Thallium line at 2.6 MeV. This reduced also the Compton shoulder at medium energies. Figure 4 depicts the region around the 40K peak. The full shield reduced the 40K line by 95% and the continuous background in this region by 80%.
All the measurements presented in this paper were performed using a DGF Pixie-4 multichannel data acquisition system produced by XIA  . The software used to communicate with the Pixie-4 modules was a graphical user interface written using IGOR from Wavemetrics  .
Figure 3. Comparison of background energy spectra obtained with three dif- ferent shield configurations. Top: without shield; middle: with the open shield, i.e. lead bricks only as side walls; bottom with the closed shield.
Table 1. List of the specifications for the REGe and the XtRa detector.
Figure 4. Comparison of background energy spectra obtained with three dif- ferent shield configurations in the region of the 40K peak. Top: without shield; middle: with the open shield, i.e. lead bricks only as side walls; bottom with the closed shield.
3. Data Taking
Tables 2-5 provide an overview of all the food measurements discussed in this paper. The relevant quantities for each measurement are its lifetime and the mass of the sample. The date of the measurement is also given as the samples were purchased shortly before the measurements and variations with production time are probable.
The measurements listed in different tables were done with slightly different setups. A calibration measurement was done for each of these configurations. All samples were weighed to a precision of 0.1 g. The samples were placed in plastic containers. The container and its relative position to the detector were the same for all the measurements in a given configuration to guarantee the same geometrical acceptance.
Figure 5 shows pictures of different samples as used in the measurements. The first phase of the study (see Table 2) focused on fruit and nuts (F_12). In addition, some chocolates (DC_12) were investigated. In the DC_12 subset the % values indicate the cacao content of the dark chocolates as listed by the producer Lindt. Fruit and chocolate samples were cut in small pieces to fill the container as shown in Figure 5. The second and the third phases (see Table 3 and Table 4) focused on chocolate. The DC_13 and DC_14 subsets contain dark chocolates with different nominal cacao content produced by Vivani and Lindt, respectively. Using these samples, together with the DC_12 samples, the connection between the cacao percentage and the potassium content was investigated. The CC_13 and CC_14 subsets contain dark chocolate with the same nominal cacao content (75%) but different places of origin of the cacao. Using these samples, the effect of the place of origin of the cacao on the potassium content was investigated. The chocolate samples listed in Table 3 and Table 4 were chopped in order to obtain a uniform distribution of the chocolate in the container. This was not done in phase one. However, the effect of chopping was measured to be below 5%  and was neglected. The fourth phase (see Table 5) focused on legumes (L_14). Those were the samples for which the highest potassium content was expected. The effect of the self-absorption was estimated in a Monte Carlo study to be below 5% for the KCl, i.e. for the highest density of all the samples, and therefore neglected.
4. Analysis Method
4.1. Extraction of the Peak Area
The number of counts attributed to 40K was obtained by fitting the peak with a function F (E, σ, b) representing a skewed Gaussian plus a step function   :
Table 2. List of measurements from summer 2012 using the XtRa detector: KCL_12 contains calibration measurements, F_12 contains measurements of fruit and nuts and DC_12 contains measurements of samples, DCN_12, of dark chocolate with N% of cacao. Background measurements are listed in B_12.
Table 3. List of measurements from autumn 2013 using the XtRa detector: KCL_13 contains calibration measurements, DC_13 contains measurements of samples, DCN_13, of dark chocolate with N% of cacao and CC_13 contains measure- ments of samples of dark chocolate of 75% of cacao with different places of origin. Background measurements are listed in B_13.
Table 4. List of measurements from summer 2014 using the REGe detector: KCL_14 contains calibration measurements, DC_14 contains measurements of samples, DCN_14, of dark chocolate with N% of cacao and CC_14 contains measure- ments of samples of dark chocolate of 75% of cacao with different places of origin.Background measurements are listed in B_14.
Table 5. List of measurements from summer 2014 using the REGe detector: KCL_14 contains calibration measurements done with KCl, L_14 contains measurements of samples of legumes. Background measurements are listed in B_14.
Figure 5. Pictures of strawberries, pistachios and chocolate as used in the mea- surements.
(1a) describes the Gaussian component of the peak, where:
E is the energy;
A is the area underneath the peak;
σ is the standard deviation of the Gaussian;
μ is the center of the peak;
(1b) represents an exponentially modified Gaussian which describes the tail in the lower energy side of the peak, where:
Fskew is the skewness factor;
b = Wskew σ is the relaxation time of the exponential component;
(1c) represents the step function which describes the Compton scattering component, where:
I is the amplitude of the step function.
The uncertainty on A was obtained from the fitting procedure as implemented in the ROOT  analysis framework. Figure 6 shows an example of a fit performed for the chocolate sample CC06_14 listed in Table 4.
4.2. Background Subtraction
The lead shield does not completely remove the effect of the environmental potassium. The number of gammas coming from the sample was calculated as:
As the environmental background varies over time it is important to not only perform long enough background measurements, but to do them also in temporal proximity with the sample measurements.
4.3. Detector Calibration
The calibration of each setup was performed with KCl. The calibration factor was defined as:
Figure 6. Energy spectrum (histogram) of the chocolate sample CC06_14 from Table 4 around the 40K peak and the result of the fit performed using the func- tion in Equation 1 plus a first order polynomial represented by the parameters p0 and p1 (line).
where is the potassium content in the KCl sample in mg per 100 g, as specified in its data sheet, and is the area underneath the 40K peak in the KCl energy spectrum after background subtraction. The uncertainty on can be asymmetric. It is evaluated varying and according to their uncertainties, in order to obtain the one sigma deviations and. These uncertainties are taken as systematic uncertainties on the potassium content of the food samples.
4.4. Extraction of Potassium Content
Using the calibration factor, , the potassium content in mg per 100 g, CK, was calculated as:
where from Equation 2 is the area underneath the 40K peak measured for the food sample and is its mass. The statistical uncertainty is calculated from the fit uncertainty on and the uncertainty of 0.1 g on. The systematic uncertainty results from the uncertainty on ǫKCl (see Secition 4.3). The uncertainties are computed as asymmetric, but in many cases are symmetric and quoted as and.
5. Results and Discussion
The measured potassium content, , for all samples listed in Tables 2-5 are given in Tables 6-9. Also given are expectations,  , and statistical () and systematic () uncertainties. The statistical uncertainties dominate. The expected values  are reported without uncertainties.
Figure 7 shows the F_12 results on fruit and nuts listed in Table 6. The agreement between the expected and measured values is good for fresh fruit like bananas, kiwis and strawberries. The disagreement observed for processed food points to the influence of the processing techniques. The nuts show a potassium content lower than expected, while the dried fruit shows a higher one. The more dehydrated the fruit, the higher is the relative potassium content. The ratios between expectation and measurement are not uniform but the trends are very clear.
The subset DC_12 of Tables 6-8 list the results for chocolate. Given the cacao content, C (%), of a bar of chocolate of 100 g, the expected potassium content is calculated as:
Table 6. List of the expected () and measured () potassium content of the fruit, nuts and chocolates listed in Table 2. The measured value is quoted together with the statistical uncertainty, , the systematic uncertainty, , and their quadratic sum.
Table 7. List of the expected () and measured () potassium content of the samples of chocolate listed in Table 3. Other details as in Table 6.
where is the potassium content of pure cacao. Potassium is assumed to be contained only in
cacao. Figure 8 shows the correlation between the nominal content of cacao and the measured potassium content for samples produced by Lindt (DC_12 and DC_14, purchased in summer 2012 and summer 2014 respectively) and Vivani (DC_13), purchased in autumn 2013. The expected potassium content (Equation (5)) are also shown. Both Lindt chocolates have higher potassium content for higher percentage of cacao. However, the proportionality is not as expected: for lower (higher) nominal cacao content the measured potassium content is lower (higher) than expected. For Vivani chocolates, the potassium content is basically independent of the nominal % of cacao. This behavior could point to inaccurately quoted percentages of cacao or to production techniques influencing the potassium content of the product.
Figure 7. Comparison between the expected (squares) and the measured (triangles) values of the potassium content of selected fruit and nuts as listed in Table 6.
Table 8. List of the expected () and measured () potassium content of the samples of chocolate listed in Table 4. Other details as in Table 6.
Table 9. List of the expected () and measured () potassium content of the legumes listed in Table 5. Other details as in Table 6.
Figure 8. Correlation between the measured potassium content and the quoted percentage of cacao in bars of chocolate for different brands: squares show the expected potassium content while up triangles show Vivani samples from autumn 2013, down triangles Lindt samples from summer 2012 and circles Lindt samples from summer 2014 (see Tables 6-8).
Figure 9. Comparison between the expected (squares) and the measured (tri- angles) values of the potassium content in bars of dark-chocolate (75% of cacao), produced by Bonnat, with cacao coming from different countries as listed in Table 7 and Table 8. No measurement for Puerto Cabello was performed in autumn 2013.
Figure 9 shows a comparison between the expected and the measured potassium content in bars of chocolate produced by Bonnat with cacao coming from different countries. Measurements for each place of origin were performed with two slightly different setups and using two samples of chocolate (75% of cacao) purchased in different years. The results are from autumn 2013 (CC_13 in Table 7) and from summer 2014 (CC_14 in Table 8). The two sets of measurements agree within the uncertainties. However, the measurements show a potassium content significantly lower than expected. This behavior could point either to a cacao content lower than declared or to specific production techniques that might influence the potassium content. No significant variations due to the place of origin are observed.
Table 9 shows the results for legumes. Figure 10 shows the comparison between the expected and the measured potassium content. Both, red lentils and green peas, have a potassium content higher than expected, while the white beans autumn short. However, the white beans still have the highest measured potassium content.
Figure 10. Comparison between the expected (squares) and the measured (triangles) values of the potassium content for different kinds of legumes (see Table 9).
In this paper, a simple and clean method to measure the potassium content of samples of food with germanium detectors was presented. The results verify the high potassium content expected for fruit, nuts, chocolate and legumes. However, the expectations are not always met quantitatively. Especially the expected dependence of the potassium content of chocolate on the listed percentage of cacao could not be confirmed. No dependence on the origin of the samples could be established with statistical significance. The highest measured potassium contents are 1705 mg/100g from dark chocolate with 90% of cacao and 1630 mg/100g from white beans. Eating either 100 g of 90% cacao chocolate or 100 g of white beans provide approximately 50% of the daily recommendation of 3.5 g  of potassium.
We thank the entire GeDet group of the Max-Planck-Institute for Physics in Munich for their help and technical support. We would like specially to thank the students A. Burger, M. Duda and M. Schuster who helped carrying out some of the measurements. Part of the material presented in this paper is based upon work supported by the LMU-Mentoring founding program (autumn 2013).
 (2015) Potassium Isotopes. www.chemicalelements.com/elements/k.html
 Vaessen, H.A.M.G. and van de Kamp, C.G. (1989) Sodium and Potassium Assay of Foods and Biological Substrates by Atomic Absorption Spectroscopy (AAS). Pure and Applied Chemistry, 61, 113-120.
 Acquaha, G.F. and Pokoo-Aikinsb, M. (2013) Foodstuffs and Cancer: Analysis of Radionuclides and Its Radiation Levels in Common Ghanaian Maize. International Journal of Sciences: Basic and Applied Research (IJSBAR), 12.
 (2015) Wavemetrics. www.wavemetrics.com
 Phillips, G.W. and Marlow, K.W. (1976) Automatic Analysis of Gamma-Ray Spectra from Germanium Detectors. Nuclear Instruments and Methods, 137, 525-536.
 Mcnelles, L.A. and Campbell, J.L. (1975) Analytic Approximations to Peak Shapes Produced by Ge(Li) and Si(Li) Spectrometers. Nuclear Instruments and Methods, 127, 73-81.
 (2015) ROOT Data Analysis Framework. http://root.cern.ch/drupal/