Se170063

PAA_Selenium/Soil_Experiments#Selenium_Sample_Analysis

Sample Description

The sample was placed in an aluminum cylinder that was to be irradiated. The target components consisted of a nickel foil on the front of the cylinder with 2 pure selenium pellets under the foil, but still outside the cylinder. Inside the target there was burnt sagebrush ash, which was burned with a blowtorch, and selenium. Below are the masses of the components

Nickel Foil: 0.2783g

Outer Se Pellets: 0.0971g

Sage Ash: 0.5111g

Inner Se Pellets: 0.0523g

Energy

LB May Calibration 2017

The Calibration for Detector A was done on the morning of 5/23/17 with the MPA software using the thorium rods (as the calibration was fairly close already) and the correction values were found to be Det A Intercept = -12.208800 slope =1.021270

Efficiency

LB May 2017 Det A Efficiency

Activity and Half Life

The window used during this investigation was [110,118]. I used the script on daq1 /data/IAC/Se/May2017/5_25_17/Eff.C which fits lines of interest to a Gaussian function plus a constant. Once this function is found, the constant is integrated across the window of interested, then subtracted out of the original equation. After that the Gaussian is integrated across the window which gives the total number of counts for the peak of interest (in this case it is the 93 keV line of Se-81). To find the error in the counts, I increased the value for the standard deviation by a factor of 3 and integrated that function across the region of interest and subtracted that value from the original number of counts. Below is an example of a spectrum seen when using Eff.C

Once the raw number of counts and the error has been obtained, before entering the value into a .dat file a few modifications were made to that number. The first and simplest modification is to convert the number of counts into activity simply by dividing by the time duration of the measurement. In other words

[math]A(t) = \frac{N(t)}{t} [/math]

Since the radioactive decay equation is given by

[math]A(t) = A_0e^{-\lambda*t}[/math]

So to get a linear equation, simply take the natural log of both sides, which gives

(1) [math]\ln{(A(t))} = \ln{(A_0)} - \lambda*t [/math]

So we take the natural log of the activity and input that value into the .dat file. To propagate the error we take the error in the number of counts and divide by the time to put the error in units of Hertz. After that, since we are taking the natural log, we must use the standard error propagation formula, which is given by

(2) [math]\sigma_{\ln{A}} = \sqrt{(\frac{\partial \ln{A}}{\partial A})^2 * \sigma_A^2 } = \frac{1}{A}*\sigma_A = \frac{t}{N}*\frac{\sigma_N}{t} = \frac{\sigma_N}{N} [/math]

Once the activity and its associated error have been computed, put the values into a .dat file with 3 columns. The first column is the time elapsed for the measurement, the second column is the natural log of the activity, and finally the third column is the associated error found as outlined above.

Once the .dat file has been created, simply change the script /data/IAC/Se/May2017/5_25_17/LB_WeightedFit.C to read in the new .dat file and run the script on ROOT. Doing so will give a plot such as the one below

This plot gives parameters of the linear equation (1). The next step in the process is to find the half life of the sample. To do so recall that

[math] t_{\frac{1}{2}} = \frac{\ln{2}}{\lambda} = \frac{\ln{2}}{0.000207} = 3,348.54 Seconds \rightarrow 55.8090 Minutes [/math]

Now we must propagate the error in the half life. Even though the error in the decay constant says that it is zero in the plot, by looking at the root output, it can be seen that the error is actually

[math] 5.91816*10^{-9} s^{-1} [/math] (See if this can be fixed)

To propagate the error, use the standard propagation of error formula

[math] \sigma_{t_{\frac{1}{2}}} = \sqrt{(\frac{\partial t_{\frac{1}{2}}}{\partial \lambda})^2 * \sigma_\lambda ^2} = \frac{\ln{2}}{\lambda^2}*\sigma_\lambda = 0.095446 s^{-1} \rightarrow 0.001591 Minutes[/math]

which gives us a half life of

[math] t_{\frac{1}{2}} = 55.8090 \pm 0.0016 [/math] Minutes

The accepted value of the half life of Se-81 is

[math] 57.28 \pm 0.02 Minutes [/math]

Note that these two values do not overlap. This means that there could be some systematic error within the measurement systems.

Now, to find the activity of the sample when it was measured, we use the other parameter given by ROOT. Note that the value given is

[math] \ln{A_{measure}} [/math]

This means we must exponentiate to get the initial activity parameter to continue

[math] e^{\ln{A_{measure}}} = e^{5.34637} = 209.845 Hz [/math]

The associated error to the activity at the time of measurement can be found as well. Since we are exponentiating and the exponential function is its own derivative, the error can be written as

[math] \sigma_{A_t} = e^{\ln{A_t}}*\sigma_{\ln{A_t}} = 0.003 Hz [/math]

which means that the activity of the Pure Selenium sample at the beginning of the measurement was

[math] 209.845 \pm 0.003 Hz [/math]

Now note that the pure sample was measured second in the split run file. This means that the activity at the beginning of the measurement must be backtracked to when the mixture of ash and selenium was measured, which was 400 seconds prior. To do this, solve the radioactive decay equation for the initial activity.

[math] A_0 = A_t*e^{\lambda*t} = (209.845)*e^{0.00207*400} = 227.960 Hz [/math]

Now to propagate the error using the standard error propagation formula

[math] \sigma_{A_0} = \sqrt{(\frac{\partial A_0}{\partial A_t})^2\sigma_{A_t}^2 + (\frac{\partial A_0}{\lambda})^2\sigma_{\lambda}^2 + (\frac{\partial A_0}{\partial t})^2\sigma_{t}} = \sqrt{e^{2\lambda t}\sigma_{A_t}^2 + A_{t}^2t^2e^{2\lambda t}\sigma_{\lambda}^2 + A_t^2\lambda^2e^{2\lambda t}\sigma_t^2} = [/math]

Runlist

Table with dates and filename and locations on daq1

LB_May_2017_Irradiation_Day#10.25_Se_Soil_Mix

PAA_Selenium/Soil_Experiments#Selenium_Sample_Analysis