Beta Ray Spectroscopy
From Physics 111-Lab Wiki
To Print this write-up, go to the wiki printable version, at the left side below.
The BETA RAY Experiment (BRA)
- Note that there is NO eating or drinking in the 111-Lab anywhere, except in rooms 282 & 286 LeConte on the bench with the BLUE stripe around it. Thank You the Staff.
Radioactive decay of nuclei by emission of beta rays (electrons) is a fundamental phenomenon of nuclear physics. In this experiment you will measure the energies and momenta of beta rays emitted by radioactive cesium. A study of the number of electrons emitted vs. the momentum yields information about decay processes.
Using a LabView Program window, you will setup the experiment settings, then a computer histograms the data. If you know what you are doing and make no mistakes, you can complete the data-taking in 4 days. Data analysis is not trivial, and you will need to plot a number of graphs, and analyze them. You will need to learn the nuclear physics of beta decay and internal conversion. The level is junior / senior, and you will need to read a fair amount about nuclear physics and radioactive decay. Get the reprints from the Physics Library. You will also learn about particle motion in magnetic fields, mass spectrometers, amplifiers, single-channel analyzers, digital scalers, counting statistics, computer control of data-taking, and how to handle radioactive materials safely. Once you set up the electronics, data -taking is computer controlled. You use Matlab to analyze the data. There are Matlab scripts written to help in the anlysis ( see Matlab fitting).
- Pre-requisites: None
- Days Allotted for the Experiment: 6
- Consecutive days: No
This lab will be graded 30% on theory, 30% on technique, and 40% on analysis. For more information, see the Advanced Lab Syllabus.
Comments: E-mail Don Orlando
The Beta Ray Spectroscopy Experiment Photos
Beta Ray setup with the Computer Click here to see larger picture
Beta Ray PMT Base Diagram Click here to see larger picture
Beta Ray Rack Experiment Click here to see larger picture
Before the Lab and SOP for Beta Ray
Complete the following before your experiment's scheduled start date:
- View the Beta-Ray Video.
- Radiation Safety SOP View the 'Radiation Safety Video'. After watching the video in the 111-Lab, get a pink Radiation Safety form from a 111-Lab staff person. Fill it out & sign the form for getting a Radiation Ring.
- Now complete the Radiation Safety Training Radiation_Safety After completion of Training turn in all forms to Don Orlando.
- Complete the BRA Pre Lab and Evaluation sheets. Print and fill it out. The Pre-Lab must be printed separately. Discuss the experiment and pre-lab questions with any faculty member or GSI and get it signed off by that faculty member or GSI. Turn in the signed pre-lab sheet with your lab report.
- View the Introduction to Error Analysis video and Error Analysis Notes.
- Read the Standard Operating Procedures (SOP) for this lab before starting SOP_3271_Cs-137_Na-22_Co-60_Mn-54_Am-241_Fe-55_2014.
- Tables for the Analysis of Beta Spectra; U.S. National Bureau of Standards‡
- Y. Yoshizawa, "Beta and Gamma Ray Spectroscopy of Cs137",‡ Nuclear Physics 5, (1958), 122-140.
- H. A. Bethe, Elementary Nuclear Theory.
- W. C. Haxton, B. R. Holstein, "[http://ajp.aapt.org/resource/1/ajpias/v68/i1/p15_s1 Neutrino Physics", Am. Jour. Phys. 68, 15 (2000).
- R. D. Evans, The Atomic Nucleus. McGraw Hill (1972).
‡ Contained information on the Fermi-Kurie plots and calculations. The reprints are all available on-line from the Physics 111-Lab Library Site.
You should keep a laboratory notebook. The notebook should contain a detailed record of everything that was done and how/why it was done, as well as all of the data and analysis, also with plenty of how/why entries. This will aid you when you write your report.
Some isotopes of atomic nuclei can decay spontaneously and may emit a combination of electrons, positrons, neutrinos, and gamma rays. Nuclei with atomic numbers greater than 80 and a few light isotopes can also emit alpha particles. In this experiment we study the momentum and energy of electrons, called beta particles or beta rays, emitted when radioactive Cesium-137 decays into Barium-137.
There are two modes of decay because the barium nucleus has two different energy states to which transitions are allowed. The nuclear equations, kinetic energies released, and the branching ratios (relative frequency of occurrence) of these two beta decays are written as follows:
Additionally, the excited state of Barium *, daughter product in (1), can subsequently decay as follows:
In (3), the decays by emitting a 0.662 MeV γ-ray, In (4), the excited barium nucleus decays by interacting with a nearby bound electron in a process called internal conversion.
The two beta-decays give a continuous energy spectrum of the emitted electrons, while the internal conversion reaction gives sharp k-, l- and m-peaks in energy (page 125 in the reprints). See Fig. 2a for a representative spectrum. It is a plot of the number of electrons N vs. momentum p. The curve is the sum of the three decays described above that emit electrons. The peak on the low energy side of the k-peak is an apparatus artifact.
What is Hysteresis?
In ferromagnetic materials, microscopic regions are separated into into domains. Within each domain, all the atoms have their magnetic moments aligned in one direction. Adjacent domains have their magnetic moments pointing in random directions with respect to their neighboring domains . A large scale sum of their magnetic moments is nearly zero, and the material produces only a small external macroscopic magnetic field.
Domains are separated from each other by "domain walls" which are 100 to 1000 atoms wide. Within the boundary of a domain wall, the individual atomic magnetic moments change directions from that existing within one domain, to the direction existing within the adjacent domain. When an external magnetic field is applied to the ferromagnetic material, the walls move and the domains with magnetization in the direction of the field increase in size and become macroscopic at the expense of adjacent domains that get smaller and disappear. The combined external field and the field of the ferromagnetic material can be orders of magnitude larger than the external field alone, when all the domains are aligned. This external field is generally supplied by a current-carrying coil wound around the ferromagnetic material.
The degree of magnetization as measured by the size of the domains is non-linearly proportional to the applied external magnetic field. Also, for a given external field, the magnetization of the material depends on the past history of the magnetization. For example, start from zero large scale magnetization, and zero external field; increase the external field to any particular value and then reduce it back to zero. The magnetization does NOT return to zero, but remains at some finite value. A reverse external field must be applied to bring the magnetization to zero. The magnetization lags behind the applied field. This effect is known as hysteresis.
The typical but exaggerate hysteresis behavior of a ferromagnetic material is shown in the figure below, where the Magnetic Induction B is plotted vs. the magnetic field intensity H, which is directly proportional to the current in the surrounding coil. The horizontal axis could equally well be labeled ``I" instead of ``H". "0" is the non-magnetized or demagnetized state (B=0, I=0) which can be reached by applying an alternating current and slowly reducing its amplitude to zero. As the current increases from the demagnetized state, the B value lags behind what one might expect and follows the curve to (a') and toward saturation at (a). With large variations in I, the behavior follows a major loop (abcdefa...); while for small variations in I, the behavior follows a minor loop (a'b'c'd'e'f'a'...) or (a'b'a'...) if the current I never becomes negative. Certain points of a hysteresis loop have names: Br is the remanent value at I = 0 after the material is saturated at (a.) Hc is the coercive field or the coercivity ( i.e. the H or I required to reduce B to zero).
The spectrometer is a 'C'-shaped vacuum chamber with magnet coils on top and bottom to produce a uniform magnetic field perpendicular to the radius of the chamber and perpendicular to the plane of the 'C'. A source of 137Cs is permanently mounted inside the chamber at one end of a rod see Figure 3a in the 'C' diagram, and the detector (see details in Figure 3b) is mounted at the other end through a slit (0.375 high x 0.125 wide) and a 2 mil mylar aluminum coated foil over a 1 inch diameter plastic scintillator with 1-1/2 inch diameter light pipe connected to a 6655A PMT. The overall chamber size is 0.4375 x 3.375 inches bent in a half circle shape.
Electrons emitted by the source are subject to the Lorentz force
and travel in a circle to the detector. The relevant equation is:
The radius r is fixed, and the field B is approximately proportional to the current I in the coils (see Appendix on hysteresis). We can therefore write:
We want to measure the momentum spectrum of the electrons, and we can do this by measuring the count N or number of electrons received at the detector in a given time interval, as a function of current I.
However, a correction must be applied to N because of the finite size of the detector aperture: The exit aperture is fixed in width and admits a given range of radii, Δr = constant. These radii cover a range of momenta Δp. The magnitude of Δp is not constant, but instead increases linearly with p, or the current, since p and I are proportional). The plot of N vs. p gives a true representation of the spectrum only if the increasing size of p is compensated or corrected for, by dividing each N by the current I. An illustration of what happens to the data is shown in Figures. 2a and 2b.
The computer controls the current in the spectrometer coils, and hence the field. It starts the current from approximately zero and increases it to maximum in discrete steps called bins or channels. For each bin, the computer counts the number of particles that strike the detector within a period of time and saves these data as a spectrum.
The detector consists of a plastic scintillator that produces a light pulse when an energetic electron strikes it. This light is detected by a photomultiplier tube (PMT) placed near the scintillator, as shown in Figure 3. The PMT produces an electrical pulse whose peak-voltage is proportional to the energy of the incident beta particle. The pulse is sent through a preamplifier, then a linear amplifier, to a single-channel-analyzer (SCA). The SCA generates a pulse for each detector signal that lies between the adjustable lower and upper thresholds. Because the SCA discriminates against noise pulses, it is also called a discriminator. The pulses are sent to the computer to be recorded in the proper channel of the spectrum. The signal can also be displayed on an oscilloscope.
You should know how the photo-multiplier tube works. Photo-multiplier tube Handbook‡ (Note this is a separate reprint booklet)
Include in your write-up an explanation of the PMT and what the pre-amplifier circuit might look like.
Photomultiplier Tube diagram
(See the RCA Photomultiplier Tube Handbook in the Physics Library)Photo-multiplier tube Handbook‡ (Note this is a separate reprint booklet)
We use a 6655A PMT for photon counting.
A Photomultiplier tube is a photon detector that converts single photons into large pulses of electric current by successive multiplying stages.
First read and get familiar with the LabView Beta Ray Computer Program [] which describe the operation of the computerized data-acquisition program that communicates with the DAQ card and controls data taking. The current is supplied to the magnet by a DC power supply in series with a current DVM meter through the relay circuits which can only be set to a particular value with this computer program.
Set the signal and discriminator levels
- Familiarize yourself with the relationship between PMT voltage with a maximum of -1100 VDC, amplifier gain, and pulse height: Consult the block diagram of the apparatus (Beta Ray Apparatus) and connect the output of the linear amplifier to channel 1 of the oscilloscope. The signal path is PMT to current PREAMP (SR570 - Bias Voltage : pos; Filter Type : none; Input Offset : pos, 2 x 100 nA; Filter Freq : not used; Sensitivity : 2 x 1 microA/V) to AMP to a 2 second [DELAY-LINE BOX to SCOPE (when the scope is triggering internally, the delay-line has no effect, but you will need it later to delay the main pulse]. Start with the following SCOPE parameters of CH 1: 0.2V, 2 microsec/div, CH 1 TRIG, and the Tran-L-Amp Amplifier settings: DIFF DL .8, INT 2, GAIN 16. Referring to a sample spectrum, set the coil current (Use the bin setting on the Integrate Single Point program) to a bin corresponding to a high count rate, (the K peak). Observing the signal on the SCOPE, gradually turn up the high voltage (negative polarity) on the photo-multiplier (PMT) until an operating voltage of -800 to -1100 volts maximum is reached. You should see many pulses of various heights. Vary the PMT voltage (always between -800 and -1100V max) and the amplifier gain and observe the effects on the signal gain and noise. Change the coil current to various points on the spectrum and observe how the intensity and height of the pulses change. Note what you see. Remember that you are trying for a maximum signal-to-noise in the pulses. Set the coil current to approximately 0 amps (bin 0) and observe the signal. Theoretically, there should be no signal at zero field. What do you see, and what can you conclude about the nature of the noise in this experiment?
You will need to use the "Signal Point Scan" for this part before seeing any signals from the appartatus. Refer to section setup of the Signal Point Scan Signal Point Scan
- Familiarize yourself with the operation of the SCA: The normal operation of a PMT will result in small-amplitude noise pulses. We want to use the discriminator on the SCA to reduce or eliminate the noise, permitting the signal pulses which are larger in amplitude to be passed and recorded. You can adjust the LLD Lower Level Discriminator or Baseline control to remove the noise line below the pulses. Connect the Tran-L-Amp linear amplifier to the SCA and the rest of the circuit elements as shown in Figure 3, apparatus Beta Ray Apparatus. Switch the SCOPE to EXT trigger. Now the scope is only displaying pulses that pass the SCA, and the 2 microsecond (-DELAY-LINE compensates for the SCA's signal transit time.
- With the Single Point Program, find a bin that corresponds to a high count on the sample spectrum. This program is a one part of two of the Beta Ray Scan Program Utilities, Beta Ray Computer Program []. Refer to that section now. Start with the SCA upper-level threshold (UL) set to 10 and the lower-level (LL) set to 0. There should be little or no signal displayed on the SCOPE. Increase the LL until a signal appears on the scope, this is the minimum LL setting you may use (approx. 0.36), the SCA will not operate properly at a lower setting. Now increase LL well past this setting and observe the output on the scope
- You should see something similar to figure 4.a. In this case, most of the low voltage noise signals are passing through the discriminator and registering as counts on the computer. Now, by slowly increasing the LL (Lower Level) setting, you should be able to obtain an output similar to Figure 4.b. In this picture the lower line and a few of the lower voltage pulses beneath the brightest pulse have been eliminated. This effectively eliminates most of the PMT noise as well as some of the lower momentum pulses..
- The object is to set the LL high enough to eliminate small-amplitude PMT noise while preserving as much of the low-momentum data as you can.
- The AMP gain should be set around 16, but be careful not to clip the top off your signal. If you set the gain to high on the amplifier it will chip or saturate and look like figure 4C. Now using the Beta Ray Scan program, take a series of quick spectra (Bins 0-2500 by 10; 3 sec/bin; save time by clicking on the Stop Early Button, click only once, after it has completed taking data in the 'up' direction at about channel # 2450 (note that the Stop Button takes about 40 second to really stop the program). This should verify your SCA and amplifier settings. With higher LL settings, you should see the lower momentum portion of the spectrum erode. Compare your spectrum to the sample spectrum. If your spectra look too distorted, you may have to repeat step three. Why don't the near-zero-field data disappear with higher discriminator settings?
Observe the effects of statistical fluctuations.
- Using the Integrate Single Bin program acquire 50 or 60 observations each, at one second per integration, for two bins with greatly different count-rates. Save the observations in a file, and use them to calculate the means and standard deviations for the two channels. How do they compare to counting statistics?
- For one bin (one single setting of coil-current), observe the change in fluctuation as you change integration times in a range from 1 to 30 seconds. Based on the lowest count-rate of all channels, for how long should you integrate each channel of your spectrum to achieve a better than 1% error due to statistical factors?
Observe the effect of hysteresis
- Re-read the [hysterisis section] and the Analysis section of VI. Beta Ray Computer Programs to be sure that you understand how the scan program establishes a hysteresis loop before it begins taking data.
- Take a quick complete spectrum (Bins 0-2000 by 100; 2000-2500 by 5; 2500-4000 by 100; 2 seconds/bin). Observe the shift in the k-peak between the two directions of current. Using the knowledge that the k-peak always occurs at a given magnetic induction (not at a given current), convince yourself that the observed shift in the spectrum between current up and current down is in the correct direction.
Take data using the Beta Ray Scan Program to scan the entire beta ray spectrum
The current in the magnet is scanned up and down which equals one complete scan. The computer controls the past history of the scans so as to follow the same hysteresis curve every time every time. In adjusting scan parameters, it is more important to integrate for longer times, 60 seconds maximum, than it is to include every point in the spectrum. Print out copies of the scans. We are networked in the 111-LAB so you can choose where to print and what printer to use.
Beta Ray Computer Program
With the aid of the computer and a LabView program, you can get a rough spectrum in 44 minutes and a full spectrum overnight.
The computer is equipped with an interface DAQ card, PCI-MIO-16E-4, that enables the PC to vary the coil current in either direction from zero to a maximum of approximately 650 mA in 4096 steps. The computer controls the coil current by setting a DAC (digital-to-analog converter) in prescribed steps from 0 to 4095. The DAC outputs its voltage to the power supply that is operating in remote-voltage-controlled current mode. The interface DAQ card also has a edge counter, for counting the beta-particles and background events that strike the detector, and digital TTL-output lines, that permit the PC to switch the current direction. Taken together, these pieces form a 4096-channel particle-momentum spectrum analyzer. The program that operates it is a Labview program entitled BetaRay Utils.exe". It will open the two main programs "BetaRay Scan V2.2" and "Single Point Scan V2" but it will NOT start them. Each of them should be run alone. Wait till one stops before running the other one. If you run both together the LabView will crash. On the Beta Ray computer, this program is accessed by the shortcut on the desktop. If it is not there, look for it under Programs before contacting the 111-Staff. Data will be saved in your My Documents folder with a pre-name of "BRA Data" and "Single Point" you should add your run # after this Pre-Name. (ie; BRA Data-1 or Single Point-1)
First, it's helpful to define some terms and outline how the system operates. The program can only set the coil-current to fixed values that are determined by the DAC and the power-supply. Each value of current is called a bin or channel. The computer degausses the magnet, sets the current to a bin and then starts counting the number of pulses coming from the detector. This is called integrating a bin. After a period of time, the program stops integrating, takes the number of particles it counted, and moves on to the next bin. After determining the counts for every bin, the program produces a file with entries of Channel Number and counts which constitute a raw spectrum (raw since it is uncompensated for hysteresis effects). The Single Point Program output is Time and Counts. The full Beta Ray Scan program data output file contains Raw Data or Summed Data with the run number and counts. Single Point Scan Program: when you start the program you need to input the channel number first before clicking the start button in the program. The program will then ask you if you are saving the data and for a file name, ie; Single PointXX (you input the XX variable as a run number). Then input the other parameters needed to run the program:
- Base File Path: where your file is being saved and what name the computer gave it.
- Status Ring: information on what the program is doing at the time.
- Integration Time per Bin: 60 is the maximum number.
- Time/Scan Direction: is the time will take to complete the requested up & down scan.
- Estimated Time Remaining: it is updated at the end of each up/down scan with the remaining time.
Note that data is only saved at the completion of each up and down scan. You will see the summed data plot after this time.
This is a good time to point out how your spectra will differ depending on the direction of the changing current. Look at the figure 5B below.
The K-peak of the beta spectrum occurs at a specific field strength, marked Bk. If B were directly proportional to the coil-current, the K-peak would always be detected at current ik. Instead, as you take data with the current increasing, the peak occurs at iku, and as the current decreases the peak occurs at ikd. This is true for every point of the spectrum, which means that both spectra you take will be shifted away (in opposite directions of I) from the spectrum you would get if there were no hysteresis. Figure 5C is closer to the real hysteresis curve for the beta spectrometer magnet.
The Beta Ray Scan Program:
Beta Ray Scan is a program that acquires, saves, and plots the number of counts vs. current (bin number), according to the user's input parameters.
Your data will be saved on the 111-Lab server backed up every night by the campus. Data is saved in your own home directory in \\Atlas2\redirect$\your_Directory_name\My Documents\LabVIEW Data\then either "BRA DataXX" or Single PointXX" .
Note you should use Matlab for all your data anlaysis. We have written scripts to help you in you anlaysis. Keep your file names short and desciptive, eg; dataup-1 ending in (.dat), then your data can be read by many other program editors for any changes needed.
The two programs are located on the your desktop. Now at this time please do not run and excute both prgrams at the same time. They will stop working. Do not close them either, just stop them.
After you have started Beta-Ray Scan, you will be asked if you want to save the data that you will be taking. If you choose not to save the data, the data will only be plotted. If you choose to save the data, you will be asked for a file name with the Base File Path as LabVIEW Data\BRA Data or Single Point. The program will save raw data points as they are gathered with the file paths:
[Base File Path] DOWN data (Raw Data Run # [Run #]).xls
[Base File Path] UP data (Raw Data Run # [Run #]).xls
For the Up and Down scan data. The Summed Data will be saved intermittently to the same file with the paths:
[Base File Path] DOWN data (Summed Data).xls
[Base File Path] UP data (Summed Data).xls
For the SUM UP and SUM DOWN data. (Summed data are just that-they are the direct sum of all the previous runs. Think of adding vectors where each component represents a bin number, and the magnitude of that component represents the total number of counts for that bin number).
If you forget where you chose to save the data, don't worry. The Base File Path indicator displays the Base File Path.
Once you have decided whether or not to save the data, you will then have to input the relevant parameters. Here's a list of them and what they do:
- Integration Time per Bin(60s Max): Controls the amount of time (in seconds) the computer will sit on each bin tallying the number of counts; e.g., if you choose 5 seconds, the computer will count at each bin for five seconds. Maximum of 60 second for this variable.
- Delta Bin: The number of bins the computer will skip (won't integrate) between bins that it integrates.
- Number of Up/Down Scans to Make: Each scan consists of data taken as the current increases (UP data) and data taken as the current decreases (DOWN data). And one SUM Total Data file. You should make as many Scans as possible for good data statistics within your time constraints.
- Scan Quality: Determines whether or not the computer will wait for the current to settle as it changes bins (while it's taking data). If you select Scan Quality to be HIGH, the computer will wait 4 seconds for the current to settle before it starts integrating the new bin. If you select Scan Quality to be LOW, it won't wait at all. If you want good data, choose a HIGH scan quality.
- Scan Start Bin #: Controls the bin number which the scan will begin taking data at. (The same hysteresis curve will be swept out regardless of the Scan Start Bin.)
- Scan End Bin #: Controls the bin number which the scan will stop taking data at. The actual end bin number is the closest one to the chosen Scan End Bin # (rounding down) so that the quantity
is an integer. (The same hysteresis curve will be swept out regardless of the Scan End bin #.)
- Time/Scan Direction: Displays the projected time required for each scan (up and down) to be completed (based on the user inputs) in the format hh:mm:ss. For example: 1:23:36 requires one hour, twenty-three minutes, thirty-six seconds to complete a scan.
- Estimated Time Remaining: Displays the projected time required to complete the entire scan sequence (i.e. Time/Scan x # Scans to Make). This indicator updates itself at the end of each scan completed, so it will display the projected time required to finish the total scan sequence (not counting the time required to complete the current scan).
NOTE: Ten scans at 5 sec/bin is better than one scan at 50 sec/bin. (Sample values would are below:
Integration Time per Bin (Max of 60 seconds) = 5
Delta Bin = 5
Scans to Make = 7 (could be 6)
Scan Quality = LOW or HIGH add (1/2 second per bin integrated)
Scan Start Bin = 100
Scan End Bin = 4095
Time/Scan = 01:07:01
Total Time = 15:38:14 this gives you an overnight run.
Once all the desired parameters have been entered, click on the flashing [START TAKING DATA] button. Now everything should go as planned. To know what the computer is doing at any given moment, keep an eye on the Status Ring.
The basic outline for the procedure is as follows:
- The current is set to 0
- The current is set to its negative maximum.
- The current is returned to 0.
- The current is ramped up to that corresponding to the Scan Start Bin #.
- The UP data are taken, stopping at the Scan End Bin # or thereabouts.
- The current is increased to its maximum value.
- The current is decreased to that corresponding to the Scan End Bin # or thereabouts whatever the last integrated bin was in the UP scan.
- The DOWN data are taken, stopping at the Scan Start Bin #.
- The current is returned to zero.
- The raw UP, DOWN, SUM UP and SUM DOWN data are saved/updated and plotted.
- The entire procedure is repeated until all the scans have been made.
During the scan process, you may stop the program pre-maturely by clicking on the [STOP EARLY] button (it will appear after the data run(s) begin). Just click it once. If it seems like it's stuck and won't click, that's probably because the computer is busy integrating a bin and hasn't had time to click it. When it's done integrating that bin, the [STOP EARLY] button will latch. The program may integrate one more bin, then it will stop and bring the magnet through the end of the hysteresis loop once or twice (depending on how many runs were left to complete it takes about 78 seconds to stop). You must let it stop normally or the equipment will be left in an unstable condition. Data taken during the interrupted run will be graphed but not be saved or added to the sum data. (the graph of interrupted down data isn't to be trusted, as the bin numbers won't actually correspond to the points integrated).
Note: Since statistical error decreases with larger samples, longer integration times will give you smoother graphs. Below are two examples of the effect of integration-time on the measured spectrum:
You can see from these graphs (even though their small vertical size slightly obscures their resolution) that the longer the integration time, the smoother the graph. On the other hand, it is not always desirable to take too long an integration time. Sometimes instrumental conditions change over time, and if the last point taken is too far in time from the first point, then the curve does not truly represent the experiment. Ten scans at 5 sec/bin is better than one scan at 50 sec/bin. Take many scans at a lower integration time.
The integrate single point program:
Beta-Ray Single Point Integration is a program that will sit at a single bin number, magnet current, (determined by the user) and periodically integrate that bin number for a certain amount of time. With this program, one can see how much (if any) instrumental drift there is over long periods of time.
Starting from the computer desktop:
- Double click on the BetaRay Utils (Both VI programs will open, but NOT start running; Single point Integration V2 and Beta Scan V2.2)
- Click on the RUN button at the top of the LabView toolbar after you have entered the start bin #.
The following are sample parameters to use:
- Program VI = Beta Ray Single Point Integration
- Integration Time per point = 10 (seconds
- Time Between Points = 1
- Integrate Bin # = 2200 ( the K-peak, enter it before start Run button is clicked)
- Number of Points to Take = 60
Click on enter ( upper left corner) Click on the Run button Question: to SAVE your data or Not to SAVE, your choice.
TOTAL TIME = 00:11:33 (11 minutes and 33 seconds) time enough for you to set the LLD on the discriminator.
To Stop the Program Early, double click on the "Stop Early" Button (Note: Button will change to Stopping, wait about 40 seconds for it to stop) .
After you have successfully started Beta-Ray Single Point Integration, you will be asked if you want to save the data that you will be taking. If you choose not to save the data, the data will only be plotted. If you choose to save the data, you will be asked for a Base File Path Name. The program will save the data after the run with the following path:
[Base File Path] Single Point-# Bin Number 2200.xls
If you forget where you chose to save the data, don't worry. The Base File Path indicator displays the Base File Path. Once you have decided whether or not to save the data (and where), you will then have to enter the scan parameters. Note that all data saving defaults to the your \\~\My Documents\LabVIEW Data\ directory.
- Integration Time per Point: Sets the amount of time the computer integrates for each data point.
- # Points to Take: Sets the total number of points that will be taken.
- Time Between Points: Determines how long the computer will wait after integrating one point before integrating the next point.
- Integrate Bin #: Selects which bin number will be integrated.
- Total Time: Displays the total time the procedure will take in the hh:mm:ss format.
For example: Total Time = 11 minutes and 25 seconds ( 11:25)
- Stop Early: If this button is clicked on, the scan will stop in the middle of the data gathering part (at the end of its integrate/wait sub-cycle), the hysteresis loop will be swept out, and the data will be saved.
The Status Ring will tell you what the computer is currently doing. The hysteresis cycle that the computer sweeps out is identical to that in the Beta-Ray Scan program. Data will be plotted as they are taken.
Please see the Matlab program for analysis and the Matlab scripts that have been written for this experiment. The scripts and the executable of Matlab should be in the same working sub-directory so that they can be called inside Matlab; otherwise, you have to use the Matlab command ADDPATH to point to the location of the scripts. Below are descriptions of some calculations.
- Please see reference National Bureau of Standards about the Fermi-Kurie Transform analysis.
You should read and understand about the Fermi-Kurie-Transform analysis. The reprint National Standards will give you all the needed information to understand this operation. See the faculty for more information.
This command takes two data sets, assumed to be beta spectra, and shifts them toward each other in order to account for hysteresis effects in the 111-Lab beta-spectrometer. You will be prompted for the 'UP' and 'DOWN' spectrum, at which time you press the alt-key numbers that correspond to the proper data sets. Then you will be asked for the amount to shift. The spectra are then shifted toward each other by means of x-axis rescaling (see that section of this manual for more info). That is, the x-axis is rescaled in the same way as the Calibrate x-axis low and Calibrate x-axis high commands. The shift command is cumulative, if you shift 50 bins and then shift -4 more, you will have a total shift of 46 bins. And the effects of the shift command may be removed with the Reset x-axis scaling command applied to each of the two spectra.
Compensate DataSet creates a new DataSet from the current one by dividing the y-value of each point by its calibrated (scaled) x-value, i.e.
This is useful for processing raw momentum data from the semicircular spectrometer, in which the bin-width is proportional to the bin.
Fermi-Kurie Transform creates a new DataSet that is the Fermi-Kurie plot of the current one [see NBS for complete treatment]. The transform is a theoretical one that changes the DataSet into a straight line. The transformation is:
where is the Quantum-mechanical correction for coulomb effects on the beta-particle as it leaves the parent atom. is calculated using the Bethe-Bacher approximation but the screening effect [ref 3 §12] is neglected (to see why, you should calculate the size of this effect for 137Cs). This transformation assumes that the x-axis is momentum and has been calibrated in units of mec. The resulting x-axis is then kinetic energy in units of mec2. If the input DataSet begins with negative momenta, Matlab will handle this information and begin processing data.
Inverse Fermi-Kurie Transform
The inverse Fermi-Kurie feature inverts the previous algorithm producing:
Beta Ray Data Analysis
What does artifact mean?
If you compare the spectrum taken in this lab to ones taken elsewhere, you will find in your spectrum, an extra peak on the lower momentum side of the k-peak (see figure 2.c "artifact"). Before we automated the scan process with a computer, the current in the magnet was adjusted by hand, and it was impossible to obtain the resolution necessary to see this anomaly. But now, with the advent of the computer, intrinsic errors in the beta ray spectrometer have become apparent.
To understand this, consider a source of finite width s emitting a monochromatic beam of β-rays subject to a uniform transverse magnetic field. Rays (electrons) are emitted in all directions. However, due to the Lorenz force they are confined to travel along arcs of the same radius r. Each electron travels along identical arcs, but not the same distance. Thus a detector of width Δx "sees"
an intensity distribution similar to that of Figure 6b. This implies that in the experiment, when a polychromatic beam is emitted and one "sits" at a bin and counts electrons, the count is smaller than it should be, because some electrons are falling into neighboring bins, due to the width of the intensity distribution. In essence our raw spectrum is the result of what we should see (the true spectrum) with this distribution (the error function) smeared into it. Mathematically this is known as convolution. We would write
G = T * E
So can we mathematically deconvolute our data and find the true spectrum? Yes in some special cases, but not here because the error function is not a constant width as we vary B and take our data. Also, the data are not accurate enough to warrant a detailed mathematical treatment. Instead, we smooth the curve by hand and eye in the region where there is an extra bump on the side of the peak, and proceed with the analysis.
So now we can deconvolute our data and find our true Beta spectrum Figure 6D. Resultant Beta spectrum, see the K-peak at Bin \# 2200 and L-peak to the right of it.
Data Analysis and Fermi-Kurie Plots
This section use the Matlab scripts The Fermi-Kurie plot is simply an algorithm that transforms your spectrum to another view of the energy spectrum. It uses the formula [B10] in the National Bureau of Standards reference. The expectation of this experiment is to experimentally determine Q1 and Q2 , respectively the high and low energy 'X' intercept points. FK (Fermi-Kurie ) transform makes it easier to determine Q1 and Q2 . First you must start with your full data spectrum from channel 0 to 4095 summed up and down data shifted together and then compensated, calibrated, and background-subtracted. Make a Folder in your My Documents folder for your data saving. Save your data ending in \*.DAT. The Excel files need to be opened and re-saved as "Text delimited with Tab" files ending in \*.DAT. After you have acquired your data, using the summed up and summed down files save them as the Fileup1.dat and filedn1.dat in excel as "Text delimited with Tab" files. Shift them together by the your calculated number, (see Matalab scripts) if they look okay then go on with the analysis.
If they are not perfect with respect to the hysteresis effects, then you will have to adjust the data yourself. You will need to shift each Y value by the different between the Y's at each X. In some instances you will have to average the Y's.
Then take your data and subtract background, Watch for negative numbers adjust if needed. Calibrate your X-axis, low and high points for momenta and compensate it "Q" in Matlab.
Save this dataset file as a new dataset in Matlab. You will want to save the files as you go to each step. Good practice in case of computer crashes, etc. Note that you can save files in ASCII or Binary format.
Now you are ready for generating the Fermi-Kurie plot;
Crop out unwanted bad data points like negative numbers, etc with the use of markers. The Fermi-Kurie plot that you should get is something like the below graph.
Fit a line to this curve and get the end-point on the X-axis, Q1. Pick a region above the noise and above the K-peak where it is linear. Now do a inverse Fermi-Kurie plot. You will get something like the following graph.
You should repeat the steps above to get the low-energy component and the low-energy end-point Q2.
Your result should be the low-energy component as seen below.
Your combined graphs should be like the following;
Reduce your data by following the numbered steps listed below, but first read Ref. 1 for detailed explanations of how beta-ray spectra are analyzed. Also read about data analysis.
In your Fermi-Kurie analysis, use units of and for momentum and energy respectively.
As you work on the data, keep in mind what your ultimate goals are: To observe the beta decay spectrum; to model the phenomenon with a Fermi-Kurie plot that incorporates the theory of beta decay; to determine the maximum energy available in the decay process. The above Figure shows four separate phenomena. Your data are the sum of these. We would like to separate this sum into the individual curves shown above, where they are shown overlapping but not added, as they are in the experimental data. See Appendix D about the Fermi-Kurie plot
- From the masses of 137Ba and 137Cs, calculate the total energy available to the electron mass energy + electron kinetic energy + neutrino energy in the beta decay. What is the maximum kinetic energy with which the 137Ba can recoil? Must this kinetic energy be taken into account in your experiment?
- Use Matlab for combining your MYDATUP1.DAT and MYDATDN1.DAT files to reduce the effects of hysteresis and to yield a single data file which you will use in all your analyses. Should you add the two sweeps? Should you displace the two spectra by a constant amount or a variable amount? Should you shift them right or left an amount that depends on channel number? Some thought and justification are needed here. Plot your two raw spectra and your combined spectrum.
- The momentum resolution Δp is proportional to the momentum p itself (you should derive this), and therefore, Δp is proportional to the current, since I is proportional to p. Divide the number of counts in each channel of your combined spectrum by I (or a number proportional to I). This is your compensated momentum spectrum Nc. Plot Nc vs. η.
- Calibrate your momentum scale in units of η by using the accepted value of the k-peak momentum (, or ), the bin in which the k-peak appears, and the bin number that corresponds to zero B-field (and hence zero momentum). Without accounting or correcting for hysteresis effects, how large an error in momentum would you have unwittingly obtained for the k-conversion electron momentum?
- Use the region of your spectrum that lies above the highest-momentum beta-particle to estimate the background radiation (you may assume the background is linear) and subtract it from your data. Plot the resultant spectrum.
- Make a Fermi-Kurie plot of your compensated spectrum. If the data were perfect, and the theory exact, then the F-K plot would be a straight line. Using a least-squares-fit line, determine the end-point of the high-energy beta component. How would the end-point plot change if the neutrino did not have zero rest mass? This need not be a mathematical explanation and Appendix F, for comments about error analysis.
- By taking the inverse-Fermi-Kurie plot of the best-fit line obtained above you will have the contribution of the high-energy beta-component to your compensated spectrum (Nc). Subtract the higher energy component from your compensated spectrum to leave only the lower-energy component.
- Make a Fermi-Kurie plot of the subtracted spectrum in the vicinity of the end-point of the lower-energy beta spectrum. Determine the end-point of the low-energy component.
- Compare the measured line width and the expected spectrometer resolution for the conversion electron peak.
- Check out Error Analysis Notes and the video on Error Analysis.
Copyright 2010 The Regents of the University of California. All rights reserved.
- M. Siegbahn, Beta and Gamma Spectroscopy, Ch. 3.‡
- K. S. Krane,“Beta Decay.” Introductory Nuclear Physics, Chapter 9, pp.272-288. John Wiles and Sons., 1987.
- Photo-multiplier tube Handbook‡ (Note this is a separate reprint booklet)
- I. Feister, "Numerical Evaluation of the Fermi Beta-Distribution Function", Physical Review 78, 4 (1950) ‡\*.
- L.R.B. Elton, "Chapter 9: Beta Ray Theory", Introductory Nuclear Theory, 2nd ed.; 1966, W.B. Saunders Company,PA;
‡Contained information on the Fermi-Kurie plots and calculations.
The reprints and equipment manuals are all available on-line from the Physics 111-Lab Library Site.