THE UNIVERSITY OF MICHIGAN INDUSTRY PROGRAM OF THE COLLEGE OF ENGINEERING NUCLEAR DECAY SCHEME ANALYSIS AND CHARACTERIZATION STUDIES OF (.d,alpha) REACTION PRODUCTS Donald Glenn Gardner ~ ~ ~. - ~ This dissertation was submitted in partial fulfillment of the requirements for the degree of Doctor of Philosophy in the University of Michigan, 1957. May, 1957 IP-2SO

W/17,A - r<s~t

Doctoral Committee1 Assistant Professor W. Wayne Meinke, Chairman -Professor Philip J. Elving Associate Professor Peter A. S. Smith Assistant Professor Milton Tamres Associate Professor Marcellus L. Wiedenbeck ii

ACKNOWLEDGEMENTS The author wishes to express his sincere appreciation for the guidance and encouragement continuously given by Professor W. Wayne Meinke throughout the course of this work. The many helpful discussions with Professor M. L. Wiedenbeck are especially valued. The section of this work concerned with the Fourier analysis of radioactive decay curves has benefited from the advice of many persons. In particular, the author is indebted to Professor J. W. Carr III who first suggested that the Fourier approach be tried, and to Professor A. J. Perlis, of the Carnegie Institute of Technology, who aided in many details of the solution. The careful reading and criticism of this section by Dr. B. Galler is appreciated. It is a pleasure to acknowledge the assistance of Mr. Ronald W. Shideler in the design and construction of the photographic coincidence apparatus, and in the maintenance of the electronic equipment. Thanks are also due to Mrs. Rosemary S. Maddock who assisted in various phases of this work. The cooperation obtained from the crews of the University of Michigan cyclotronj the Argonne National Laboratory cyclotron, and the University of California (at Berkeley) cyclotron is appreciated. The author is indebted to the U. S. Atomic Energy Commission for partial financial support of this research, and to the Eastman Kodak Company for the receipt of a research fellowship. Both the enriched and radioactive isotopes used in this work were supplied by the Isotopes Division of the U. S. Atomic Energy Commission. Thanks are due to the Industry Program of the College of Engineering for the typing of this manuscript. Finally, the author is truly grateful to his wife, Jeanne, who made many valuable contributions to all phases of this work related to data analysis and computers, and who painstakingly helped in the preparation of the manuscript. iii

TABLE OF CONTENTS Page ACKNCWLEDGEMENTS............................................ iii LIST OF TABLES..................................... viii LIST OF FIGURES................................................. x ABSTRACT................................................. xv CHAPTER I INTRODUCTION.................................... 1 II NUCLEAR THEORY.................................. 9 A. Nuclear Shell Structure..................... 9 1. The Shell Model of the Nucleus......... 9 2. Ground States.......................... 14 3. Excited States.......................... 18 -4. Shell Structure and Beta Decay................................... 21 B. Beta Decay................................. 21 1. Forbidden Transitions and Selection Rules......................... 25 2. Kurie Plots............................. 28 3. Comparative Half-life............... 30 III GAMMA SPECTRAL ANALYSIS......................... 32 A. Scintillation Coincidence Spectrometer...... 33 1. Components.............................. 33 a. Scintillation Detectors............. 33 b. High-Voltage Power Supply........... 36 c. Linear Amplifiers................. 38 d. Pulse-Height Analyzers.... 38 e. Count Rate Meter................... 39 f. Recorder.....*.................. 39 g. Oscilloscope........................ 4 h. Cameras..................*....... 40 iv

TABLE OF CONTENTS (CONT.) Page CHAPTER III (Cont.) i. Coincidence Analyzer.................. 40 j. Scaler................................ 41 2. Calibration of Spectrometer................ 41 B. Photographic Coincidence Method............... 42 IV BETA SPECTRAL ANALYSIS............................ 50 A. Description of Hollow Scintillator P-Ray Spectrometer............................ 54 1. Phosphors................................. 54 2. Operation................................. 58 3. Calibration.............................. 61 B. Performance of Hollow Scintillator P-Ray Spectrometer............................ 64 1. Effect of Detector Size................... 64 2. Effect of Sample Backing.................. 71 3. Backscattering from Sample Holder......... 73 4. Backscattering Shields and Collimators.... 73 5. Thick Samples............................ 75 6. Gamma Attenuation in Absorber............. 75 7. Amplifier Overload and Drift.............. 77 8. Range of Spectrometer..................... 78 9. Forbidden Spectra......................... 82 10. Complex Spectra.......................... 87 11. Positron Emitters......................... 90 12. Conversion Coefficients................... 94 13. Comparison with Flat Detectors........ 97 14. Correction for Instrument Resolution..... 99 15. Summary................................... 105 V DATA ANALYSIS: APPLICATION OF DIGITAL COMPUTERS.. 110 A. Kurie Analysis................................ 111 1. Computer Program for Kurie Computation.... 111 2. Simple Least-Squares Line Fit............ 113 B. Decay Curve Analysis.......................... 115 v

TABLE OF CONTENTS (CONT.) Page CHAPTER V (Cont.) 1. A Review of Previous Methods........... 117 a. Graphical Method.................... 117 b. Least-Squares Method................ 118 c. Prony's Method....................... 129 d. Householder's Method............... 131 e. Cornell's Method..................... 132 2. Method of Fourier Transforms............. 133 a. Development of Initial Problem into the Form of a Laplace Integral Equation..................... 134 b. Solution of Integral Equations by Method of Transforms.............. 136 c. Solution of the Laplace Integral Equation..................... 138 d. Details on the Method of Solution............................. 140 e. The Numerical Integration Scheme..... 150 f. Details on the Numerical Solution.... 153 g. Summary and Conclusions.............. 156 VI NUCLEAR DECAY SCHEME STUDIES................ 159 A. The Isotope Co2........................ 159 1. Chemical Separation.................... 159 2. Experimental Methods..................... 166 3. Results................................ 166 4. Conclusions.............................. 171 B. The Isotope Y92.......................... 174 1. Chemical Separation...........~........ 177 2. Experimental Methods............... 183 3. Results............................... 183 4. Conclusions.1................. 191 vi

TABLE OF CONTENTS (CONT.) Page CHAPTER VI (Cont.) C. The Isotope Ir96.................... 194 1. Chemical Separation...................... 194 2. Experimental Procedure................... 199 3. Results Obtained from Natural Platinum Bombardments................... 202 4. Results Obtained from Enriched Platinum Bombardments.................... 207 5. Conclusions.............................. 209 D. The Isotope Sc47........................... 211 1. Chemical Separation...................... 212 2. Bombardment Results...................... 216 APPENDIX A Fabrication of Hollow Plastic Scintillator.......... 219 APPENDIX B Energy Uhits.................... 222 BIBLIOGRAPHY.................................................... 223 vii

LIST OF TABLES Table Page I Distribution of Stable Nuclides........................ 10 II Possible Total Angular Momenta in the Configuration (j)K, Where K Equals the Number of Nucleons in the State.................................. 16 III Selection Rules in Beta Decay.................... 26 IV Unique Shape Factors............................. 27 V Log (ft) Ranges for Beta Transitions.................. 31 VI Gamma-Ray Calibration Source....................... 41 VII Characteristics of Some Scintillators.................. 54 VIII Response of Cs137 Conversion Peak to Pulse Height Analyzer Window Setting......................... 59 IX Conversion Electron Peaks Used in Calibration.......... 61 X Resolution of the Three Phosphors for the Cs137 Conversion Electron Peak............................... 68 XI Low-Energy Response of Some Phosphors.................. 82 XII Possible Multipolarities for Gamma Radiations.......... 95 XIII Comparison of Flat and Hollow Detectors............... 97 XIV Radiation Emitted During the Decay of Co62 and Ni62...................................... 162 XV Analysis of Enriched Nickel Isotopes................... 163 XVI Cobalt Separation Chemistry............................ 165 XVII Radiations Emitted During the Decay of y92 and Zr92.. 175 XVIII Chemical Impurities in Typical Zirconium Target Foil as Determined by Spectrographic Analysis.......... 179 viii

LIST OF TABLES (CONT.) Table Page XIX Yttrium Separation Chemistry........................ 181 XX Coincidence Results for Isotope 92........ 190 XXI Spectrographic Analysis of Typical Platinum Foil.................................................. 196 XXII Iridium Separation Chemistry...............19..........7 XXIII Cross-sections for the (d, a) Reaction on Platinum........................................... 204 XXIV Beta-Ray Spectrometer Data from Natural Platinum Bombardments................................. 204 XXV Gamma-Rays from 8.3-Day Iridium Activity.............. 205 XXVI Analysis of Enriched Platinum Isotopes................ 207 XXVII Gamma-Ray Energies in Mev from Ir89 and Ir... 210 XXVIII Radiations Emitted During the Decay of Sc4..... 212 XXIX Scandium Separation Chemistry......................... 214 XXX Gamma Rays from Scandium Fraction...................... 217 ix

LIST OF FIGURES Figure Page 1 Decay Scheme of Co6.......................... 2 2 Nuclear Level Sequence Proposed by Mayer.........13 3 Shape of Momentum Distribution for Beta Particles..... 23 4 Block Diagram of Coincidence Spectrometer.......... 34 5 Scintillation Coincidence Spectrometer............... 35 6 Photomultiplier and Detector in Lead Housing......... 35 7 Scintillation Detector Preamplifier and Photomultiplier Voltage Divider......................... 37 8 Gamma-ray Spectra of Cs 37..... 43 9 Gamma-ray Spectrum of Co.................44 10 Block Diagram of the Photographic Coincidence Arrangement............................... 46 11 Circuit Diagram of the Multivibrator Unblanking Gate................................................. 47 Gate.47 12 Gamma-ray Spectra of Co6............................. 49 A. Single Spectrum B. Spectrum in Coincidence with 1.17-Mev Gamma Ray C. Spectrum in Coincidence with 1.33-Mev Gamma Ray 13 Cross-sections of the Three Hollow Plastic Scintillators........................................ 56 14 The Three Hollow Plastic Detectors.................57 15 Phototube Holder and Lead Shield for Hollow Scintillator Beta-ray Spectrometer................... 57 x

LIST OF FIGURES (CONT.) Figure Page 16 Cross-section of Beta-ray Detector Mounting........... 58 17 Effect of Window Width on Resolution of Cs137 Conversion Peak........................................ o 18 Rate of Energy Loss as a Function of Energy for Electrons in Aluminum and Air......................... 63 19 Gamma-ray Backgrounds in the Three Plastic Detectors.. 6 20 Kurie Plots of Cs137 from Each Detector................ 67 21 Kurie Plots of Pm147 from the l.5-Mev and the 2.2Mev Detectors.................. 69 22 Kurie Plots of p32 from the 2.2-Mev and the 3.6-Mev Detectors...........................70 23 Effect of Sample Backing on Shape of the Kurie Plot of Cs137.......................................... 72 24 Effect of Backscattering from Sample Holder on the Kurie Plot of p32.....................74 25 Effect of Sample Thickness on the Kurie Plot of Cs137.............................. 76 26 Effect of Amplifier Overload and Drift on the Kurie Plot of p32.....................................79 27 Kurie Plot of S35...................,... 80 28 Beta Spectrum and Kurie Plot of Inl.............. 81 29 Kurie Plot of Cs137 Showing Effect of Forbidden Correction Term.83 Correction Term..................................... 83 30 Kurie Plot of Rb86 Showing Effect of Forbidden Correction Term................................. 84 31 Kurie Plot of Y90 Showing Effect of Forbidden Correction Term.............85.................... 85 xi

LIST OF FIGURES (CONT.) Figure Page 32 Kurie Plot of y91 Showing Effect of Forbidden Correction Term................................. 86 33 Kurie Plot of Cs137 + pm47 Mixture...... 88 34 Kurie Plot of Pm147 Resolved from Cs137 +Pm147 Mixture......................................... 89 35 Kurie Plot of Rb8..............91 36 Three-Channel Coincidence Arrangement for Positron Emitters...................................... 93 37 Comparison of Flat and Hollow Detectors........98 38 Kurie Plots for P32 from Both Flat and Hollow Detectors........... 100 39 Effect of Resolution Correction Near the Endpoint of Beta Spectra........................ 102 40 Resolution Correction by Bracewell....................103 41 Form of Experimental Line Shape Determined by Freedmen, et al........104 42 Experimental Conversion Peak Profile for Cs137... 106 43 Effect of Non-Gaussian Resolution Correction on Cs137 Kurie Plot................................ 107 44 Effect of Non-Gaussian Resolution Correction on y91 Kurie Plot................................. 108 45 Flow Diagram of Kurie Plot Calculation by the IBM 650 Computer............................... 114 46 Two-Component Decay Curves........................ 123 47 Sketch of f(x) as a Function of x.................... 14 xii

LIST OF FIGURES (CONT.) Figure Page 48 Sketch of g(X) as a Function of X.................... 143 49 Sketch of g(X) Showing Effect of Cut-off Error........ 145 50 Key to the Charts of the Isotopes............. 16 51 Chart of Isotopes of Cobalt, Nickel and Copper....... 161 52 Decay Curve of Cobalt Fraction....................... 167 62 53 Gamma-ray Spectrum of Co62 168 54 Kurie Plot of Co62........................ 170 55 Decay Scheme of Co62...... 172 56 Decay Scheme for y92 _ Zr92 - Nb92 as Given by K. Way,.et al............ 176 57 Chart of the Isotopes of Yttrium, Zirconium, and Niobium.................... 17 58 Decay Curve of Y-~ with Least-Squares Line Fit..... 84 59 Gamma-ray Spectrum of Y92 185 60 A. Total y92 Gamma-ray Spectrum B. Gamma-ray Spectrum in Coincidence with 0.47-Mev Gamma ray C. Gamma-ray Spectrum in Coincidence with 0.21-Mev Gamma ray........................ 186 61 Beta-ray Spectrum of 92............................ 188 62 Kurie Plot of y92 Ordinate Scale for Top Curve Differs from Scale for Lower Curves............... 189 63 Decay Scheme of y92 - Zr92 - Nb92.................... 192 64 Chart of the Isotopes of Iridium, Platinum and Gold.. 195 xiii

LIST OF FIGURES (CONT.) Figure Page 65 Foil Arrangement for the Three Natural Platinum Bombardments........................... 201 66 Beta-ray Spectrum of Iridium Activities Eight Days After Bombardment.......................... 206 67 Relation Between Mev and Relativistic Units........... 222 xiv

ABSTRACT The objectives of this work were to develop the necessary instrumentation, chemical procedures and aids in data analysis so that nuclear characterization studies of certain (d, a) reaction products could be made. A study was made to assess the value of a fluorescent plastic detector for use in P-ray spectroscopy. This led to the design, construction and evaluation of a hollow-type scintillation P-ray spectrometer. The performance of the device was evaluated with respect to such aspects as back-scattering from sample backing and sample holder, thick samples, geometry of detector, complex and forbidden spectra, and the useful energy-range. A comparison was made with the conventional flat-type scintillator. When corrected for instrument resolution, the hollow-type spectrometer enables 3-spectral distributions to be accurately studied down to energies of less than 0.1 Mev. Other instrumentation improvements include a photographic method of obtaining coincidence information. Existing methods for the analysis of complex radioactive decay curves were studied, and a resolution program was written for the IBM 650 digital computer using the method of weighted least squares. A new method of solution, based on a Fourier analysis, is proposed which seems to have advantages over all previous approaches. Here prior knowledge of the number of components in the decay curve is not needed since it automatically appears in the solution, and the calculation is not endangered by the presence of half-lives of similar value. Regarding n-spectral analysis, xv

one computer program was written to provide a Kurie analysis, and a second for general linear curve fitting using the method of unweighted least squares. Characterization studies were made of the isotopes Co62, y92, Ir196, and Sc47 which were produced by the (d, a) reaction in cyclotron-bombarded targets. Isotopically enriched samples were used with the Co62 and Ir196 studies. Chemical separations were required in all cases, and carrier-free procedures were developed for cobalt, yttrium, and scandium. New decay schemes are proposed for Co6 and Y92, and are substantiated by arguments based on nuclear shell structure theory. The previous assignment of an - 8-day activity of iridium to Irl96 is shown to be incorrect, the activity being due rather to Ir189 and Ir190. New 7-ray information is given for Ir89 and Irl90, together with excitation functions for the (d, a) reaction on t194 and pt196 for deuteron energies in the range of 9.6 to 20.4 Mevo The 7-ray information obtained from the scandium work verifies the data reported recently by other workers. xvi

CHAPTER I INTRODUCTION The subject fundamental to this work is the analysis of nuclear decay schemes. Basically, a decay scheme is a schematic representation of the decay process which includes the available information concerning the energy levels of the nuclei involved. In some respects, then, the nuclear decay scheme may be compared to the energy level diagrams associated with atomic or extra-nuclear excitation as in the case of the atom or molecule. The quantized excited states in an atom with orbital angular momentum and spin have their analogue in the nucleus, where protons and neutrons are believed to be in somewhat similar energy shells. Even what may be called "rotational" states occur in certain nuclei. As an atom gains energy, the number of available quantum states increases while the energy difference between the higher states decreases until finally ionizaticn is reached. A similar situation occurs in the nucleus, only here the "ionization" would correspond to the loss of a nuclear particle when its binding energy was exceeded. An isolated excited atom may lose its energy only by the process of photon emission when ionization cannot occur. It is here that the atomic and nuclear cases basically differ. An excited nucleus will also tend to lose its energy by the emission of photons or r rays. However, under certain conditions, the de-excitation may proceed by particle emission or capture, with the formation of a new element. This new element may or may not be formed in an excited -1

-2state. The above radioactive decay will proceed until eventually the ground state of a stable isotope is reached. The process of decay scheme analysis involves observing the radiation emitted during the decay and deducing properties of the energy levels such as spin and parity in the associated nuclei. In Figure 1 the decay scheme of Co60 is given, omitting some of the finer details. Co6 is observed to decay by the emission of a 0.3-Mev 3 ray to an excited state of Ni60 which then drops to the ground state by the emission of two 7 rays. The half-life of the P-ray transition is 5.2 years. After the number, energies, and relative intensities of the P-ray and 7-ray transitions have been determined, the data must then be assembled in the form of an energy Co60 5.2 YRS 0.31 MEV/X 1.17 MEV Y 2.50 MEV I 1.33 MEV I STABLE Ni60 Figure 1. Decay scheme of Co60

-3level diagram. One technique of value here is called "coincidence" analysis. Using suitable electronic equipment it is possible to determine which radiations produced in a given decay occur simultaneously or "in coincidence," within the resolving time of the equipment. In Figure 1 it would be observed that the P ray and both y rays were in coincidence. However, if the dotted y-ray transition could occur, it would only be in coincidence with the P ray. This is because once the second excited state of Ni60 has been reached, it may lose energy either by the emission of the 2.50-Mev y ray or else by the emission of the two lower energy 7 rays. In a given decay only one of the two paths may be followed. Coincidence information is of great value in the construction of decay schemes, but in general other information must be obtained before spin and parity assignments may be made. The term "spin" will be used here to mean total angular momentum. Some other types of information that are a measure of the spin and parity changes are: 1. The half-life and maximum energy of the P-ray transition. 2. The form or "shape" of the 3-spectral distribution. 3. The presence of meta-stable or isomeric states that may decay by the y-ray and/or 3-ray emission. 4. The presence of monoenergetic peaks, called conversion electron peaks, in the P-ray distribution. A great deal of the evidence (91) points to the possibility that the nucleus has, at least in low excitation states, a structure

-4similar to that of the electrons in an atom. It seems that certain numbers of nucleons produce a more stable nucleus, analogous to the completion of electronic shells. All but one of the isotopes studied here are within + 1 nucleon of one of these stable shells. This was done to facilitate the theoretical discussion of each nuclide. Furthermore, all but one of the isotopes contained an odd number of protons and an odd number of neutrons. These "odd-odd" nuclei have been in general less well characterized than other classes of nuclei, primarily because chemical separations are usually required. Finally, odd-odd nuclei are of interest because the arrangement of their energy levels differs on occasion from the arrangement found in other types of nuclei. The isotopes that were studied were all produced by the (d, a) reaction in the 42-inch cyclotron at the University of Michigan, the Argonne National Laboratory 60-inch cyclotron, or the 60-inch cyclotron located at the University of California at Berkeley. The product nuclei had therefore an a^"mic number of one less, and a mass number of two less than the target nuclei. The (d, a) reaction, at the energies used here, is a less favorable reaction by a factor of roughly 100 to 1000 than the competing (d, n) and (d, p) reactions. Hence it was necessary to devise a chemical separation procedure in all cases. The investigation was facilitated in two cases by the use of isotopically enriched targets. For example, Co was produced by the (d, a) reaction on a nickel target. The natural abundance of Ni64 is about 1%. With the enriched target this was raised to about 97%, which made possible the study of this isotope's reaction product.

-5Chemistry has become invaluable in present day nuclear scientific research. In many cases nuclear information cannot be obtained unless a prior chemical separation has been accomplished. The restrictions on the chemical separation are usually quite stringent. Frequently the separation must be rapid, as is obvious when half-lives in the order of minutes are involved. Furthermore, chemical separation factors may be required that are much greater than normally found in analytical chemistry. It is not uncommon to require separation factors of 106 or greater from all other elements present. It is here that such techniques as, for example, ion exchange (57, 76) and solvent extraction (29) are used to good advantage. Still further restrictions may be placed on the separation, such as the requirement that no inactive carrier may be used. This is quite common when a-spectral measurements are involved, since any weighable sample will distort the P-ray spectrum to a certain extent. In dealing with such minute amounts of materials as, say, 1010 atoms or less, new and interesting chemical problems arise that are not found in the macroscopic case. The contents of this thesis may be roughly divided into four general classifications, although this division is not adhered to in the text. Separations chemistry is one of these classificaticas. Carrier-free chemistries with good decontaminations were devised for the studies of Co62, 92 and Sc47. The chemical separation used to separate iridium from platinum and gold did not need to be carrier-free, since the iridium activity was found to decay by electron capture.

-6Another category is instrumentation. Much of the equipment used in this work, in particular the two-channel coincidence scintillation spectrometer, has been in use in the nuclear chemistry laboratories of the University of Michigan for several years (20). In many cases, however, it was necessary to modify the previous equipment or to design and construct new devices in order to cope with the specific problems found in this work. One of the new devices is a scintillation 3-ray spectrometer whose detector is a specially constructed hollow plastic scintillator. This spectrometer has been thoroughly investigated, and has replaced the older 1800 magnetic spectrometer (97) which now finds its best use in separating positron from negatron spectra. Several new scintillation detectors were obtained for 7-ray and x-ray work. As new and better phototubes have become available, these were incorporated into the equipment. A new automatic sweep device was constructed in the laboratory that will, in connection with the count-rate meter and Speedomax strip-chart recorder, continuously scan and record spectra over energy ranges governed by completely variable limits. A photographic coincidence technique was developed to handle short-lived activities. A third category involves data analysis. In the case of P counting, techniques have progressed far enough so that using 4n proportional counters, essentially every 3 particle down to energies of the order of 1 key or less may be detected. Very accurate decay curves are routinely obtained. And yet the half-lives, particularly in the case of complex decay curves, are usually measured graphically. Mention is made of some mathematical methods that have been tried

-7in the past for analysing decay curves, and a method making use of Fourier analysis is proposed. Due to the size of the problem this method has not been programmed for the IBM 650 computer. An alternative, although admittedly less elegant, method involving a weighted least-squares scheme has been programmed for the IBM 650. Two computer programs were written for use in the analysis of P-ray spectra. The first is a Kurie analysis program and the second is a simple least-squares program. The latter has been used to determine maximum P-ray energies from Kurie plots, but may also be used for any straight-line curve fitting problem. The Kurie analysis program was considered a necessity in the evaluation of the hollow scintillation 3-ray spectrometer since over sixty 3 spectra were analyzed during the evaluation. Several months would have been required had the analysis been done on a desk calculator, while once the data on a p-ray spectrum has been fed into the computer the computational time required for a thirty-point spectrum is less than onehalf minute. The last category is decay schemes. Using the techniques developed here the decay schemes of Co62 and y92 were determined, and are discussed from the point of view of nuclear shell structure. Spin and parity assignments were given to each energy level. An investigation of an activity previously assigned to Ir96 was made using both natural platinum and also platinum enriched in the isotope Pt198. The results showed that the activity did not belong to Ir96 but to other iridium isotopes. Excitation functions for the (d, a) reaction on two platinum isotopes were obtained. Finally, a preliminary investigation

-8of the decay of Sc was made, although the project was not completed since several excellent studies of this isotope were published while this wori was underway.

CHAPTER II NUCLEAR THEORY In the following sections, the subjects of nuclear shell structure and d decay are briefly discussed. An understanding of these subjects is vital in the analysis of nuclear decay schemes. A. Nuclear Shell Structure The development of the theory of nuclear shell structure must be considered to be only in a formative stage. It is not as yet possible to predict in advance energy levels and spins in any given nuclei. Actually, more experimental evidence must be obtained before the present theory and its modifications may be improved. When used as a guide, the present theory will at least limit the possible choices as to nuclear states, and in some cases will allow unambiguous assignments of spins and parities to be made. 1. The Shell Model of the Nucleus Since the earliest days of nuclear physics, attempts have been made to detect regularity in the structure of the nucleus analogous to the construction principle and electron catalogue (102) for atomic electrons. The distribution of stable nuclei is shown in Table I with respect to even and odd proton, neutron, and total mass numbers (52). Table I suggests that like nucleons are paired off with opposite spins. This is further confirmed by the fact that according to all known nuclear spins, nuclides with odd mass numbers have halfintegral spins, and nuclides with even mass numbers have integral spins. In particular, nuclides with A and Z both even have zero spin. Binding -9

-10TABLE I. DISTRIBUTION OF STABLE NUCLIDES Mass number, A Even Even Odd Odd Proton number, Z Even Odd Even Odd Neutron number, N Even Odd Odd Even Cases 162 4 55 5 energy considerations also indicate paired nucleons. However, some obvious shape irregularities are found in the binding energy curve at certain proton and neutron numbers. If one goes further, it is possible to correlate additional facts by assuming that the paired nucleons are arranged in shells after the fashion of the planetary electrons. Neutral atoms with electronic configurations of 2, 8, 18, etc. are exceptionally stable, and closed electron shells are spoken of, Atoms with one more or one less electron than that needed to complete a shell are exceptionally active chemically. While in the nuclear case, the effects aren't so dramatic, nuclei with 2, 8, 20, 28, 50, 82, and 126 nucleons of the same kind are particularly stable. These numbers are termed "magic numbers". A great deal of empirical evidence has been brought forth, particularly by M. G. Mayer, to support the concept of "magic" numbers (40, 52, 88, 90, 91). Several systems have been proposed to predict in a logical way the magic numbers by a proper quantum mechanical treatment usually based on the one-body model. The one by Mayer (89, 90, 91) seems to be the most satisfactory, and will be described qualitatively here. Other treatments (39, 40, 102) including a modification of Mayer's scheme will not be discussed here.

-11The first assumption is that each nucleon moves in a spherically symmetrical field of force, which is independent of the exact instantaneous position of all of the other nucleons. Since proton and neutron magic numbers are the same, it may also be assumed that the nucleons inside the nucleus retain their individuality to a fair degree. The quantum numbers (n, X, j, mj) characterize the nuclear orbits, and the corresponding energy levels are filled in succession by neutrons and protons. Proton and neutron shells are separate and are filled independently. Proton and neutron levels with the same quantum numbers do not coincide since the protons are subject to Coulomb forces which shift their levels to higher values. This Coulomb effect may even change the order of the proton levels slightly from that of the neutrons, because the lower I orbits penetrate more deeply and are consequently less favored for the protons. The nuclear potential differs considerably from the atomic force field, and the nuclear levels have very little correspondence with their atomic analogues. The actual potential lies somewhere between a square-well and a harmonic-oscillator potential. A rough estimate of the level order may be obtained by interpolating between these limits. Since the character of the nuclear forces is essentially unknown, and the level order is greatly dependent upon the potential chosen, there is little use in making exact calculations at this stage. What is done is to try to determine the exact order of the levels from evidence that lead to spin and parity values. Since the parity of the eigenfunctions is (-1)A it follows that all eigenfunctions belonging to one shell will have the same parity.

-12The above approximations show that up to N or Z = 20 a wide level spacing occurs beyond the occupation numbers 2, 8, and 20. However, there is no indication of gaps at 28, 50, 82, and 126. Mayer has resolved this difficulty by considering the spin-orbit coupling. In the atomic case, this coupling causes a splitting of the levels j = I + 1/2. This splitting is relatively small compared to the distance between levels of different A. In the nuclear case, however, this coupling has been shown to be quite pronounced (56), and the splitting considerably larger than for electrons. The energy splitting is proportional to (2A + 1), and also depends upon the radial quantum number n. Proton scattering experiments show that the level j = A + 1/2 has a lower energy than the j = A - 1/2 level, which is reversed in the atomic case. Since the level splitting decreases with increasing n, this means that the split will be the greatest when an angular momentum A appears for the first time. This is precisely where the higher "magic" numbers occur. Because of the fact that the energy of a given level depends strongly cn the way the spin and orbit are aligned, a level must be characterized by its J value as well as its A value. The levels will be filled for both neutrons and protoas such that the 2J + 1 states of j = ~ + 1/2 will be filled before the 21 states with J = X - 1/2. The level sequence proposed by Mayer, with strong spinorbit coupling is shewn in Fig. 2. On the left are the levels first found by extrapolating between the square well and the harmonic oscillator. These are connected by dotted lines to the resultant pairs of split levels. The degeneracy of the level is given in the parentheses,

-13I"-lj2 (1 6)-[184] —-184 — s,-.3d 23/2 (4)~-48 ~-'~ P —48k2~ (2)- 3d -', 2g/2 (8)6 hw I' 12 -(12)even - ~' 3d 5~ (6)* /' — 2g, ~ (10)-....,. li[i2 —(14) -[126] 126 (2)3P -— 3pp (4)-2f _____~2f (6)5 hi - 2fS, (8)-[100] odd / - lh~2 (10)~ ~\ -_____ _ h1, ~ (12)-[82] 82 3 —' s 3s~2 (2)2d0i (4)4 i- J 2d5/2 (6)-[64] even / 1 g2 (8)-.... g19- - (10)- [50] --- 50,- 2p' (2)- [40] 3hw r ~2p -" ~. lf. (6)-[38] 3 h ^ __. { - 2P3/2 (4)odd If. < lf~ f/2 (8)-[28] 28 2hw -2s~' ld 32 (4) - [20] 20 even - -ld -- 28 (2)- [16] 2'- ~ ~- ld~5 - (6)- -[143 1hw 1 ___ — lp - (2)- [8] 8 odd — lp2 (4)-[6] 0 - ls - ls (2)-[2] 2 Figure 2. Nuclear Level Sequence Proposed by Mayer

-14while the total occupation number of all levels up to and including the level under consideration is given in square brackets. It is seen that the nuclear shells are not the same as the oscillator shells because of the coupling. The precise order of the levels cannot be predicted without knowing the precise form of the potential well and the magnitude of the spin-orbit coupling. Furthermore, the order apparently differs slightly for protons and neutrons. The magnitude of the splitting, for the lg and lh levels at least, appears to be about 2 Mev. Essentially then, Mayer (90) proposes that: (1) the succession of energies of single particle orbits is approximately that of a square well with strong spin-orbit coupling giving rise to inverted doublets, (2) for a given ~ the levels P = I + 1/2 have invariably lower energy and will be filled before j = I - 1/2 levels, (3) pairs of spin levels within one shell, which arise from adjacent orbital levels in the square well in such a way that spin-orbit coupling tends to bring their energies closer together can, and very often will, cross, (4) an even number of identical nucleons in any orbit will always couple to give a spin of zero and hence will not contribute to the magnetic moment, and (5) an odd number of identical nucleons in a state j will usually couple to give a total spin of j and a magnetic moment equal to that of a single particle in that state. Very rarely, the total j will be one less than the single particle J. 2. Ground States In dealing with nuclear ground states the simplest case, other than when both neutron and proton numbers are even, occurs when

-15 - there is either one nucleon more or less than a closed shell in odd-A nuclei. Here the spin is due to the odd particle or hole and is just the j value for the level. When there are several neutrons and protons in incompletely filled shells the situation becomes more complex. In order to describe nuclear configurations the terminology.of Mayer (91) will be used. The number of neutrons and protons in a level will be denoted by K and X, respectively. The levels will be designated by the radial quantum number n, orbital angular momentum 1, and total angular momentum j, (nZj). n will be used for protons and v for neutrons. For example, (ld5/2) is the first neutron level with orbital angular momentum & = 2 and parallel intrinsic spin. If there are K neutrons in this state and they couple to give a resultant spin of I, K and will be used as superscript and subscript. Suppose three neutrons are in the ld5/2 level, and the resultant spin is 5/2, the configuration of the state would be (vld5/2)5/ Table II (91) lists the allowed levels of total angular momentum J for populations up to 2j + 1. For larger populations, the spin states are the same as for (2j + 1) - K(X). If there are X protons as well as v neutrons in incompletely filled levels, all vector sums are possible: IJV - JJrl J Jv + JJt. Mayer (91) states the following rules for ground state coupling: 1. The ground states of all nuclei with an even number of protons and of neutrons have zero angular momentum and even parity.

-16TABLE II. POSSIBLE TOTAL ANGULAR MOMENTA IN THE CONFIGURATION (j)K WHERE K EQUALS THE NUMBER OF NUCLEONS IN THE STATE j = 3/2 K= 1 3/2 K=2 0, 2 j = 5/2 K = 5 /2 K=2 0, 2, 4 K = 3 3/2, 5/2, 9/2 J 7/2 K = 7/2 K = 2 0, 2, 4, 6 K = 3 3/2, 5/2, 7/2, 9/2, 11/2, 15/2 K = 4 0, 2(twice), 4(twice), 5, 6, 8 J = 9/2 K =1 11/2 K = 2 0, 2, 4,6, 8 K = 3 3/2, 5/2, 7/2, 9/2(twice), 11/2, 15/2, 17/2, 24/2 K = 4 O(twice), 2(twice), 3, 4(three times), 5, 6(three times), 7, 8, 9, 10, 12 K = 5 1/2, 3/2, 5/2(twice), 7/2(twice), 9/2(three times), 11/2 (twice), 13/2(twice), 15/2(twice), 17/2(twice), 19/2, 21/2, 25/2

-172. In nuclei consisting of an even number of one type of nucleon, the ground state properties are determined by the odd type nucleons alone. 3. In a nucleus of odd-A, the nucleon of the type which is present in odd number will usually couple so that the total angular momentum is that of the last partially filled orbit, j. Occasionally, J = -1 is encountered. Rule 1 implies that pairs of nucleons with zero angular momentum are energetically favored. This means that the binding energy of the nucleus is increased by an amount P for each such pair. It has been observed that this "pairing energy" P increases with increasing J. When closely adjacent levels with different J values occur, the sequence of filling of these levels may be greatly affected by the difference in pairing energies for the two levels. Odd-A nuclei have been of prime importance in evaluating the shell theory presented above. That such a simple picture should be surprisingly accurate is attested to by the following facts. There are 101 odd-A nuclei with N or Z less than 84, and with experimentally determined spins. Of these, 93 fit the theory perfectly, 4 have spins of 1 less than the orbit of the last particle, and 4 show minor irregularities in the level order. In the same Z or N range there are 89 entries of data on magnetic moments. In 86 cases the magnetic moments agree with prediction, 2 are ambiguous, and 1 is a definite contradiction. Thus, for odd-A nuclei in which the odd nucleon is less than the 84 of its kind, the ground state characteristics are very nicely predicted by the shell model.

-183. Excited States It is not in general possible to predict the energy of excited nuclear states from shell theory alone. In a great many cases, spin predictions must be considered only as tentative, particularly in the higher excited states. However, in certain cases, shell theory can be used to limit at least the possible spin choices for a given state, and occasionally to give definite assignments. When experimental evidence such as beta decay, angular correlation, and/or conversion coefficient information is available, then spin and parity assignments of a more concrete nature may be made. One correlation of the excited states of even-even nuclei is called the "interval rule." This rule states that for nuclei with easily deformable cores, the lowest excited states will be due to intrinsic nuclear distortions similar to rotating surface waves. These states will all have even parity and carry an angular momentum J = 2 in the first, J = 4 in the second, and J = 6 in the third excited states. The corresponding energies will be a function of the J values as follows: E o< J(J- + 1) For most nuclei where the cores are not easily deformable the singleparticle excitations and energy splitting for different couplings are smaller than the states of rotating surface distortions, and such states are hard to identify among the lower excitation states of eveneven nuclei. In systems of "mirror" nuclei, that is, in pairs of nuclei such as,1~ - 6Cc which differ only in an exchange of proton and

-19neutron numbers, the level systems are expected to be similar. This is because the forces between proton pairs and neutron pairs are the same (as a first approximation). When all of one type of nucleon are in closed levels and the other nucleons are in an incompletely filled level, the lower excited states will be due to the latter type of nucleons. The available spin states are listed in Table II. For odd-A nuclei the order of these states with respect to increasing energy cannot be obtained from the table. For even-even nuclei, the first excited state is almost always 2+, that is, J = 2 and even parity. In fact, out of 69 known cases there are only 3 exceptions to this rule. The second excited state is also 2+ in about 1/5 of the cases, 4+ in another third, and the remaining third have various spins and occasionally odd parity. If the excited states arise from the same configuration as the ground state they must have the same parity, namely even. When both the proton and the neutron levels are completely filled the excited states arise from the promotion of one or more nucleons to a higher level. In this case, states of odd parity can occur. In a paper by Scharff-Goldhaber (123) the energies of the first excited states of even-even nuclei are plotted against the number of neutrons or abscissa. General trends, such as low excitation energies when both proton and neutron numbers are in the middle between two magic numbers, may be obtained from this curve. Another fact concerning even-even nuclei with incompletely filled shells is that there is always at least one excited state with J = 2, both for protons and neutrons, and often more than one.

-20Odd-odd nuclei present a very complicated case. Here both the proton and neutron shells have spins which couple together. The parity of the state is even if the values in and jp are either both even or both odd. If one is odd and the other even, then the parity is odd. In excited states, jn and jp may change in accordance with Table II. The new j values may then couple vectorally iJn - jpl J i jn + jp. Nordheim (103) has formulated two rules for coupling in odd-odd nuclear ground states. These are: 1. If both proton and neutron are in levels in which intrinsic spin and orbital angular momentum are parallel, as in Jn = 2n + 1/2, jp = p + 1/2, or if both are in levels where Jn = n - 1/2 Jp = 8p - 1/2, then jn and jp tend to add but not necessarily to the highest value J = in + jp. 2. If one of the nucleons is in a level with intrinsic spin and orbital angular momentum parallel and the other in an antiparallel level, that is, Jn = In + 1/2, jp = p - 1/2, or in = In - 1/2, Jp = 1p + 1/2, then the total angular momentum of the ground state is the smallest possible one, J = in - pl. These rules are not as successful when applied to excited states, as they are in ground state predictions. Another complicating fact associated with excited states of odd-odd nuclei is that the order of adjacent levels may not be exactly the same as in odd-A nuclei. When there is one proton and one neutron hole in the same level, j (except j = 1/2), they will couple to give a J = 2j - 1 or 1 less than the maximum. In the case where the last nucleon of one type is an S 1/2 level one expects to find an excited

-21level close to the ground state. This is because the S 1/2 nucleon's spin can only be added or subtracted from the angular momentum of the configuration of the other nucleons, and it is generally assumed that the spin-dependent forces are only about 20% of the total forces. 4. Shell Structure and Beta Decay Information derived from the study of 3 spectra is of great value both in testing the theory of shell structure and also in gaining information about levels where shell predictions are ambiguous. Myer (91) lists 200 cases of B decay from odd-A nuclei. Only 2 of these transitions definitely do not agree with shell model predictions. In the complicated case of odd-odd nuclei little direct experimental evidence is available concerning spins and parities since all but four of these nuclei are unstable. Beta-spectral analysis, however, can give information concerning the spin and parity change that occurred during the 3 transitions. When this information is combined with the y-decay data of the even-even daughter, in many cases the odd-odd ground state properties may be deduced. A discussion of the theory of 3 decay is given in the following section. B. Beta Decay There are several types of radioactive decay which are classified as { processes. These are (1) negatron emission, (2) positron emission, and (3) orbital electron capture. Positron and negatron emission may be represented by ZPA,(Z 1)Ai + e~ + v (1)

-22and orbital electron capture by ZA + e -(Z -1)A + (2) where v represents the possibility of accompanying neutral particles. One of the first difficulties encountered in the theory of 3 decay was that unlike Ca and y rays the electrons seemed to be emitted in a continuum of energy states from zero to a maximum energy. This would indicate that the total energy of the system was not being conserved. This, along with the fact that linear momentum, angular momentum, and statistics (130) were also not being conserved, lead Pauli in 1927 (52) to postulate a new particle, the neutrino, with the following characteristics: (1) very small or zero rest mass, (2) spin of 1/2, () obeys Fermi statistics, and (4) no charge. It is then the simultaneous emission of the 3 particle and neutrino, both sharing the available energy, which produces the continuous B-ray spectra. In 1934 Fermi set forth his theory of 3 emission. This theory, set up in direct analogy to the Dirac theory of 7 emission, predicts the shape of the 3 spectrum and a relationship between the maximum energy and the half-life. In it, Fermi postulated an interaction between the initial and final nucleus and the electronneutrino field, analogous to the coupling between an atom and its electromagnetic field wherein photons are created. The electrons and neutrinos do not exist in the nucleus any more than photons are thought to exist in excited atoms. Both are created during the emission process. This type of interaction seems to apply not only to 3 decay but to all known types of processes involving fermions. The known fermions

-23are: the nucleon (N or P), the electron (e+), the neutrino (v), and the mu meson (+-). A derivation of the equation for the momentum distribution of a P spectrum will not be given here since it is available in several references (12, 75, 128). However, neglecting for the moment the Coulomb interaction between the emitted 3 particle and the residual nucleus, the momentum distribution of the P particles can be shown to be N(p)dp p p2(Eo - E)2dp (3) where Eo is the maximum energy of the transition. This expression qualitatively describes the shape of observed 3-ray spectra. A plot of this expression is given in Figure 3. N(p) Figure 3. Shape of momentum distribution for beta particles

-24For large momenta N(p)dp is proportional to (Eo- Eo)2, for small 2. momenta it is proportional to p. The distribution vanishes at both limits and has a maximum in the middle. So far we have ignored the Coulomb force between the 3 particle and the daughter nucleus. This force can only be neglected for the lightest nuclei, and high electron energies. The effect of the force is simply to accelerate the positive electrons and to decelerate the negative ones. Hence, the spectrum of a positron emitter will show relatively few slow particles, while the negatron spectrum shows more slow particles than it would in the absence of the Coulomb effect. An approximate expression for F(Z, E), the Coulomb correction factor is F(Z, E) t 2 (4) 1 - e-2 where =+ze2 and (5) fiv v = Velocity The upper sign in r is for negatrons, and the lower for positrons. This is given in Reference 12 and was obtained by treating the electon non -relativistically. The relativistic treatment is also given, but the simple form of the above expression shows more clearly how F(Z, E) varies with energy. With no Coulomb force the momentum distribution is proportional to p at low energies. For negatron emission with low momenta, F(Z, E) becomes inversely proportional to p, and the distribution is proportional

-25to p rather than p2. For positron emission F(Z, E) becomes proportional to 1 e-c/p (c is a constant), and the distribution is proportional to pe-C/P. It is seen then that the positron spectrum drops to zero much faster than the negatron spectrum. At high energies and low Z, where rT is small and changes slowly with energy, F(Z, E) becomes more like a constant and doesn't affect the shape of the spectrum as much. We now have the relationship N(p)dpoC F(Z, E)p2(E - E)2dp (6) 1. Forbidden Transitions and Selection Rules Consider a nuclear factor represented by the symbol [M2] This is a sum of terms, some of which are energy dependent and others are not. Whether [M2] is a function of energy depends on the magnitude of the various terms. Let the orbital angular momentum of the electron and neutrino be Be and jv, respectively. Now the wave function 4e for an electron emitted with a wave number ke = Pe/h and orbital angular momentum Be is proportional to (kere)e for small distances re from the center of the nucleus. It would appear then that re at the center of the nucleus would vanish except when le = 0. The same holds for the neutrino. If the above approximations were true, then only electrons and neutrinos with zero orbital angular momenta would be emitted. Actually, the particles may be emitted with le and lv different from zero; however, the probability for such emission is much smaller. Therefore, transitions where Se = v = 0 are called

-26"allowed" transitions. "Forbidden" transitions are those where one or both of the particles have ~ I 0. Both the neutrino and the electron have spins of 1/2. These may combine to give a total spin of 0 or 1. As the orbital angular momentum carried off during emission is zero for allowed transitions, the change in the total angular momentum At may be 0 or i 1. The parity of the system is determined by the sum of the individual I's and is even if L is even and odd if L is odd. It follows that for allowed transitions there is no change in parity. Selection rules for some P transitions are given below. TABLE III. SELECTION RULES IN BETA DECAY Transition Parity Change L = "Spin" Change Allowed no 0, 1 1st Forbidden yes 0, 1, 2 2nd Forbidden no 2, 3 3rd Forbidden yes 3, 4 For allowed transitions [M2] is approximately equal to one and independent of energy. Each successive degree of forbiddenness decreases [M2] by at least a factor of 100. In forbidden transitions [M2] is usually a function of energy, but many cases are known in which [M2] varies only slightly.

-27Thus, we are lead to the result: N(p)dp z F(Z, E)p2(Eo - E)2[M2]dp (7) The spectrometer to be described in a later section measures energy (E) rather than momentum (p), and so it is convenient to rewrite equation (7) in terms of energy. Instead of using E which is measured in terms of Mev (million electron volts), it is more convenient to use the term W in relativistic or rest mass units. See Appendix B. Here W= E + 1 (8) mec where mc2 is the rest mass of the electron. Therefore, N(W)dW CpW(Wo -W)2F(Z, W)[M2d1 (9) In general, forbidden transitions are difficult to study because of the complexity of the factor [M2]. There is, however, one class of forbidden spectra which may be analyzed rather easily. These are said to be "unique" forbidden spectra, and occur when the spin change is AJ = n + 1, where n is the order of forbiddenness. These spectra were theoretically predicted to be unique since only one nuclear matrix element fails to vanish. In this case [M2] can be represented by the shape factor Cn(W) where n is again the order of forbiddenness. The three shape factors for first, second, and third forbidden spectra are listed below. TABLE IV. UNIQUE SHAPE FACTORS. (W2 - -....... - - - (w) = (W2 - 1)+ (Wo -W)2 2(w) = (W2 - 1)2 + (Wo - w)4 + 10/3(W3 - l)(Wo - w)2 a(W) = (w2 - 1)3 + (Wo -w)6 + (W2 - l)(Wo - W)2(W2 - 1) + (W -W)2]

-282. Kurie Plots Referring to equation (9) it may be seen that in the case of allowed transitions where [M2] is close to unity and independent of energy, pWF(Z, W) If the observed energy spectrum is used to compute F N(W) 1 1/2 pWF(Z. W)j which is then plotted against W, an allowed spectrum will yield a straight line intersecting the abscissa at Wo. This plot is referred to as a Kurie plot after the author who suggested this method of analysis (80). Tables of pWF(Z, W), called the Fermi function, are available from the National Bureau of Standards (99). The Fermi function must be corrected for the effect of the outer atomic electrons, but this may be done using Table 6 of Reference 99. A necessary parameter is the "outer screening potential," Vo, which is available from Figure 13 of the same article. The effect of the screening correction is to reduce the value of the Fermi function for 3" particles, since the energy of the negative P particles which must penetrate the outer screen of electrons is not as high as it would have been in the absence of the screen. Conversely, a positron would have a higher kinetic energy. The beauty of the Kurie representation is that the endpoint value need not be known. Because the more accurately determinable points in the body of the spectrum may be used in a linear extrapolation to give the value at Emax, it is not necessary to experimentally determine points in the

-29low counting rate region near Emax. The plot is also a simple means of testing a spectrum to see if it has the allowed shape. The presence of more than one group of 3 rays can sometimes be recognized immediately from the shape of the spectrum, as in the case of Sr90 - Y90. More frequently it is only by means of a Kurie analysis that the presence of more than one group will be recognized. Since an allowed Kurie plot is linear, a straight-line extrapolation of the highest energy group can be made. This is subtracted from the total curve, just as one resolves half-lives in the case of complex decay curves, and a new Kurie plot made from the difference. Since the number of electrons in each group is proportional to the area under each straight-line curve, one can easily obtain the relative abundance of each 3-ray group. Because of the errors involved in the subtraction procedure, the shapes of the lower energy components may be inconclusive. More accurate information concerning the lower groups may sometimes be obtained using coincidence methods. Where forbidden transitions are involved, one must have information on the factor [M2]. If this can be obtained, as in the case of unique transitions, the corrected Kurie plot may be obtained by dividing the uncorrected plot by [M2]1/2. Kurie plots will all show an upward curvature at the low energy end due to source thickness, scattering from the backing material, etc. Still, if the spectrum does not contain lower energy groups there will be an appreciable straight section. If there is not, then this is a strong indication that the transition is forbidden. Merely because the Kurie plot is straight does not, however, mean that the transition

-30is allowed. For example, most first forbidden spectra show allowed shapes. In fact, until 1949, all 3 spectra investigated appeared to have an allowed shape except that of RaE. This is because the energy dependent terms are overshadowed by large energy-independent terms. All second and third forbidden spectra so far measured have shapes that differ from the allowed shape. If the Kurie plot is straight, the selection rules for allowed transititions fulfilled, and/or a function called log ft (to be discussed in the next section) is in the allowed range, one would conclude that the transition is allowed. 3. Comparative Half-life Another item of interest is called the "comparative halflife." If the equation for the energy distribution is integrated from zero to the maximum energy, one obtains the total probability per unit time that a P-active nucleus will decay. This probability defines the half-life. 3" decay probability WO X = ln2 = / N(W)dW a f_(Z, W) t1/2 1 where (11) W f (z, W) = f F(Z, W)pW(Wo -W)2[M2]dW 1 Thus, f-(Z, Wo) times tl2, or simply "ft," is a constant of the disintegration. The value of ft is characteristic of the type of transition, i.e. allowed transitions have a value of ft smaller than the 1st forbidden transitions, and so on. Since the values of ft for kncown transitions are large and tend to scatter, what is usually done is to compare log10(ft) values. It is found that there are no distinct separations

-31between the various orders of forbiddenness, only trends as is shown below in Table V. TABLE V. LOG( ft) RANGES FOR BETA TRANSITIONS Transition Log ft Supperallowed 3.0 - 3.7 Allowed 4.4 - 5.9 1st Forbidden, 7 2nd Forbidden 8.0 - 11.0 The "superallowed" transitions are a favored case, as in transitions between "mirror" nuclei. In spite of the fact that the ranges overlap, log ft values are of use in characterizing transitions if other information such as parity change or spectral shape are known. The values of f(Z, WO) have been extensively tabulated by Feenberg and Trigg (41).

CHAPTER III GAMMA SPECTRAL ANALYSIS When a y ray passes through a substance it may lose its energy by three processes: 1. Photo-electric effect. The entire energy of the y ray is given to an electron in the substance, causing ionization. 2. Compton scattering. Here the initial y ray gives only part of its energy to an electron, and the secondary, lower energy y ray is scattered from the point of collision. 3. Pair-production. When the initial y-ray energy exceeds 1.02 Mev the creation of a positron-negatron pair may take place. If the substance struck by the 7 rays is fluorescent, such as a T1 activated NaI crystal, then the light out-put is a measure of the 7-ray energy. It is most convenient to use the monoenergetic peak caused by the photoelectric effect, since here the light intensity is directly proportional to the 7-ray energy. In the Compton scattering case, the energy of the secondary y ray is a function of the energy of the primary y ray and the photon scattering angle. Therefore, the energy transferred to the substance will vary continuously from 0 to a maximum energy, known as the Compton edge, where the photon scattering angle is 180~. The resulting 7-ray spectrum will then appear as a continuous distribution from 0 to the Compton edge followed by the photo-electric peak. The relative intensities of y rays were determined by integrating under the respective photo-peaks. Due to the fact that the -32

-33absorption coefficient for the photo-electric process is dependent upon the energy of the incident y ray, the size of the detector, and atomic number of the atoms in the detector, a correction must be made on the measured relative intensities. McLaughlin and O'Kelly (92) discuss how the correction should be made. The application of the methods of coincidence counting to experiments in nuclear physics has been described by Dunworth (37), Wiedenbeck and Chu (147) and others (85, 98), and so this subject will not be discussed here. A. Scintillation Coincidence Spectrometer The dual channel scintillation coincidence spectrometer used in this work is -a slight modification of the apparatus originally built and described by W. A. Cassatt (20). An overall block diagram of the system is shown in Figure 4, while a photograph of the spectrometer appears in Figure 5. 1. Components a. Scintillation Detectors. Three sodium iodide, thalliumactivated scintillation detectors are available for use with y rays. The crystals are right circular cylinders 1-1/2" in diameter by 1" thick mounted in hermetically-sealed aluminum housings, using magnesium oxide as a light reflector. One of these detectors was purchased already assembled from Harshaw Chemical Company, Cleveland, Ohio. The other two had been purchased in rough form from Harshaw and polished and assembled by Cassatt. One of the latter two detectors was remounted by the author when the crystal lost its optical coupling with the glass window cemented

COINCIDENCE SPECTROMETER Phototube ~ Phototube and and Amplifier Preamp Amplifier AmlfrPreamp ^Preamp Regulated Power Supply Regulated Power Supply Single Channel CSoincidence ingle Channel Pulse Height A Pulse Height Analyzer Analyzer Tro ficounter Automatic Prnt Type Scaler Scaler Linear Sweep Recorder Control SF. 4. B k d o co aler Fig. 4. Block diagram of coincidence spectrometer.

35-.......... Di i-......................... l ni!!!!..............." l ~ i~`i'i;j ~ ~ ~ F;I~~: iii~~~~~~~~~~~~I1Eltlflfllr~~~~~~~~~~~~~~~~dt~~~~f~~~~'I — s 4~~~~~~~~~~~81B~~~~~~~i.~~~~~,%Li~~~~~~~~;h~~~~~.....................~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~'..........iii ii ii............................ ~ ~ ~ ii~ii~~j~~jC~~i........................... ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~~~ii~i~.................................~~~~~~~iii............................en e S ec rom te:::ix~ ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~v:.................................~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~. -i:;__i: —— _-ii::iiiliiiii-:ili~ii-liiilil-:::i-..........:::~ ~:::_:::::;-:::::-:::::i~i~i~i:........................~~~~~~~~~~~~~~~~~~~ii::'l:iiiii.":....:::~:: _: -:i —:-::-:: -:-:i-::-i::iiii~i:i:ii- _ii-.......................i~-: l~ii~iiii...............pir nd Deeco.......... ausin

-36to one end of the aluminum housing. For details on polishing and mounting crystals see Reference 20. One x-ray detector was obtained in assembled form from National Radiac, Inc., Newark, New Jersey. This consists of a thin slice of thallium-activated sodium-iodide crystal in the form of a right circular cylinder 1-1/2" in diameter by 2 mm thick. It is thin enough to be reasonably transparent to y rays and yet thick enough for use with x-rays. Several different 0-scintillation detectors are available. One is a 3/16" thick by 1-1/2" diameter trans-stilbene crystal. Others are flat cylinders of Sintilon-brand plastic scintillator manufactured by National Radiac. These are available in several diameters and thicknesses. One type of specially constructed plastic a-ray detector will be described in detail in a later section. A detector is connected to the spectrometer by optically coupling it to the flat photocathode of a Dumont 6292 multiplierphototube with Dow-Corning 200 oilicone fluid. A mu metal shield is then placed around the phototube to reduce the effects of magnetic fields, and the phototube is enclosed in a light tight housing which also contains a cathodefollower preamplifier. A lead enclosure houses two such detector-preamplifier units placed at 1800 to each other. Figure 6 shows the detector arrangement, while Figure 7 is the preamplifier schematic. Low voltage is supplied to the preamplifiers from the linear amplifiers. b. High Voltage Power Supply. High voltage is supplied to the preamplifiers by an Atomic Model 312 Superstable High Voltage Power

VOLTAGE DIVIDER AND PREAMPLIFIER CIRCUIT R16 51^~~,Icl~~~~?eI I H.V. 5 + 300v ^//<^^^^^4~~~10 10 DC~^^ ^^^^ Vi~~~~~~~~ ^^13 ^p.~.3v Rlt ~i~~ "r' ~o "CI 1 I R3 CL\X,~ ~ r ~R13 I - 24 I I R15 R1 5~~~~~~~~~~~~~~~~~~~~~~~~~~~~(A R = 2.2 Meg, 1/4 w C.001, 600 v V = Dumont 6292 1 R2 -R =1 Meg, 1/4 w C2.001, 600 v V = 6J6 R13 =1 Me g, 1/4 w 2 3 1 R = I Meg, 1/4 v C5=.05, 16oo00 v P = Amphenol AN 3102-145-25 R14= 4.7 K, I w C4 =.001, 3500 v MICA P2 = Amphenol 831R R15 = 30 K, 1/2w C5 =.05, 1600 v R16 = 1 Meg, 1/4 w Figure 7. Scintillation detector preamplifier and photomultiplier voltage divider

-38Supply. This replaces the power supply originally constructed for the spectrometer. The new power supply has increased the resolution of the spectrometer by reducing the ripple in the output, and shows much better resistance to variations in line voltage and long-time drift. In normal operation a potential of 830 volts is used. Since the gain of a photomultiplier is very sensitive to variations in voltage a Sorensen Model 2000S voltage regulator is used between the power supply and the 110 a.c. line. This particular regulator is a saturable reactor type, and gives a regulation accuracy of ~ 0.1%o over an input range of 95 to 130 volts. Most of the other units of the spectrometer also have their power supplied by the Sorensen. c. Linear Amplifiers. The linear amplifiers are both Atomic's Model 204B. The gain may be varied from 3,500 to 15,000 by means of a course and a fine gain control. The original variable controls were replaced by two 11-position step switches for the sake of better reproducibility. d. Pulse-Height Analyzers. The single-channel differential pulse-height analyzers are Atomic's Model 510. Ten turn heliopots are provided for setting the base-line descriminator from 0 to 100 volts, and for setting the channel width from 0 to 7 volts. At a window setting of more than 7 volts there is no upper limit to the window, and consequently all input pulses larger than the base line setting will produce an output pulse. The analyzers have been modified so that the base-line height can be varied either manually or by an automatic sweep mechanism. A new sweep drive has been constructed which will allow a spectrum to be swept

-39continuously, back and forth while a "slide-rule" dial indicates the position of the sweep at all times. Finally, the positions of the reversing switches are completely variable so that any given portion of the spectrum may be scanned at any of five different sweep speeds. In addition, a step sweep is provided which will advance the base line in discrete steps. A preset counter will advance the step sweep mechanism upon the accumulation.of a given number of counts from the singlechannel analyzer. e. Count-Rate Meter. The count-rate meter, used in connection with the linear sweep, receives its pulses from the pulse height analyzers. The output of the meter is fed directly into the chart recorder. A selection of eight linear scales ranging in full scale counting rates from 100 to 20,000 counts-per-minute, and four values of the charging condenser from 1 to 16 microfarad are available. A 7,200 counts-per-minute pulse generator is built in for calibration purposes, and a meter in the output circuit may be used for visual monitoring. f. Recorder. A Leeds and Northrup # 60352-M02-191-Spec Model S Speedomax, four point, printing chart recorder with a 10 millivolt full scale sensitivity is used. Other features include: multicolor printing, dots only; strip chart tear-off device without indices; three seconds per point with two second balancing speed; and a standard cell against which the instrument automatically calibrates itself periodically. The recorder is connected so that during a measurement every fourth dot printed indicates the baseline voltage. The other three dots record the output of the counting-rate meter.

-40g. Oscilloscope. The Tetronix scope Model 514D is used. The scope is provided with a blanking circuit which cuts off the trace except when a pulse is being viewed. This prevents an objectionable spot from appearing at the zero position of the trace between pulses. h. Cameras. Two cameras are available. The first is a DuMont Model 295 Polaroid Land Camera. With it oscilloscope pulse traces can be photographed directly. When suitably calibrated, the picture may be measured to determine y-ray energies (35). By means of a sweep expander and a neutral gray wedge filter (10, 22, 122) the relative intensities of the y rays may be determined. Due to the relatively low contrast of polaroid film, a DuMont 35 mm oscillograph camera is sometimes substituted for the Polaroid Land Camera. Kodak Linagraph Ortho film is used for these pictures. The Land Camera is very good for taking rapid survey pictures, since the finished positives are ready one minute after exposure. i. Coincidence Analyzer. The two twin single channels, each consisting of detector, preamplifier, linear amplifier, and pulse-height analyzer, are connected to an Atomic Model 502A coincidence analyzer. Actually, the instrument has four input channels. Any combination of the channels may be operated in coincidence, or any set of the other three channels may be used in anticoincidence with channel one. The resolving time of each channel may be set at 1/4, 1/2, 1, or 2 microseconds. In normal use coincidences are sought between certain pulses appearing in one detector and the entire spectrum, scanned in discrete steps, as it appears in the other detector. B-y or y-y coincidences are observed with the proper choice of detector and

-41absorbers. When pulses generated in each detector arrive at the coincidence analyzer within the resolving time of the instrument, an output "coincidence" pulse is obtained. j. Scaler. The gross output of each channel, along with the coincidence output, is fed to a specially constructed three-channel scaler. The scaler was built by Radiation Instrument Development Laboratory, Chicago, Illinois, especially to fit our need. Each scaler channel has a scale of 100, composed of Berkeley decimal counting units, and an electric reset mechanical register. Provision has also been made for making the pulses from the scale-of-100 decimal stage available at the rear of the scaler. 2. Calibration of Spectrometer The coincidence spectrometer was checked for linearity and calibrated in the usual manner using radioactive samples emitting y rays of known energy. Some of these sources and their energies are shown in Table VI. TABLE VI. GAMMA-RAY CALIBRATION SOURCES Source Principle 7-Ray Lines (Mev) I114 O. 9 in 0.190 Cr51 0.32 Sn113 O.394 Cs137 0.662 Znr65 1.12 Co6O 1.17, 1.33

-42Two typical 7-ray sweeps and pictures are shown in Figures 8 and 9 for Cs137 and Co6, respectively. The coincidence apparatus was checked using isotopes with known decay schemes such as Cs137 and Co60. The resolving times were checked by the method described by Bleuler and Goldsmith (13), and found to agree with the dial readings within experimental error. B. Photographic Coincidence Method When it is necessary to obtain coincidence data on shortlived materials, and when multichannel pulse-height analyzers are not available, some type of photographic technique is generally used. A point-by-point manual sweep of the coincidence spectrum would necessarily take many half-lives to complete if good statistical results wereto be achieved. A photographic method, on the other hand, is much faster since an entire spectrum may be viewed at one time for coincidences with the other channel, rather than only a narrow window width which sweeps the spectrum in discrete steps. Also, since all parts of the spectrum are viewed simultaneously, no correction for decay is necessary when relative peak intensities are to be measured. Several photographic methods have been described in the past (10, 18, 22, 33, 49, 67, 122). Some of these require motion picture cameras, while others require specially modified oscilliscopes. All require additional circuitry of various degrees of complexity to produce the desired results. The photographic method described here requires a minimum of additional circuitry, is simple to construct and operate, and is easily adaptable to Gray-Wedge analysis (see Section lBh on Cameras). The results are simple to interpret since the coincidence spectrum may be directly compared to the total or "single" spectrum.

-43GAMMA-RAY SPECTRUM OF Cs 137 24,000 LLU H 16,000 LU o 8,000 0 0.2 0.4 0.6 0.8 1.0 ENERGY (MEV) Figure 8. Gamma-ray spectra of Cs137

a — GAMMA- RAY SPECTRUM OF Co60 16,000 ) 12,000 4,000 ( 4 10 0 0.2 0.4 0.6 0.8 1.0 1.2 ENERGY (MEV) Figure 9. Gamma-ray spectrum of Co60

-45Figure 10 is a block diagram illustrating the coincidence set-up. The detecting phosphors are chosen for P-y or y-y coincidences. Channel II is the interdicting channel and is set on the pulses with which coincidences are sought. The base line of the channel I pulse height analyzer is set below the lowest energy point in the spectrum above which information is desired, and the window opened completely so as to view all of the spectrum above that point. Usually the baseline is set just above the noise level. Occasionally, when there is an intense y ray which is obscuring y rays immediately above it by overexposing the photographs, the intense y ray may be eliminated merely by setting the baseline just above it. Pulses from the channel I linear amplifier are fed to the normal oscilliscope input with internal trigger after first passing through a 4-L sec. delay line terminated by a reostat of 2500 ohms. When a coincidence occurs, the coincidence analyzer's output pulse triggers a monostable multivibrator unblanking gate whose output is an \,20 p sec. square wave. Figure 11 is the schematic of the gate circuit. The square wave pulse is applied across the oscilliscope's CRT cathode terminal and ground. With the brightness control on the oscilliscope turned all the way down, and the pulse height of the square wave adjusted to - 60 volts, a coincidence event will then illuminate the trace for \ 20 p sec. and reveal the pulse that gave rise to the coincidence. At least 106 pulses per minute from channel I may be applied to the scope before the chance occurence of two pulses appearing within a few p sec. of one another becomes objectionable.

CHANNEL I DETECTOR PULSE LI NEAR AND ~ ~ HEIGHT AMPLIFIER PREAMP ANALYZER GATE COI DELAY COINCIDECE OSCILLISCOPE PULSE LINE ANALYZER M.V. DETECTOR LINEAR PULSE AND ~ HEIGHT AMPLIFI ER PREAMP ANALYZER CHANNEL EL Figure 10. Block diagram of the photographic coincidence arrangement

,o B+ 300V. IM 27K IM I ri^^^^ ^~25K 2M^ / /,^. 2< OUTPUTISK 470 ] __ ^ ^ ^ ^^^ Hi ^v^ ^ ) ~~12AT7 _JA INPUT l|4 IOK 500 K Figure 11. Circuit diagram of the multivibrator unblanking gate

-48Still pictures are then taken of the ocsilliscope's presentation using either a Land polaroid camera or else a 35-mm camera. If desired, the pulses may be stretched and Gray-Wedge pictures taken. Since the gating pulse is quite flat on top the trace intensity remains constant and independent of when the coincidence pulse appears. The singles spectrum may be observed merely by switching off channel II on the coincidence analyzer, thereby producing an output pulse for each channel I pulse. Figure 12 shows the results obtained for Co6. Good results have also been obtained in B-y coincidence work. Gamma rays in coincidence with portions of the d spectrum have been measured, using a plastic phosphor as the 3 detector.

-496o Figure 12. Gamma-ray Spectra of Co. (A) Singles Spectrum. (B) Spectrum in Coincidence with 1.17-Mev Gamma Ray. (C) Spectrum in Coincidence with 1.33-4ev Gamma Ray

CHAPTER IV BETA SPECTRAL ANALYSIS In order to utilize the potentialities of the Kurie analysis previously discussed, the spectrum measured by a P-ray spectrometer must either be essentially free of distortion, or else a suitable procedure for correcting the spectrum must be available. Certain types of spectral distortion are common to most types of spectrometers. These include distortions due to thick P-ray sources and backscattering from the source mounting and surrounding equipment, which tend primarily to distort the lower energy section of the curves. These errors cannot in general be removed by subsequent data treatment. Another source of distortion is due to instrument resolution. The extent of this effect depends upon the type of spectrometer involved. Certain magnetic spectrometers are in use today which have excellent resolution, in the order of a few percent or less. While undoubtedly these spectrometers will provide one of the best spectral determinations available, they do in general present certain drawbacks. The excellent resolution is counteracted by the requirement of high sample activity, of the order of tens of millicuries, due to low transmission. In many instances the time required to pump the spectrometer down to vacuum after a sample has been placed in the system sets a definite limit on the half-life of the source. Further restrictions include the cost of the spectrometer, and the occasional difficulties encountered in constructing 3-7 coincidence equipment. The ability to separate positrons from negatrons is quite a definite advantage, however. -50

-51Another type of spectrometer, whose resolution is inherently poorer than the magnetic type, utilizes a fluorescent detector after the fashion of the scintillation 7-ray spectrometer. One method (63, 113) involves the use of a thin flat phosphor covered with a reflecting material such as aluminum foil. This arrangement, however, tends to distort the low-energy portion of the 3 spectrum since some of the B particles will be scattered out of the phosphor leaving only part of their energy behind. In many cases only about the highest third of the spectrum is undistorted. A method designed to circumvent the scattering effect (129) uses two separate phosphors placed very close together with the source in between. The light is taken out of the phosphor edges and is viewed by one photomultiplier. This gives in effect 4- geometry. The method is quite good for pure A emitters with none of the 3 transitions in coincidence, since electrons scattered out of one phosphor will be stopped by the other phosphor. The method is restricted since the 7 background cannot be subtracted, and also since coincident radiations will tend to add, producing a pulse representing the sum. In 1950, P. R. Bell suggested an arrangement known as the Hollow Crystal Spectrometer (69, 108, 109, 129). In this method, the 3 rays are collimated into the bottom of a conical hole in the phosphor. The scattered P rays will then have small probability of escaping by way of the entrance hole. Bell found that considerable scattering occurred in the air inside the crystal. This was eliminated by either evacuating the equipment or else using atH2 atmosphere. Hydrogen was more satisfactory since there was little difference between the pulse

height distribution in H2 and in vacuum, and the organic crystals such as anthracene and stilbene volatize slowly in a vacuum. While the method seems excellent, only one article was found in the literature up to April, 1956 that used the hollow crystal in a spectrometer (111) other than review articles covering the initial work done by Bell and co45 workers on the spectrum of Ca4 Since the 180~ magnetic spectrometer (51) in this laboratory was not designed for accurate spectral measurements, a new spectrometer was required. It was deemed necessary that the spectrometer have the following features: 1. Since the isotopes to be investigated will be cyclotron produced, the instrument must be able to work with small amounts of activity so as not to require excessively long bombardment times. 2. The instrument should be simple to construct and relatively inexpensive. 3. It must be possible to run a spectrum immediately after the source is prepared, as some of the isotopes of interest will have half-lives of the order of minutes. 4. The instrument should not require special source preparation such as vacuum plating, which would be time-consuming. 5. The restriction to purte emitters must not be required. 6. The instrument should be easily adaptable to the measurements of d spectra in coincidence with r rays. 7. The resolving time of the instrument must be less than 107 sec.

-538. It must be possible to change samples rapidly, so that calibration of the instrument may be dmne before and after the spectrum is run, if necessary. 9. It should be possible to examine conversion electron peaks with some degree of accuracy. 10. It must be possible to differentiate between allowed and forbidden spectral shapes. 11. Accurate end point values, and relative intensities of different B groups in a complex spectrum must be obtainable. While scintillation spectrometers have relatively poor resolution, this may be corrected for, as is shown in the later section. The 3 distribution is a smoothly varying function, and only the measurement of conversion electron peaks will suffer greatly due to poor resolution. Since a good two-channel 7-ray spectrometer was available, it appeared that a P-scintillation spectrometer would be best since only the phosphor would be needed. And of the scintillation techniques available, it seemed that a reinvestigation and modification of Bell's hollow crystal spectrometer would be the best approach. It was decided to have the spectrometer operate in an air atmosphere rather than using hydrogen gas or a vacuum because of the time requirements mentioned above. A few experiments were run using a helium atmosphere, and the results showed little improvement over air. The only serious drawback in this method is that positron spectra cannot be accurately analyzed due to the annihilation radiation unless a special equipment set-up is made.

A. Description of Hollow Scintillator P-Ray Spectrometer 1. Phosphors. The first problem is the choice of the phosphor. It is desired to have as large a pulse height per Mev lost in the phosphor as possible, since the resolution is approximately inversely proportional to the square root of the number of photoelectrons in the pulse. Resolution, as used here, is defined as the full width of the pulse distribution at half-maximum counting rate for monoenergetic electrons. Table VII lists various features of a few scintillators. TABLE VII. CHARACTERISTICS OF SOME SCINTILLATORS Material Rel. Light Decay Melting Yield Const.(sec) Point ( Sodium Iodide 2.0 25.0 x 108 3.67 651 -8 Anthracene 1.0 3.0 x 10 1.25 217 Stilbene 0.65 0.8 x 0o 1.16 126 Sintilon Plastic 0.6 0.8 x 10-8 1.05 10 While NaI crystals produce the largest pulses, scattering of electrons increases rapidly with atomic number and only about 10-20% of the electrons impinging on a flat NaI crystal will lose their entire energy in the crystal. This may be compared to 90% for anthracene (129). However, organic crystals are expensive, difficult to machine, and subject to shock. Therefore, it was decided -54

-55to try a plastic phosphor. The plastic phosphor chosen was Sintilon, produced by National Radiac Company, since it appeared to have properties at least as good as other plastics available. A detailed comparison of various scintillators is given by Dhar (36) and others (11, 32). A discussion of the fabrication and assembly of the hollow scintillators used here is given in Appendix A. The result of this work was the construction of three detectors to be used with n-ray energies up to maximums of 1.5, 2.2, and 3.6 Mev, respectively. Each phosphor was constructed in two parts which consisted of a flat, solid, right circular cylinder used as the base, and a top piece consisting of a hollowed-out cylinder or truncated cone. Figure 13 shows the drawings for the three phosphors which appear in Figure 14. A phototube holder and lead shield constructed to hold the detectors is shown in Figure 15. The radioactive samples whose spectra are to be examined should be carrier free, or else have as high a specific activity as possible. The samples are prepared as follows: 1. An aluminum disc, 1" in diameter with a 7/16" hole punched in the center, is used. The disc is usually v 10-20 mils thick. A piece of l/4-mil Mylar coated with a thin aluminum layer is used as a backing film. This is glued so that the organic film is next to the disc. 2. A solution of the radioactive substance is evaporated onto the organic side of the film in the center of the disc. 3. The sample is covered with a layer of 1/4-mil Mylar or Teflon film. Figure 16 shows a cross-section view of a mounted B detector. The purpose of the aluminum film evaporated onto the Mylar film is to prevent

1.5 MEV 2,2 MEV 3.6 MEV DETECTOR DETECTOR DETECTOR 76" 9i" I 8 ^ ^8" _ 1"L I 9 1 3nl1 9 1 l fix///////i//// 7 Hx 1 Ell 5,, i4 i6 "Q 16 8 __ Figue1.Cosscin ftetrehlo-lsi tnit_ 41^ Irr I I rr t- l321 -i ^ f' I, I 4' 5-2 16'5-2~ 8 __ I I ~ K-/ ----- I I Fr 1 Cr s- o th 2ti 2.2~^ 2 Figure 13. Cross-sections of the three hllow plastic s-intillators

-57~~~~~~~~r CI: ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ ~~~~~~~~~~~~~~ "~~~~~~~Msispi Figure 14. The Three Hollow Plastic Detectors | _ E | | | |. S t~~~~~~~~~~~~~~~~~m sS>.,s it-Ss' sS.9. -;-ig~~~~~~~~~~~~~~~~~~~~~~~it-:igg i:: - E g i E -'iiii':i~~a~O';i:: 01 03Wat 1 | li~~~~iiii~lX~ ss-i&>Sz Se Ss zsssv5s..'s. i's.'i~gE'ssE':S'SSEsstssS''s'ltiusss''older:"nd s~ead za.SahtizaSS;, -a.. B i~z ss si! i igi s'l sE s..8z sasz~sBea.>sot8 t Sg.8ztzlaBaaae $Z$88Z$......... XB>*B B~~~~~gs'? s i Es sssssa~~~~~~~~~s? s. i z'. z B.az~~~~~~~~~z~~ss. aS~~~~c? aY. YZ............. B....s. ss... ss~~~~~~~~~~~~~~~~~s. g 3 a aa eS e Bs a B za*~~~~~~~~~~~~........ S iz?.^;? SB~~~~~~~~~~g. Ss SSS.-. SZ S s- 5 s- S - s s - s? B Z>,,, aa?,>,aesi s~~~~~~~~~~~~-...I...-... —........ ess.>.Z;><sissii~~~~~s eB~~ss~~ss~s?*s9 i sss~~~~s ssi~~z~~ss?~z~~s e ZZ a; a.Ss;a s.?.;.S..se is.,B~~~~~~~~~sB..........s..................B 3 -ss~~~~~ss~~~s~~i s. ss. sz~~~~~~~~~s. ss i..ss.?.ss gs. sa ZY s a zss...YizB >. s~~~~~~~.............................. Zz~~~e~~sbs;?s.'skaes 32'-g —'S' —-,l~~~~~~~~~~~~~~~~~a>.S.."2;S s'ES SYM-................... -. }.;.;..;; iii..... Z..? >SB.Ss r..B..;.S;EBBN..aB&................................... Figare 1>. Phototube Holder and Lead Shield for~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~....... Hollow Scintillator Beta-ray Spectrometer~~~~~~~~~~~~~~~~~~~~~~~~.......

-58light from entering or leaving the detector through the hole in the top. Hollow Plastic 0.0015"AI Foil Phosphor Reflector AI Collar Steel Phototube Housing Black -Phototu be Felt Figure 16. Cross-section of beta-ray detector mounting 2. Operation If coincidence spectra are to be run, the { detector is simply mounted on a phototube in the housing described in the section on the 7-ray spectrometer. Otherwise, the detector is mounted as shown in the previous section. The lucite sample holder is then centered over

-59the phosphor and placed in contact with it, with the locking nuts holding it tightly in position. A sample is put into position with the aluminum coated film away from the hole over which is placed the lucite retaining ring. A light-tight cap, made of rubber sheet 1/4" thick which also serves to shield the phosphor from electrons knocked out of the lead shield, is placed over the top of the phototube housing. The top and bottom doors of the shield are then closed. The spectrum may now'be taken. It has been found that increasing the high voltage from 650 to 1300 volts does not affect the resolution of the spectrometer appreciably. Therefore, in order to reduce the noise in the phototube, a voltage of 830 volts is commonly used. The pulse-height analyzer channel width used depends upon the counting rate of the sample. The "window" may be opened to - 1.5 volts before the resolution, position of peak, or peak to valley ratio for the conversion peak of Cs137 are affected. With a wider window the resolution and peak to valley ratio become poorer. This is illustrated in the set of curves shown in Figure 17. Table VIII shows the peak to valley ratio and apparent resolution at the various window settings. TABLE VIII. RESPONSE OF C137 CONVERSION PEAK TO PULSE HEIGHT ANALYZER WINDOW SETTING,. i.,,..,., j,...: _, Window (volts) Peak/Valley Resolution % 0.10 5/1 14.0 0.20 6/1 14.0 O. 40 5/1 13.5 1.00 5.5/1 14.0 1.75 5/1 15.0

-603000I i RESOLUTION AS FUNCTION OF WINDOW WIDTH Cs'37 2500 2500 WINDOW WIDTH PULSE HEIGHT UNITS A 175 0 100 9nnn~2000 40 8 2000 20 o I10 1500 _0 o 1000 500 600 700 800 900 1000 PULSE HEIGHT Figure 17. Effect of window width on resolution of Cs137 conversion peak. Ordinate scale for top two curves is 2.5 times indicated scale.

-613. Calibration. Since the linear amplifiers are dependably linear only from 10 to X 80 volts, the gain must be adjusted so that the spectrum to be analyzed, along with the calibration points of In, Sn, and Cs137 falls in this range. The spectrometer is first calibrated by means of the conversion electron peaks of known energy shown in Table IX. TABLE IX. CONVERSION ELECTRON PEAKS USED IN CALIBRATION Isotope Energy (Mev) n14 0.162 Sni13 0.364 Cs137 0.624 The stated energies are not, however, the energies of the conversion electrons striking the phosphor. Corrections must be made for the amount of energy lost in the 1/4-mil Mylar film covering the sample, and also in the air in between the sample and the phosphor. This can be done if the rate at which energy is lost, as a function of energy, as the D particle passes through an absorber is known. Using Equation (12), dE 2icNe Z log dx 2 2' \ (l + ul + log,[u(u2 + 2u)] - fl -( ][ log2+ 1 } (12)

-62and information given in Reference 5, the curves in Figure 18 were calculated. -dE/dx is in units of Mev/mg cm-2, and u = E/mc2. The rest of the symbols have their usual meaning. While the curves in Figure 18 are only shown to 1.0 Mev, they may be extrapolated out to X 10 Mev with little error. No attempt was made to calculate a mean path length in the air within each phosphor. It was simply assumed that a good average distance was X 1.6 cm. The value of 1.18 mg/cm3 was used as the density of air. For the Mylar film a surface density of 0.635 mg/cm2 (of aluminum) was used. The pulse height analyzer was calibrated by assuming a linear relation between the dial reading and the observed energy of the conversion peaks in the range from 0.16 to 0.62 Mev. It has been found that a linear extrapolation to higher energies is valid at least to 3.6 Mev, the endpoint of the most energetic spectrum measured in this work. For lower energies a linear extrapolation is made to 0.1 Mev, and then a smooth curve is drawn from the end of the line to the origin. It must be remembered that this is at best a guess, and this portion of the spectrum cannot be relied upon. However, in most cases this portion of the curve falls in the non-linear lower range of the amplifier and should be discarded anyway. As a further check on the calibration, spectra with known maximum energies and shapes were obtained either before or after an unknown is run. After each run the energy assigned to each pulse height amplifier reading the curves in Figure 18 to obtain "true" energies.

-63RATE OF ENERGY LOSS FOR ELECTRONS IN 7 AIR AND ALUMINUM 6^~ ^ AIR o ALUMINUM ro 5 g: LI 0 0.1 0,2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 MEV Figure 18. Rate of energy loss as a function of energy for electrons in aluminum and air

A word should be said about drift in the equipment. As long as the ambient temperature remains reasonably constant, say within 5~ or so, no drift of over' 1.5% was found during three continuous weeks of operation. B. Performance of Hollow Scintillator P-Ray Spectrometer An investigation was made to determine the capabilities and limitations of the hollow scintillator spectrometer. Near the end of this section a comparison is made with a flat detector constructed from the same scintillating plastic. The last subsection deals with the problem of instrument resolution. None of the spectra presented up to this last section have been corrected for instrument resolution except near their endpoints. 1. Effect of Detector Size As mentioned earlier, the three detectors available are for maximum 3 energies of 1.5, 2.2, and 3.6 Mev. Any y radiation accompanying the 3 particles must be subtracted from the gross curve to yield the 3 component. Since error will invariably be introduced by this subtraction procedure, it is preferable to have the 7 component as small as possible. Figure 19 shows the 7-ray background in each detector with respect to a normalized Cs137 P-ray spectrum. While in general quite small, the 7-ray background is comparable to the P-ray spectrum near Ema. It is to be understood that the "7 background" refers to the gross spectrum minus the 3 spectrum. Hence, the "y background" includes also the cosmic and background radiation, along with the noise in the amplifier and phototube. No distinction will be made between the true

10,000 GAMMA BACKGROUND COMPARISON * BETA-RAY SPECTRUM 8 000~ Cs37 GAMMA-RAY BACKGROUNDS Ui^~~~~~'~o -1.5 MEV DETECTOR H 0-2.2 MEV DETECTOR Z A-3.6 MEV DETECTOR 0 6 00 LU~ H 4 000 z:D 0 2q00 0.1.2.3.4.5.6.7 ENERGY (MEV) Figure 19. Gamma-ray backgrounds in the three plastic detectors

-667 background and the term as used above. It is interesting to note the detection efficiency of the phosphor for ( rays and y rays. Assuming an equal number of 3 rays and 7 rays in Cs137, the ratio of the y to 3 efficiencies in the usual experimental set-up is roughly 0.23 for the 3.6-Mev phosphor, and 0.15 and 0.11 for the 2.2 and 1.5-Mev phosphors, respectively. Figure 19 suggests that for maximum 3-ray energies less than 1.5 Mev the smallest detector should be used. There is, however, another point to keep in mind, i.e. the phosphor shape. Two important effects are: (1) the thickness of the phosphor, and (2) the amount of optical coupling between the top and bottom pieces of the phosphor. As the size of the phosphor increases the resolving power and the light output decreases. The change in resolving power is only in the order of 1 or 2% for the phosphors used here, but the decrease in light output may be as much as 5%. As the volume of the phosphor decreases, for a given base diameter the walls of the hollow top of the phosphor become thinner. This means that the optical interface becomes smaller. This interface should be as large as possible to reduce the scattering in the phosphor. As it is, the interface is located over the periphery of the photo-cathode and reducing the wall thickness will just make things worse. Since these two effects are in opposite directions, the 2.2-Mev phosphor is used for energies up to X 2 Mev. The 1.5-Mev phosphor is used only when the small difference in 7 backgrounds becomes important. In Figure 20 the Kurie plots for the three Cs137 spectra connected with Figure 19 are shown. All three spectra have essentially

-67COMPARISON OF - o SCINTILLATORS o Cs137 0'A N0 0o 0o A-.5MEV ~, B-3.6 MEV C-2.2 MEV 0. 0. 2 0.3 0. 4 0,5 ENERGY (MEV) Figure 20. Kurie plots of Cs137 from each detector

the same endpoint value of 0.51 Mev as compared with the literature value of 0.517 Mev. As expected, the 2.2-Mev detector yields the best Kurie plot, which is linear down to about 0.18 Mev. Each curve was corrected by the unique first forbidden correction factor a1(W), where Wo =2.01. To find out how the linearity behaves as a function of Ems, a similar set of curves were run on Pm147 and p32. Pm147 yields a single allowed P-ray spectrum with Ea = 0.224 Mev. Figure 21 shows the resulting Kurie plots. There seems to be little difference between the two detectors. p32 also has a single allowed 8 transition with Emax = 1.7 Mev. No y rays are present. Figure 22 shows the Kurie plots. The source contained P33 which caused part of the upward curvature at lower energies. As a final means of comparison, the resolution of the three phosphors for the Cs137 conversion electron peak are listed in Table X. TABLE X. RESOLUTION OF THE THREE PHOSPHORS FOR THE Cs137 CONVERSION ELECTRON PEAK Phosphor % Resolution Uncorr. 1.5 Mev 14.7 2.2 Mev 14.0 3.6 Mev 15.8 These resolutions have not been corrected for the higher energy 3 ray in Cs137~ nor for the smearing of the low energy edge of the peak. The true % resolution would be somewhat smaller in each case.

-69I I I I \COMPARISON OF - o\ SCINTILLATORS - 0\ Pm147 NF 0-1.5 MEV -2.2 MEV 0.02 0,06 0.10 0.14 0.18 0.22 ENERGY (MEV) Figure 21. Kurie plots of Pm147 from the 1.5-Mev and the 2.2-Mev detectors

-70o I I 1 1 COMPARISON OF SCINTILLATORS o p32 0 0 0 0 0r~ * -2.2MEV ~ \o 0 -3.6 MEV F* I 0.5 1.0 1.5 2.0 ENERGY (MEV) Figure 22. Kurie plots of P32 from the 2.2-Mev and the 3.6-Mev dete ctors

-71The evidence presented above indicates that the 2.2-Mev phosphor is the best for general use for B rays up to 2.2 Mev. The advantage is only slight, however, and good results may be obtained with both of the detectors. 2. Effect of Sample Backing Essentially any amount of backing upon which the sample rests will cause some distortion of the a spectrum due to backscattering. It is also of interest to compare light reflecting backings with transparent ones. It would naturally be supposed that a reflecting backing would be better since any light loss would distort the spectrum and decrease the resolving power. But the magnitude of this effect was not known. In Figure 23 the Kurie plots of Cs137 are shown for the following backing materials: 1. 1/4-mil Mylar (0.64 mg/cm2). 2. 1/4-mil Mylar coated with aluminum (0.7 mg/cm2). 3. One 6-mil aluminum disc behind l/4-mil Mylar (38.1 mg/cm2). 4. Nine 6-mil aluminum discs behind 1/4-mil Mylar ("infinite" backing). Comparing the transparent Mylar backing aloane with the same backing covered with one aluminum disc, it would seem that the distortion caused by light loss in the first case is about equal to the backscattering distortion in the second case. The aluminum-coated Myla is by far the best backing.

-72* o BACK-SCATTERING * ~ \o EFFECT - 0 Cs137 o To 0 Elcr. D \I aF O-MYLAR+AI FILM *-9 Al DISCS a-1 AI DISC -MYLAR (TRANSPARENT) | 1.0 1.2 1.4 1.6 1.8 2.0 ENERGY (RELATIVISTIC UNITS) Figure 23. Effect of sample backing on shape of the Kurie plot of cS8137. For energy in Mev see Appendix B.

-733. Backscattering From Sample Holder. When the phototube assembly was first constructed it was not known if backscattering from the sample holder would be a problem. Two sample holders were built, one of aluminum and the other of lucite, both of the same design. The design was such as to minimize backscattering as much as possible. Figure 24 shows a P32 spectrum obtained with each sample holder. As may be seen, the aluminum sample holder caused enough backscattering to make the Kurie plot linear from v 0.65 Mev to Emax, whereas the other Kurie plot is linear down to, 0.45 Mev. 4. Backscattering Shields and Collimators. A disc of 18-mil copper with a 7/16" hole in its center was placed over the phosphor to prevent electrons from striking the phosphor except through the entrance hole. The resulting Kurie plots showed very little if any improvement. Hence, this phase of backscattering appears to be a negligible effect in the geometry used. Another possible detrimental effect is that the optical interface between the top and the bottom of the crystal might cause some scattering and loss of light. This would tend to cause the Kurie plot to bend up at low energies. To investigate this effect a collimator was made from a disc of 18-mil copper with a 3/16" hole in its center. This was placed between the source and the detector so that more of the electrons would strike the base of the detector. Again, the Kurie plots showed little if any improvement. This would indicate that the optical interface was good, and that light loss in the phosphors was not important.

-74p32 0 SAMPLE HOLDER O O Al 0 * LUCITE \ 0 O I 0 * * \ N 0 0 0.5 1.0 1.5 2.0 0 0.5 1.0 1.5 2.0 ENERGY (MEV) Figure 24. Effect of backscattering from sample holder on the Kurie plot of P32

-755. Thick Samples. Occasionally carrier-free 3 sources are not available, due perhaps to difficulties in the chemical separation or more commonly when the radioactive isotopes produced in a nuclear reaction have the same atomic number as the target material. The best one can hope for then is to have a high specific activity. Since any weighable sample will distort the 3 spectrum, it is of interest to determine the magnitude of this effect. Finely powdered BaC03 was mixed with a carrier-free solution of Cs7. A more or less uniform sample. 3.6 mg/cm in surface density was then made by evaporating a portion of the BaCO3 suspension onto the usual aluminum coated Mylar film. The sample was covered with a layer of 1/4-mil Mylar. Figure 25 compares the Kurie plot for the thick sample with a plot obtained for a carrier-free sample mounted in the same way. It can be seen that while the thick sample shows more distortion than the other, it still yields a good curve down to about 1.6 relativistic units (W = E + 1) or X 0.3 Mev. 6. Gamma Attenuation in Absorber. To determine the y background an absorber is placed between the source and the detector. Since the absorber normally used is 4.53 x 103 mg/cm2 thick, it was thought that the attenuation of lowenergy y rays, say up to 0.1 Mev, might become important. A comparison between the normal absorber and an18-mil copper disc showed no detectable difference in the Kurie plots for Cs37. The normal absorber is therefore the one that is routinely used.

o EFFECT OF o SAMPLE THICKNESS o Cs'37 0 0o 0 3.6 mg/cm2 Lo *I CARRIER FREE 0 0 VaF.... 0 1.0 1.2 1.4 1.6 1.8 2.0 2.2 ENERGY (RELATIVISTIC UNITS) Figure 25. Effect of sample thickness on the Kurie plot of Cs137 For energy in Mev see Appendix B, Figure 67.

-777. Amplifier Overload and Drift. Occasionally it is of interest to examine only the lower portion of a spectrum, such as in the case of a complex spectrum. For the sake of accuracy, as much of the amplifier range should be used as possible. For example, if a spectrum contains a 2.0 Mev and a 0.5-Mev P ray, and the amplifier is adjusted so that the entire spectrum is viewed, then only the lower quarter of the range contains the weaker 3 ray. Due to the non-linearity of the amplifier, about 1/3 of the lower quarter will be distorted. If the amplification were changed so that the 0.5-Mev portion occupied one-half the available range, then only about 1/5 of the lower energy spectrum would be distorted. In the latter case, the amplifier would overload at the higher end. It is of interest to know how much distortion this overload would cause, and if it is still possible to obtain useable Kurie plots. Normally only about 80%0 of the amplifier range is used since the upper 20% is somewhat non-linear. A P spectrum was run where the amplification was such that the endpoint of the spectrum fell outside the full range of the amplifier by X 15%. The Kurie plot is shown in Figure 26. It can be seen that while the upper one-third of the spectrum is distorted a straight Kurie plot is still obtainable with an endpoint within 2% of the correct value. Amplifier overload is a relatively easy type of distortion to recognize. It is not readily confused with a forbidden spectral shape. A more subtle type of distortion occurs when drift occurs in the system. Sometimes the drift is obvious as shown in the second Kurie plot in Figure 26. Here a power failure occurred for approximately

-78five minutes. Although the equipment was allowed to stabilize for over an hour, the gain of the system had definitely changed. By far the most common effect of drift is to yield a Kurie plot that is almost but not quite straight, with the wrong endpoint value. An example is the top plot in Figure 26. The solid line shown is only one of many straight lines, each of which would fit the data about equally as well. Corresponding endpoint values would vary by 7% or more. Probably the worst danger in the analysis of a driftdistorted spectrum is the temptation to "resolve" the curve into several fictitious components. For example, the top plot in Figure 26 might be "resolved" to give a weaker component with a forbidden shape and an endpoint of ~ 0.9 Mev. Checking the energy calibration before and after a spectrum is taken seems to be the only answer. 8. Range of Spectrometer The highest energy n-ray spectrum examined in this work belonged to Y92 Ema = 3.60 Mev, and will be discussed in a later section. The weakest 3 transition analyzed was that of S35 where Emax = 0.17 Mev. The resulting Kurie plot using the 2.2-Mev phosphor is shown in Figure 27. An intermediate range transition of use in calibration is that of Inll where Eax = 1.98 Mev. This allowed spectrum and Kurie plot is shown in Figure 28. For most purposes the size of the phosphor is the only limitation on the maximum D energy that may be examined (below X 10 Mev) with this spectrometer. This is not true for the lower energies. All phosphors including NaI tend to be non-linear at low energies. Table XI lists a few phosphors and the energy below which their response is nonlinear (11, 32, 36).

-79o AMPLIFIER OVERLOAD AND DRIFT o p32 0 00 0 A-Drift F F o ~ 0 0 B- Power failure C-Overload effect 0.5 1.0 1.5 2.0 ENERGY (MEV) Figure 26. Effect of amplifier overload and drift on the Kurie plot of P32.

-80KURIE PLOT OF S35 0 1/o 0 VF 0,05.1.15.20 ENERGY (MEV) Figure 27. Kurie plot of S35

200O~ I n"114 ~KURIE PLOT OF BETA SPECTRUM I I114 LJU D 15,00 Z F 10000 ZD 4,.08 1.2 16 2.0 o 0.4 MEV ENERGY (MEV) 51000 CURVE x 4 0 0 10 20 30 40 50 60 70 80 VOLTS Figure 28. Beta spectrum and Kurie plot of In114

-82TABLE XI. LOW ENERGY RESPONSE OF SOME PHOSPHORS Phosphor Lowest Energy in Linear Portion of Response Curve (kev) NaI 1 Anthracene 1.00 x 102 Stilbene 1.25 x 102 Terphenyl 2.40 x 102 Sintilon 1 x 102 9. Forbidden Spectra. It is imperative that a spectrometer be accurate enough to distinguish between the shapes of allowed and forbidden spectra. With this information a great deal more may be said about the nuclear transition generating the spectrum. Without this information the resolution of complex spectra into components may be difficult if not impossible. To evaluate the spectrometer used here a number of "unique" forbidden spectra have been examined. In the cases where the spectra are complex, only the highest energy component will be shown. The resolution of complex spectra will be discussed in the next section. 137 In Figures 29, 30, 31, and 32 the Kurie plots for Cs Rb, Y9, and Y1 are presented. Both the plot which is produced when the spectrum is assumed to be allowed, and also the plot corrected by the unique first forbidden correction term are shown. The pronounced difference in shape, i.e. the skewing of the spectrum toward higher

-83EFFECT OF FORBIDDEN CORRECTION FACTOR Cs137 0 0 0- Uncorrected 0 *-Corrected 0 0 0 ~1/F 0 00 *xen 00 0.0 1..... 1.0 1.2 1.4 1.6 1.8 2.0 ENERGY (RELATIVISTIC UNITS) Figure 29. Kurie plot Cs137 showing effect of forbidden correction term. For energy in Mev see Appendix B, Figure 67

-84EFFECT OF FORBIDDEN CORRECTION FACTOR Rb86 - U ncorrected - Corrected 0 0 0 0 NR o0o VF 0 0o I X 2.5 3.0 3.5 4.0 4.5 ENERGY (RELATIVISTIC UNITS) Figure 30. Kurie plot of Rb8 showing effect of forbidden correction term. For energy in Mev see Appendix B, Figure 67.

I I! EFFECT OF FORBIDDEN CORRECTION FACTOR y 90 0 -Uncorrected 00 * -Corrected 0 _0 r0 N K 0 aF 0 _ 0 0 c 00000 2.20 0 2.2 3.0 4.0 5.0 ENERGY (RELATIVISTIC UNITS) Figure 31. Kurie Plot of Y90 showing effect of forbidden correction term. For energy in Mev see Appendix B, Figure 67.

EFFECT OF FORBIDDEN o CORRECTION FACTOR o y91 0 0-Uncorrected 0 0 - Corrected 0 IN Oo0 F 0 1.5 2.0 3.0 4.0 ENERGY (RELATIVISTIC UNITS) Figure 32. Kurie plot of y91 showing effect of forbidden correction term. For energy in Mev see Appendix B, Figure 67.

-87energies, is typical of forbidden spectra. In the case of a complex spectrum where one or more of the 3 groups is forbidden, it is obvious that the spectrum could not be resolved with any accuracy at all unless the forbiddenness is taken into consideration. 10. Complex Spectra. Besides knowing the Eax and the forbiddenness of a transition, it is also important to know the branching ratios in a complex spectrum. These ratios must be known, for example, if the log (ft) values are to be calculated. One means of measuring these ratios is to compare the areas under the respective Kurie plots. Naturally, these areas must be compared on an equal basis. If two components are present, say, and one is forbidden, then the area under the plot corrected far forbiddenness cannot be directly compared with the area under the allowed plot since the forbidden spectrum has been multiplied by a correction factor. In this case, the areas are compared when each spectrum is assumed to be allowed. One might alternately calculate back from each Kurie plot to find the true spectrum and then compare those areas. To simulate a complex spectrum, a pml47 and a Cs137 sample were examined simultaneously. The samples were on separate plates placed such that the Pm147 sample was closest to the phosphor. Here the 2.2-Mev detector was used. The spectrum was analyzed as if it were 137 due entirely to Cs. After correcting for forbiddenness the Kurie plot shown in Figure 33 was made. The Cs137 component was subtracted out, and the remaining plot corrected back to give the Pm147 spectrum. The resulting Pm147 Kurie plot appears in Figure 34. The areas under

-88KURIE PLOT OF Cs137+ Pm147 MIXTURE 0 0 0 0 0 0aF 0 0 0 0 137 Cs 1.0 1.5 2.0 ENERGY (RELATIVISTIC) Figure 33. Kurie Plot of Cs137 + Pm147 mixture

-89KURIE PLOT OF oo Pm147 FROM Cs137+ Pm147 MIXTURE II I X....05 0.10 0.15 0,20 ENERGY (MEV) Figure 34. Kurie plot of Pm147 resolved from Cs137 + Pm147 mixture

-90both Kurie plots were compared with the Kurie plots taken individually on each sample. The ratio of the areas in the mixture was within 4% of the ratio from the individual samples. This is considered to be within the experimental error since the phosphor had been remounted before the mixture was analyzed. Consider next the complex spectrum of Rb8. Of the e rays present, the first has a reported energy of 0.68 Mev and a relative abundance of \ 10%. The second E ray is ^ 90% abundant with Emax = 1.77 Mev. The branching ratio is not accurately known; some relative values for the higher energy component go as low as 80% and as high as 92%. While the higher energy component is forbidden, several investigators have found an allowed shape for the lower transition. Figure 35 shows the Kurie plots obtained with the 3.6-Mev phosphor. The branching of the lower component as determined here is u 12%. The endpoint of 0.60 Mev is considerably lower than is indicated in the literature. Since the higher energy component's endpoint is within 1% of the literature value of 1.77 Mev, it is assumed that any errors must lie in the graphical resolution. In a spectrum like this, coincidence techniques are of great value. 11. Positron Emitters. The P spectrum of a positron emitter may also be obtained using the scintillation detector if a three-channel coincidence spectrometer is available. With positron emitters the difficulty lies in the presence of annihilation radiation. When the positron reaches the end of its range and is annihilated, the resulting energy is given off primarily in the form of two 7 rays, each 0.511 Mev in energy, and emerging

-91KURIE PLOT OF 0 LOWER ENERGY COMPONENT o0 0 O 0 1.5 2.0 ~ \o (RELATIVISTIC UNITS) 0O aF Rb86 2 3 4 5 ENERGY (RELATIVISTIC UNITS) Figure 35. Kurie plot of Rb86 For energy in Mev see Appendix B, Figure 67.

-92at 180~ to each other. Since the y rays will originate in the detector itself, there is no simple way of determining the y background. Furthermore, since the y radiation is in coincidence with the 3 pulse, an unknown amount of energy, from zero up to the Compton Edge (( 0.34 Mev) of the 0.511-Mev y pulse distribution, will be added to the energy of each 3 pulse. Actually no positron spectra were measured in this work, but a method for obtaining these spectra is given below. Using the three-channel coincidence arrangement shown in Figure 36, the distortion will be for the most part eliminated. The detectors marked 1 and 2 are NaI(T1) crystals, and their associated pulse-height analyzers are set on the 0.511-Mev y energy. When the positron is annihilated in the P detector (say at the place marked with an x), the two 7 rays will have a certain probability, dependent upon the geometry, of one entering each NaI(Tl) detector as shown. Thus, disregarding chance coincidences, if both detector 1 and 2 simultaneously "see" a 0.511-Mev y ray, this means the 3 pulse is undistorted since neither y ray lost any energy going through the 3 detector. By setting all three channels in coincidence only those positrons whose energy pulses are undistorted will be detected, and the a spectrum may be scanned in the usual way. In this arrangement positron spectra may be determined in the presence of accompanying negatron radiation. As in all coincidence work relatively intense sources or long counting times are necessary in order to achieve good statistical results. Particularly so in this case, since in order to register a coincidence count both 7 rays must pass through the B detector without

-93losing energy, both must strike a y detector, and both must produce a pulse in the photopeak region. Using 1" x 1-1/2" NaI(Tl) crystals and the 2.2-Mev B detector previously described in a reasonable geometry, one might expect to find a coincidence counting rate down by a factor of - 103 from the true 3 counting rate. Na I(TI) Na I(TI) POSITRON BETA SOURCE DETECTOR Figure 36. Three-channel coincidence arrangement for positron emitters

-9412. Conversion Coefficients. An excited nucleus, which would normally make a transition to a lower state by emitting a y ray, may lose energy by another process. The alternative process, called internal conversion, involves the ejection of one of the orbital electrons. Here the electron appears with the energy of the y ray minus the binding energy of the shell from which it was ejected. A third process wherein an electron-positron pair is emitted occurs relatively seldom compared to internal conversion. It is primarily found at high energies and/or low atomic numbers. The transition energy has to be greater than 1.02 Mev. For 0 - 0 transitions no conversion or 7-ray emission is possible and pair emission becomes important. We shall define the conversion coefficient a as follows: a = aNe = K + aL + m +.. (13) N7 where Ne is the number of conversion electrons emitted, and N7 is the number of 7 rays. As defined, a is the total conversion coefficient, and is the sum of the conversions coefficients for each electron shell as indicated by the subscripts. It is understood that aL, for example, is the sum of the subshell coefficients: aL = +L + LII +, (14) Conversion coefficients are of interest because they are quite sensitive to the following parameters: 1. The atomic number of the emitter, 2. The transition energy. 3. The particular shell or subshell.

-954. The multipolarity L of the competing y radiation. 5. The character of the nuclear transition, i.e. electric or magnetic. Gamna rays may be classified by multipole orders L, according to the angular momentum carried off by each quantum. For each order there are two subclasses electric 2L pole (EL) and magnetic 2L pole (ML). These differ with respect to parity as is shown in the table below. TABLE XII. POSSIBLE MULTIPOLARITIES FOR GAMMA RADIATIONS Parity Angular Momentum Change Change 0 or 1 2 3 4 5 No M1(E2) E2(M3) M3(E4) E4(M5) M5(E6) Yes El(M2) M2(E3) E3(M4) M4(E5) E5(M6) The transition in parenthesis is usually not significant. Once L is fixed, the character of the nuclear transition uniquely fixes the parity change. Furthermore, the transition probability is a function of L and the character of the transition. As is expected, the larger the number of units of angular momentum the quantum must carry off, the smaller is the transition probability and the longer is the half-life of the state. In Cs137, where the transition is the type M4, a delayed state exists with a half-life of about 2.6 minutes. Conversion coefficients are measured in a variety of ways depending upon the particular isotope involved. If no $ rays occur the

-96number of conversion electrons are compared with the number of x-rays or y rays emitted. In other cases, the number of P rays emitted may be used to calculate the 7-ray intensity. After the coefficient is measured, it is compared with tables such as those by M. E. Rose, et al. (119) in order to deduce the type of transition involved. In most cases the instrument described here, due to its relatively poor resolution, will not be able to resolve K electrons from L or M electrons. This is unfortunate since CK, oLI, and LIII ratios aL aLII CLI are generally more sensitive in energy dependence than the absolute conversion coefficients. Nevertheless, in many cases the coefficients differ sufficiently so that only the total conversion coefficient need be used. Using Figure 40 and assuming a symmetrical conversion peak, the conversion coefficient for Cs137 was measured by integrating under the peak and the 3 spectrum. Assuming an mK ratio of 5 (26), the aK for the corrected curve was found to be 0.095. For the uncorrected curve (K was 0.082. These may be compared to 0.094 as given by Rose's tables for M4 radiation, and 0.097 as determined experimentally by Waggoner (26). The coefficients from both the corrected and uncorrected curves are good enough to determine the type of transition if the transition is assumed to be "pure." If the possibility of mixed transitions is considered, then the uncorrected value would lead to erroneous results.

-9713. Comparison with Flat Detectors. As another means of evaluating the hollow-type detector, a comparison was made with the results obtained with a flat detector. The flat detector was made by using the base of the 2.2-Mev hollow scintillator covered with a 1.5-mil sheet of aluminum foil. The detector was coupled as usual to the phototube, while the source was placed 2.2 cm from the detector along its axis. Figure 37 shows the results for Cs137. The table below presents some of the information obtained from the spectra. TABLE XIII. COMPARISON OF FLAT AND HOLLOW DETECTORS Detector Resolution Conversion Peak to Valley Ratio Flat 23% 4/1 Hollow 14% 6/1 The better resolution of the hollow-type detector proves advantageous when dealing with conversion electron peaks, endpoints of spectra, complex spectra, and other cases where the slope of the spectrum changes rapidly. Also, the better resolution facilitates the recognition of forbidden shapes which usually appear smeared out with poorer resolution. It should be noticed that in Figure 37 the flat detector spectrum is badly distorted at the lower end. This is probably due to the fact that for a given angle of incidence, the lower the energy of a 3 particle the greater is its probability of scattering out of the

-98Cs37- FLAT DETECTOR 000 0 00 0 0 0 0 000 0 0 0 0 ~W ~~~0 00 - 0 0 0 0 0 0 0 0 0 wL0. 0. 0.2 0.3 0.4 0.5 0.6 CL 0 0 0000O 0O ~D ~0 0O ^0 0 0 0 0 0 1~~- ~o0 0 0 0 0 0 0 0 0 0 L i | |I I I 000? 0 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 ENERGY (MEV) Figure 37. Comparison of flat and hollow detectors

-99detector. Figure 38 compares the Kurie plots of P32 obtained on the two types of detectors. While the flat detector produces a Kurie plot that is linear over about half the range, the hollow detector's plot is linear over more than three-quarters of the range. On both curves the sharp upward curvature at the lower end is due to P33 in the sample. 14. Correction for Instrument Resolution. Because of the low resolving power of n-scintillation spectrometers, the observed spectral shape does not correspond exactly to the true shape. In fact, at both the low energy and the high energy ends the distortion usually proves to be very significant. The experimentally measured spectrum M(E) is related to the true spectrum N(E) by the equation E M(E) = max N(E)S(E)dE (15) t EN(Ei)S(Ei)AiE where S(E) is the so-called instrument profile (112). Until recently, the usual approach to this problem was to assume that S(E) could be represented by a Gaussian distribution of the form (E -E)2 S(E)= A e 2~" (16) where A is a normalizing constant, and a2 is the variance of the pulse height distribution for a monenergetic beam of electrons of energy, E. Katelle (71), and Palmer and Laslett (113) have used this approach, based on the method of Owen and Primakoff (112). In the body of the spectrum this method has uniformly failed to eliminate upturns at low energy. At

. i I I I I I. I I I' FLAT vs. HOLLOW DETECTOR p32 0 o 0 FLAT DETECTOR o * HOLLOW DETECTOR 0 O F 0 ~_ ~ 0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 ENERGY (MEV) Figure 38. Kurie plots for P32 from both flat and hollow detectors

-101the high energy end near Ea,, however, the procedure works well. A suitable correction factor is illustrated graphically in Reference 113, and has been used to correct all of the spectra shomwn in this work up to this point. The correction factor is a function of the maximum 3 energy. Therefore, in order to determine the correction, an endpoint must first be estimated. Thus, the correction serves only to indicate whether or not the observed data are donsistent with the estimated endpoint. Figure 39 shows the correction applied to the Kurie plots of Cs137 and m147. Another method for resolution correction, due to Bracewell (17), also uses a Gaussian distribution. The method has the advantage that it is simple and easy to apply. When applied to the Cs137 spectrum, the results are shown in Figure 40. Since the method is graphical, the correction is not too exact. The resolutions calculated from the peaks in Figure 40 are: Uncorrected - 14.7% Corrected - 12.7% Probably the best approach to the resolution problem was put forth by Freedman, et al., in June, 1956 (45). They made the reasonable assumption that the instrument profile was not symmetrical, and that a monoenergetic beam of electrons will produce a profile that is measurable down to zero energy. Figure 41 illustrates this point. Four conversion peaks were measured in their work, and it was found that the relative height of the tail of the curve remained

-1020 0-Uncorrected 0~ * - Corrected Cs137 RESOLUTION CORRECTION 0o VF I 0 1.2 1.4 1.6 1.8 2.0 o,ENERGY (RELATIVISTIC UNITS) 0 O- Uncorrected * - Corrected pm147 RESOLUTION CORRECTION s Jo1O.05.10.15.20.25 ENERGY (MEV) Figure 39. Effect of resolution correction near the endpoint of beta spectra. For energy in Mev see Appendix B, Figure 67.

-103II 1- 1 EFFECT OF RESOLUTION CORRECTION Cs137 6000 0 - Uncorrected * - Corrected 5000 LLU 4000 LU 1 3000 0 2000 1000 0 200 400 600 800 PULSE HEIGHT (ARBITRARY UNITS) Figure 40. Resolution correction by Bracewell

-104essentially constant, and was 7.5% of the peak intensity for each conversion peak. >H 1.0 Li/. 0.5ENERGY Figure 41. Form of experimental line shape determined by Freedmen, et al To obtain N(E) from Equation (15), a zero-order approximation to N(E) is made: M(E)/AE, where AE is the window width. The integration is performed numerically producing a new M1(E). If S(E) has been absolutely normalized, the new spectrum will have the same area as the

-105original M(E). A first approximation to the true spectrum is now obtained. N1(E)AE = M(E) (E) - M(E)] = 2M(E) - M(E) (17) Next N1(E) is put into the integral and the numerical integration repeated to obtain a second approximation. Usually two to four approximations were needed to obtain a spectrum which, when used in the integration, reproduced the experimental curve to within 2%. The above procedure has been used here to correct the Kurie plots of two spectra, and was found to be excellent. Figure 42 shows 137 the profile of the Cs conversion peak obtained by requiring coincidence with the x-ray. In this work, it was found that the height of the tail of the profile was 5.0% of the peak intensity, for the 2.2-Mev detector. Knowing that the half-width of the peak is a function of energy (113), the following equation was used to construct profiles at any given energy once one profile was measured: 1/2 Wl/2 = CE/ (18) where W1/2 is the half-width, C is a constant, and E is energy. Figures 43 and 44 show the effect of this type of resolution correction on the Kurie plots of Cs137 and 91. In the case of Cs137 the Kurie plot is straightened out down to 0.05 Mev. For 9 the Kurie plot now is straight down to 1.2 relativistic units or about 0.1 Mev. 15. Summary. A -ray spectrometer of the scintillation type has been designed, built, and evaluated. Using a hollow plastic detector the resolution has

500 I EXPERIMENTAL D I CONVERSION PEAK 400 PROFILE FOR Cs137 LU 0 300 0 100 200 300 400 500 600 70 o I200 0 ( 0 0 100 200 300 400 500 600 700 PULSE HEIGHT Figure 42. Experimental conversion peak profile for Cs137

-107RESOLUTION CORRECTION oUSING NON-GAUSSIAN LINE PROFILE-Cs137 0 0 UNCORRECTED \ 0 CORRECTED 0 aF 0 0.1 0.2 0.3 0.4 0.5 ENERGY (MEV) Figure 43. Effe7ct of Non-Gaussian Resolution correction on Cs1 Kurie plot

RESOLUTION CORRECTION USING - NON-GAUSSIAN LINE PROFILE Y9' ~0~~~~ 0 UNCORRECTED * CORRECTED __ 0 aF 1.0 2.0 3.0 4.0 ENERGY (RELATIVISTIC UNITS) IN MEV Figure 44. Effect of non-Gaussian Resolution Correction on Y9 Kurie plot- For energy in Msv see Appendix B, Figure 67.

-109been improved by a factor of about 2 over a flat-type detector. Even more important is the fact that electrons scattered out of a flat detector cannot easily be corrected for. This effect has been greatly minimized by using a hollow detector. Finally, when necessary, the still relatively poor resolution may be corrected for using a nonGaussian instrument profile scheme.

CHAPTER V DATA ANALYSIS: APPLICATION OF DIGITAL COMPUTERS Two principal types of problems encountered in this work require lengthy mathematical calculations. These are the Kurie analysis, discussed in Chapter II, and the resolution of complex radioactive decay curves. It, therefore, seemed advisable to make use of the computational facilities availabe at the University of Michigan. The IBM 650 Digital Computer was used. It is of medium size (2000 words of storage) and has a magnetic-drum type memory. Punched cards are used for input and output. Physically, there are three separate units that comprise the computer. These are: 1. The type 553 Read-Punch Unit. It is here that the instructions are fed into the machine, and where the results are punched out. Maximum input and output rates are 200 cards per minute and 100 cards per minute, respectively. 2. The type 650 Console Unit. This unit contains the magnetic drum on which information necessary to process a problem is stored. There is also a console from which information my be read into the memory, or displayed from the memory on the display lights. 3. The type 655 Power Unit. This contains all of the power supplies for the other units. It also translates the decimal input into the machine code, and the machine code into decimal output. The set of instructions used by the machine to compute the desired results is called a "program." There are certain routines available which facilitate the writing of programs; these will be re-110

-111ferred to under each specific program. In the program itself, one instruction is written per card. After the program is read into the machine, the data cards are then entered. The form of the data cards varies, and will be given for each program. After the computation is completed the results appear on punched cards and can then be printed out on the Reproducer. A. Kurie Analysis 1. Computer Program for Kurie Computation. The program was first written in symbolic notation in which no definite drum locations or "addresses" were used. This was then assembled and optimized by means of a program known as S.O.A.P., Symbolic Optimal Assembly Program. By means of this program addresses are assigned for the instructions in such a way that a minimum of time is required to perform the program. All data are punched on cards in eight groups or "words" each containing ten digits. The first word on each data card contains the address of the first piece of data on the card along with the total number of pieces of data on the card. This is the usual seven word format used when data are to be read in by means of the seven word loading routine. The Kurie analysis program is essentially a data handling program. Seven tables are stored in the computer. The experimental data are modified and correlated by means of extensive "table lookups," interpolation within the tables, and simple arithmetic calculatians. The following tables must be supplied, along with the instructions.

-1121. Table A contains the disintegration rate for each point in the B spectrum. A maximum of 51 points may be used, and the maximum rate is 99,999.99999 counts per unit time. 2. Table B contains the energy in Mev associated with each point. Again a maximum of 51 points may be used. The maximum energy permitted is 999.9999999 Mev. 3. Table C is the list of energies that serves as the argument for Table D. This table is the same for every element. 4. Table D is the list of values of the Fermi function as listed in Reference 99. A separate table must be made for each different element to be analyzed. There are on hand at present tables for S, P, Pm, Sr, Y, Rb, In, Cs, and Hf. 5. Tables E and F contain the information necessary to correct the Fermi function for the screening effect in the base of the B- emission. The data come from Table 6 of Reference 99. 6. Table G contains the energy argument for Tables E and F, and is the same for any value of Z. 7. V1 and Zl are two numbers on a single card, different for each element. V1 is the "outer screening potential" referred to before. It is obtained from the K binding energy data given in Figure 13 of Reference 72. Z1 is a number needed to interpolate in the E and F Tables. When the output cards are read by the Reproducer, the results of the Kurie analysis appear in two columns. The first column is the list of energy values, directly opposite which, in the second column, appear the appropriate ordinates for the Kurie plot. At present, forbidden spectra are not correct in the machine.

-113While the IBM 650 computation time for an average Kurie analysis of thirty points is less than one-half minute, the overall time required to analyze a spectrum may be about one-half hour since it is necessary to punch the data on input cards. A hand calculation of the problem might require a day or longer. Over sixty computer analyses have been run to data. Figure 45 is a flow diagram of the program. 2. Simple Least-Squares Line Fit A simple least-squares program was written primarily to be used to determine more accurately the endpoints of Kurie plots. However, the code is general purpose and may be used whenever a straight line is to be fitted to data. Some examples of its use are: analyzing one or two component decay curves, interpolation in and extrapolation of tables, and accurate calibration of linearly behaving equipment. The principle of least squares states that the best representative curve that can be fitted to a set of data is that curve for which the sum of the squares of the residuals is a minimum. Using the equation of a straight line: a = a + bt (20) the residual ri for the ith measurement would be: ri = ao + bti -ai (21) The least-squares program was coded in MITILAC, an interpretive routine which simplifies coding for small problems. The data are punched in the standard MITILAC data card format. The program was designed to be able to analyze ane y number of curves without having to read in the MITILAC deck and the instructions each time.

-114START Using B table elements interpolate in tables C and D to find Fermi Functions for each energy. Using B table elements, and Zl, double interpolate in table E, F, and G to find screening correction factors. Screening correction factors times V1 % corrections. Modify Fermi Functions by correction terms to get fcor' Counting rate r cor Counting rate 1/2 fcor Punch out Energy and Kurie Plot ordinate tables. STOP Figure 45. Flow diagram of Kurie plot calculation by the IBM 650 computer

-115The program fits a straight line to a maximum of 200 points in v 1 minute. The results of the computation appear on oe card. The first word is the ordinate intercept, the second word is the slope of the line, and the third word is the abscissa intercept. B. Decay Durve Analysis Accurate knowledge of decay constants is of real value both from a theoretical and also a practical point of view. The determination of cross-sections of nuclear reactions is a good example. The number of atoms produced in a given reaction is usually determined by following the decay of the product nuclei if they are unstable, and extrapolating the decay curve back to a "zero" time. This extrapolation presumes the knowledge of the half-life, and hence the value of the cross-section will be in error by at least the error in the half-life. Log (ft) values described in a previous chapter also require a knowledge of the half-life. Information of a fundamental nature in 3-decay theory may be obtained from knowledge of transition probabilities and half-lives of meta-stable states. The identification of components in a radioactive mixture may sometimes best be done by half-life measurements if short half-lives are involved. In any case, a half-life measurement may be used to confirm an identification based on other facts, or if determined first may limit the choice so that further confirmation may be facilitated. A second and equally as important problem is the determination of the true number of components that make up the decay curve. If the actual number of components is uncertain or arbitrarily determined by the experimenter, the concept of accurately determined half-lives has no meaning.

-116Since accurately known half-lives are valuable and since it is technically feasible and even convenient to obtain good decay curves, the answer to the question of why half-lives are not measured more accurately must lie in the actual analysis of decay curve data. The decay of radioactive nuclei is a first order process, and may be represented by a differential equation of the type: - dN = NX (22) dt where X is the decay constant and N is the number of atoms of a given type that are present at time t. Integration gives: N Noe-Xt (23) where No is the number of atoms present at "zero time." In terms of observed counting rate A, Equation (23) becomes A = Aoe t (24) If several different radioactive species are present, the observed counting rate h(t) is then n h(t) = Z (Ao)ieXi (25) i=l i = integers from 1 to n The problem may now be stated as follows: a function h(t) is approximated by experimentally determining an estimate of h(t) at a finite number of values of t. From this discontinuous set of data it is desired to obtain n (total number of components) and estimates of the (Ao)is and XIBs. An error analysis is also desired, based on the errors inherent in the method of solution, as well as in the values of h(t) and t. It should be noted that in Equation (25) the exponentials are all assumed to be separate and

-117unrelated. In the physical case this means that there are no decay chains involved; each radioactive parent decays into a stable daughter product. The essential difficulties in the solution of this problem are that we are dealing with a series of non-linear equations, and that the data are only approximating the function h(t) over a finite range of t. The importance of the cut-off of h(t) at a finite t depends upon the method of solution. Since there can be no analytic solution to the problem, a numerical approach must be taken. Applications of techniques in numerical analysis make possible an approximation to the solution. The computational labor involved, which would normally be prohibitive, is now no longer a problem due to the advent of high-speed digital computers. Several approaches to the problem of decay curve analysis will be reviewed, including a method that apparently has not been successfully applied to this problem before. At present the author knows of no truly satisfactory method for solving this problem that has actually been tried. Most methods do not make use of all of the accuracy inherent in the data, impose conditions an the data that cannot realistically be fulfilled, do not necessarily converge to the proper solution, or else are incapable of yielding an error analysis. It would seem that the solution proposed will in part overcome these difficulties. 1. A Review of Previous Methods a. Grapical Method. By far the most common method used to resolve a decay curve into its components is the graphical approach. Here one plots the counting rate as a function of time on semilog paper.

-118The last portion of the curve where only one component remains will appear as a series of points to which a straight line may be fitted (see Figure 46). This is arbitrarily done using a straight edge after determining by eye where the straight line segment ends and the upward curvature begins. The straight line segment is extended back to time "zero" and the subtraction or "peeling off" procedure is begun. The ordinate values from the extended line used in the subtraction are read from the graph as is the intercept at time zero. Normally, the halflife is read directly from the graph. The peeling off process is repeated until the entire curve has been analyzed. The limitations of the method are apparent. Each straight line is fitted with no better criteria than a visual impression. Human judgment rather than some mathematical treatment is used to determine the number of components. Fictitious components may be introduced or real components overlooked in this process. For details on the errors introduced by the peeling off procedure see the next section, titled "Least-Squares Method." In this graphical method there is no error analysis, and hence no criteria on which to base a judgment as to which combination of components out of many possible combinations is the best. Finally, values are plotted on and read from a graph, which in itself is a practical limit to accuracy. The method is, however, certainly the easiest to perform. b. Least-Squares Method. The graphical method may be considerably refined by employing a least-squares curve fitting procedure. This procedure my be applied only to equations linear in the coefficients,

-119or equations that may be reduced to a form linear in the coefficients (121). Consider the form of a single component decay curve written in terms of counting rates: A = A eXt (24) o This may be transformed to give lnA = lnA - Xt (26) which is linear in the coefficients In Ao and t. Two forms of the least-squares method have been programmed for the computer. The first is the simple, or unweighted, line fit previously described. It may be used to analyze a single component decay curve of the form shown in Equation (26). Even two component curves may be treated in the following way (7): Assume that the half-lives of the two components are known and it is desired to obtain the activities A1 and A2 of the two components at a certain time to called the "zero time." This is the case where nuclear reaction cross-sections are concerned. The total observed activity is h(O) = A + A2 (27) and at any time t the expression becomes h(t) = AleXlt + A2e 2t (28) where X1 and X2 are the decay constants and t is the time interval between zero time and the time of determining h(t). Two cases present themselves: Case I. Here the half-life of A2 is very long compared to Ab and A2 remains essentially constant over the period of the decay measurements.

-120Then eX2t 2 1 and the expression becomes h(t) = Ale + A2 (29) This expression is in the form of a straight line with h(t) as the ordinate values and eXt as the abscissa values. A least-squares determination of the slope of this line gives Al while the ordinate intercept gives A2. Case 2. Here the half-life of A2 is of the same order of magnitude as. Al. The expression then becomes: -Xlt h(t) Ale + A2 (30) e-2t e -X2t When h(t) is plotted against e -( 2-1) the slope is A1 and the e-X2t ordinate intercept is A2. Since the above equation only holds when there is a finite amount of Al activity present, a divergence from linearity will occur when the ratio of the activities becomes very small. In practice the curved portion should be neglected when the data are being analyzed. The second least-squares method involves the use of weighted ordinates. Recall Equation (26): InA InAo - Xt (26) Define nA = a, nAo = a, and -X = b. The equation now becomes a = ao + bt (equation of straight line) (20) The residual for a given measurement is defined as: ri = ao + bti - a (21)

-121If the measurements are not all of equal precision, or as in the present case the residuals are of a function of a variable rather than of the variable itself, then the residuals are not all of equal weight. In this discussion it will be assumed that the variable t is not subject to error, and that the measurements of A are of equal precision. However, since the straight line is being fitted to values of InA, the squares of the residuals must be weighted by the quantity win and for lnA functions wi is proportional to A (see Scarborough, Reference 121). Hence, it is desired that: wir = Z Ai(ao + bt - a)2 = minimum (31) Taking partial derivatives with respect to ao and b and setting them both to zero the normal equations are found: ao A2 + bZ A2tZ = Aa (32) a Z A2t + bZ (A2t2) = E A2ta (33) These yield a = Z A2aZ A2t2 - Z A2at A2t = A2aZ A2t2 - Z A2atZ A2t (34) EA2 A2t (A A2t)2 and b = Z A2Z A2at - Z A2tZ A2a = LAT (35) Z nA2Z -2 ( A2t)2 LT EATAt (Z( t) Certain combinations of sums are defined as LT, LAT, and LA for convenience in writing. The sum of the squares of the deviations of the points from the fitted line divided by (n - 2), where n is the number of points, is

-122a measure of the scatter of the points and is called S2 or the "variance of the estimate" of a single ai measurement (137, 150, 151). S2 may be calculated from the following equation: S2 [ZA2Z A2t - (ZA2t) 2] [AA22a2 - (ZA2a)2] -[ZA2ZA2ta -(A2t ) (Za)] (n - 2)( A2)[ A22t2 - ( A2t)2] (36) or in abbreviated notation: s2 = [LT[LA] - [LAT2 (37) (n - 2)( I A)LT The standard deviation of the estimate ae = S2The other standard deviations of interest are: Intercept: ao = A2 (38) LT Slope: b = a /ZA2 (39) LT The following procedure was used in the IBM 650 program for the resolution of complex decay curves. Consider the two component decay curve shown in Figure 46. Let the total number of points on the curve equal n, and the number of points on the "tail" of the longestlived component equal p. First the p points, each with coordinates Ai and ti are read into the computer in order beginning with the highest t value, one point per card. From these points the slope, intercept, ae, aa, and ab are calculated. The quantities (ao + ca ) = a+ and (ao - a ) = ao are then formed. Finally, the values of the slope, ae, O ~o a o

-123n-p-I POINT o (ao)2\ II: < \o0 p= POINTS ON"TAIL" (ao)l — t Figure 46. Two component decay curve a+ + the intercept AoJ and the range limits of the intercept e o = Ao and eao = Ao are punched out. The next data card containing the coordinates for the point n-p-l is then read. Call the coordinates Ar and tr. The value of a at tr for the previous component is calculated and the quantity es = A is found. A is then subtracted from Ar giving A' or the value of r

-124A of the second component at tr. Finally the values of A' and t are r r punched out. The next data card is read in and the procedure repeated until all data cards have been read. The output cards are punched out in such form that they can be used for input cards, when the program is rerun for the next component. An automatic programming routine called S.I.R., Symbolic Interpretive Routine, was used to program the problem. It is similar to S.O.A.P., with the addition of floating point arithmetic (automatic scaling) and several additional operations such as SINE, EXP, LOG, etc. The program can handle up to approximately 300 points on the experimental curve at one time. Errors. Actually, the weight of each residual should be w = A2A where wA is the weight of the ordinate A due to the uncertainty in the determination of A. This uncertainty may be due just to the statistical nature of the counting process, or it may also contain experimental errors such as variations in the sensitivity of the counter or in counter geometry. Experimental errors of the above class are usually impossible to describe mathematically, and so the best one can do is to try to minimize these errors and to leave them out of the mathematical treatment. Since the statistical errors in the counting rates usually follow a Poisson distribution (Xt<< 1, N >>l, and (At)<< No, where No = number of atoms at t = 0 and t = time of the count) the e standard deviation in the counting rate is just 0A =1T When these deviations vary over a large range from point to point they should be taken into

-125account when the residuals are weighted. This may be done by calling WA = 1 (40) (A)2 Then the total weight is: w=A2 (41) \aAj The factor oA was not included in this scheme because in a good decay curve the value of aA is of the order of 1 or 2% and does not vary greatly from point to point. Furthermore, the inclusion of aA would complicate the computer program and place much greater limitations on the size of the problem that can be handledc Finally the results obtained using GA after the first component has been subtracted are subject to criticism from a mathematical point of view. The omission of A has two important effects. First, in dealing with short halfA lives where the decay is followed automatically by recording the number of decays that occur in a fixed time interval, the aA becomes progressively larger. If all points are considered to have equal precision the standard deviations of the fitted curve will be larger than necessary. However, in just such a case the error in the time interval may become important and may not justifiably be ignored. Therefore, the least-squares procedure employed here would no longer be truly valid, and a new one would be required which took into account the errors in both the ordinate and the abscissa (53). It should be remembered that with short halflives where the time of a single count is say 1% or more of the halflife a correction must be made for the decay that occurs during a count.

-126If N atoms are present at the beginning of a count lasting for time t, then the number of disintegrations that occur will be WN. = No -No e-Xt = No(1 -e) (42) If Ao = d = NoX, (43) dt then A = N (44) -Xt 1 - e Thus, each count may be corrected to find the disintegration rate at the beginning of the count. The second effect caused by the omission of aA is to further emphasize what is probably the chief disadvantage of this resolution scheme. The "peeling off" process, that is, the procedure of subtracting off the line fitted to the longest-lived component, automatically throws all of the error due to the scatter of the gross points into the remaining curve of the shorter-lived components. This is inherent in any peeling process, be it graphical or mathematical. Assuming no errors of the experimental type, the scatter of the gross curve is a function only of the actual number of disintegrations recorded at each point. Consider a point on a two-component curve where each component contributes to the gross point. Let us say that the counting rate and standard deviation is 10,000 + 100 counts per minute. Assume that the extrapolation of the longer-lived component gives a value of 9,000 c/m at this point. Upon subtraction the shorter component would have a value of 1,000 c/m. Yet the absolute value of the scatter of the points

-127is the same for the second corpanent as it was for the gross curve. Percentwise, the error has increased from 1% to 10%. Merely because a subtraction has taken place, the uncertainty in the second component has increased. While no smoothing of the original data based on an assumed knowledge of the error in the longer-lived component is justified, the following procedure may be used to get a better estimate of the second component. The standard deviation of each point in the gross curve is calculated and carried along with each point. After the first component is subtracted the numerical standard deviation of the points on the gross curve are given to the points on the remaining curve at the same abscissa values. These are then converted to percent errors and squared to get the variances. The value of 2 for the first curve is added to the above variances to give a new Cl. The weight of the residual for the second curve would be as before: \.A1/ W ()2/ (45) Using this weighting scheme a somewhat better line and a smaller error are obtained. Two things must be remembered: first, the weight factor must be carried to as many decimal places as the rest of the data; and second, the weighting procedure just described is not mathematically rigorous. The latter is due to the fact that the error in the second component was calculated using only those points on the gross curve where the contribution of the shorter-lived component was essentially zero. The information present in the first portion of the curve was not used to determine the longer-lived component. Furthermore, the errors in the

-128original data are homogeneous, thereby producing skewed errors in the logarithm form. The above procedure has assumed homogeneous errors in the logarithm form, which leads to a skewing of the errors in the final result. In sumnary, it may be said that there are several disadvantages to a least-squares "peeling" process. Every time a subtraction takes place the scatter of the original points automatically falls into the lower component. This is because the residuals of the entire curve are not made a minimum. Therefore, the shorter components cannot be determined with the accuracy that is inherent in the original data. The number of components has to be decided upon, it does not fall out of the method. This usually required a prior rough graphical resolution. All but the shortest-lived component are determined using fewer points than are available in the original data. This again introduces error. Also half-lives that are very close together may not be resolvable. The reasons for using the method are as follows: better results are obtained compared to the graphical method, and an error estimation is available. The method is simple to program and the problem is of such a size that it could be put on a medium size digital computer, such as the IBM 650, which was immediately available. A possible modification of the above resolution scheme that would place more of the burden of analysis on the computer would be to allow the machine to decide how many points should be used to fit a line to the long-lived tail. After this has been decided, the subtraction could take place, and the desired information pertaining to the first curve punched out. Next the machine could automatically begin to fit

-129a line to the next component, and so on until the entire curve was analyzed. This would mean no prior graphical resolution would be necessary, although it might be necessary to tell the machine how many components make up the curve. The procedure for accepting points might be as follows: Arbitrarily select a minimum number of.points that will surely fall on each component's tail; this number will depend on the nature of the curve and how many points were taken. Let us say that the minimum number of points is 3. Have the machine fit a line to these three points, and then calculate the variance of the estimate. Let the next point be accepted, fit a new line to the now 4 points, and calculate a new variance of estimate. If the new variance is less than the first variance, the fourth point was good. The process is repeated until a variance is found that is larger than the preceding one. The last point would then represent a contribution from the second component and should be rejected. This procedure presumes reasonably smooth data, and might break down if the gross points scattered greatly. c. Prony's Method. A rather elegant method of fitting exponentials was suggested by A. L. Prony (58). As usual it is desired to determine an approximation of the form h(t) l alelt + 2e 2t+... + aX e-nt (46) 2 n which is equivalent to t t t h(t) 0 a.2 + 0 +... + %0n (47) where X ei (48)

-130Prony's method requires that h(t) be determined at N equally spaced points, t = 0, 1, 2,..., N-l. This restriction is a distinct disadvantage, because decay data are usually not taken in this manner. However, this restriction is common to all of the following methods, and may be attained by an interpolation procedure of whatever degree of accuracy is Justified by the original data. Once equally spaced points are determined a simple linear transformation of the variable t will yield the required integral values of t. Now if Equation (47) were actually an equality, we would have 1 +2 +.. +n = h(0) pl 1 + a2>1.2 + *@ + Afh = h(l) N -1 N -1 N-l1 %aL +a'2 + +..+ 0N h(N - 1) If the constants C!,2'...,,' n were known, this set would comprise N linear equations in n unknowns, Ca1, aC2,, o, and could be solved exactly if N = n or approximately using the least-squares method if N >n. However, since the I's are also to be determined, a minimum of 2n equations are necessary. The set of Equations (49) is still not linear in the p's. Therefore, let P1, P2,..., in be the roots of the equation,,n _ Bn-l 1_,1n-2 - P...2 - = t -. (50)

-131In order to find the 3 coefficients, the first equation in (49) is multiplied by A, the second by Pn.1, and so on. Adding the results h(n) - lh(n - 1) -... - h(0) = 0 (51) is obtained. In a similar way the N-n-i equations are obtained h(n - 1)1 + h(n - 2)32 +... + h(0)Pn = h(n) h(n)p1 + h(n - 1)02 +.. + h(1)Pn = h(n + 1) (52) h(N - 2)p1 + h(N - 3)t2 +... + h(N - n- 1)n = h(N -1) The set of equations (52) may be solved directly to give the Its if N = 2n, or approximately using least squares, if N >2n. Knowing the's the roots p1i, 2l'.. pln of Equation (50) may be found. The p. values may be used directly to calculate the X's from Equation (48) and together with the set of equations (49) may be used in a leastsquares method to obtain the a values. While Prony's method gives the appearance of being a leastsquares fit, it is actually not since both the 3's and the h(t)'s are subject to error. Hence no error estimate is available. Also, there is no provision for weighting the values of h(t) according to their precision. Another serious drawback is that the number of components must be known. Finally, for large numbers of data, the method becomes cumbersome computationally. d. Householder's Method. A. S. Householder in 1949 proposed a modification of Prony's method (62) which permitted valid least-square estimates to be obtained. Essentially he applies Prony's

-132method to obtain initial estimates of the exponents, and then by an iterative procedure arrives at least squares estimates. A test is also included to determine how many exponentials are needed to adequately fit the data. The method was programmed for the ORACLE, a computer at Oak Ridge National Laboratory, and for a reason that is not known, the method sometimes failed to converge, and even occasionally converged to unreasonable estimates (61). e. Cornell's Method. Another procedure for the estimation of a linear combination of exponentials was put forth by R. G. Cornell in August, 1956, who tested the method on the Oak Ridge National Laboratory's computer, the ORACLE (28). In some respects the method appears to be similar to Prony's method, except that among other things it may also be applied to group means instead of only individual observations. This is quite important in investigations where an observation at a given point may be repeated. Furthermore, Cornell's method does not involve questionable least squares calculations. In its unmodified form this new method requires not only that equally spaced points be available, but also that the number of components be known and that the number of points N be equal to 2rn, where n is the number of components and r is any integer, if an error estimate is to be made. Ease of computation is comparable to Householder's method or Prony's method. Actually only the case of a single component exponential was investigated to any great extent. Here the results seem to be at least as good as most other methods provide, and better in many cases, particularly when the half-life was long and smooth data were available.

-133Since the method is still new and relatively untested a further study must be made before it can be accurately evaluated. 2. Method of Fourier Transforms. The previously listed methods all contain certain undesirable limitations, and it was felt that a search for a new method was warranted. There is little that is original in the mathematics of the method to be presented. In the main, it consists of a grouping of previously known techniques into a program which appears to have several distinct advantages over other methods. As in the case of Cornell's method, the procedure now to be presented has not been thoroughly tested, and until this is done its actual value is difficult to assess. One of the mpst important reasons for desiring a new method is that none of the previous methods, except Householder's, can estimate or determine the actual number of exponentials that comprise the gross curve. It was necessary to rely on human judgment to obtain this information. It seems to this author that it is just as important to be able to accurately estimate the number of components as it is to extimate their constants. In fact, if the number of components is uncertain then any accurate calculation of the constants involved begins to lose its meaning. This would not necessarily be true if it were only required to produce an arbitrary function that would simulate the data. But the constants involved have physical significance and their estimation is fundamental to the entire analysis. A satisfactory analytical procedure must also include the following features: 1. An error estimation. Without this there is no way of evaluating the results.

-1342. Full use of the accuracy inherent in the data. The lack of this feature was pointed out in the graphical and least-squares peeling process. 3. The replacement of human judgment by mathematical analysis. This is particularly needed in the case of the number of components. 4. Treatment of all observations as a whole. This is essential if a minimum error is to be obtained over the entire curve. (See also remarks in Least-Squares Method section) a. Development of Initial Problem into the Form of a Laplace Integral Equation. Consider first the case where nuclei of type A decays into stable nuclei of type B. Also let A = A(t) and B = B(t) be the numbers of atoms of types A and B present at time t. Assume that B(O) = 0 At time t = 0 A = A B = 0 At time t =oo A =0 B B = A At time t = T A = Aoe-AT B = Ao A = A( - e-XT) The rate of decay of A is equal to the rate of growth of B dA = dB dt dt (53) Also, dA= AA =AoxAeXAt (54) dt

-135The actual counting rate obtained experimentally may or may not be the true time derivative of A. It is assumed that the observed counting rate is proportional to the true disintegration rate. Combining (53) and (54), we may arrive at A B t - dA dt = f dB dt = A e A dt (55) A dt o dt o 0 Integrating we obtain B = - [Ae-At] = - AoeAt + Ao (56) Rearranging gives Ao B =Ae A (57) Let f(t) = Ao - B(t). Then f(t) Aoe-XAt (58) In a sample containing several radioactive substances decaying independently into stable daughters, Equation (58) can be generalized to give f(t) = A eXAit (59) i i The following conditiins are required: Aoi O, i > 0, and Xr A Xs when r s. Equation (59) is in the form of a Dirichlet series (148) and may be expressed in integral form as: ) ()e d (60) o

-136It should be noted that a plot of f(t) vs. t is actually a plot of the sum of the decay curves of the radioactive parents on an atom basis. Therefore, in the case of a single component the f(t) is actually the number of atoms of A present at time (t). The function f(t) may be obtained by integrating the disintegration rates from zero to t, for each t, and plotting the results against t to obtain B(t). This curve will approach a limit as t -~ oo of Z Ao. Each point on the curve is i i subtracted from the limiting value to give f(t). Equation (60) has the form of a Laplace Integral equation. b. Solution of Integral Equations by the Method of Transforms. The process of multiplying a function f(t) by a known function k(u, t) and integrating with respect to t from a to b to generate a new function F(u) is a type of transformation whose name depends on the kernel k(u, t). A Fourier transform is an integral equation of the form, F(u) / f(t)e dt (61) -00 iut where e is the kernel k(u, t). A generalized linear integral equation of the "first type" may be written as b f(t) = I k(X, t)g(X)dX (62) a where g(X) is the unknown function, k(X, t) is a kernel of known form, and f(t) is a known function or a tabulated set of points. In most cases, an analytic solution does not exist and few numerical methods have been worked out. The method (115) that will be used is limited to cases where k(X, t) is of the form k(t - X), k(t + X), or k(Xt).

-137A simple change of variables will transform each of the above into the form k(t - X). First, 00 f(t) = f k(t - X)g(X)dX (63) -00 Now iut F(u) = 1 f(t)e dt (64) where F(u) is the Fourier transform of f(t) and u is a dummy variable. Substituting for f(t) in Equation (64) the integral in (63) we obtain 00 00 F(u) = 1 f k(t - x)g(X)dX f eutdt J2-J and F(u) = 1 g(X)dX f k(t - X)eiutdt (65) "' " 00 - o0 Let t = x +. Then x = t - X. Therefore: F(u) = 1 Lg(X)dX k(x)eiU(x+X)dx (66) Since eiu(x+X) = eiux eiu, we have ao00 oo It is seen from Equation (61) that the first part of Equation (67) is the Fourier transform G(u) of g(X), and the second part is the Fourier transform K(u) dc k(x). Therefore, Flu) = ~2 G(u)K(u), (68)

-138and G(u) = 1 F(u) (69) Nfi K(u) Taking the inverse Fourier transform of G(u) we obtain the original function g(X) which is now a function of the transforms of f(t) and k(t, X). 00 g(X) = 1 f F(u) e-iUXdu (70) 2-1 K uJ The results obtained in Equation (70) are purely formal. For a precise statement of the conditions under which Equation (70) holds, see References 139, 141, and 115. c. Solution of the Laplace Integral Equation. We begin with Equation (60). f(t) = etg()dX60) (60) and proceed to transform the variables X and t. Define= e t = e (71). Then -e(Y) f(ex) = - g(e-Y)e-Ydy (72) -00 Multiply both sides of (72) by eX: eXf(ex) =- L e (x)g(e-y)e (- )dy (73)

-139Define f(e ) = (x), g(e ) (y) (74) Then -00 eX(x) = - y ee (y)dy (75) The Fourier transform of the left hand side of Equation (75) is F(u) =- 1 J e f(x)eiUdx (76) Let s = x - y, and forming the Fourier transform of the kernel of Equation (75) K(u) =- 1 I e e Seiusds (7) fir Finally we form g(y) = 1 f F(u) e-idu (78) In this case, K(u) can be evaluated analytically. In Equation (77) let es = Z, dZ = eSds, ds = e-SdZ = dZ/Z. Also let iu = n - 1. Therefore, n = 1+ iu (79) Substituting into Equation (7 ) we' get: K(u) = f e-ZZn-ldZ (80) -00o

-140which is by definition (65) the rfunction, r(1 + iu). Therefore, the kernel K(u) turns out to be the Euler integral for 1 (1 + iu). d. Details on the Method of Solution. Briefly, the method of solution was shown to consist of first integrating the initial disintegration rate curve and then forming a new decay curve on an atom basis. This was called f(t). Next, the Fourier transform F(u) of f(t) was needed. F(u) = 1 I eX(x)eiUdx (76) K(u) was shown to be the I function K(u) = 1 r(l + iu) (81) Next, the quotient F(u)/K(u) was formed, and the inverse Fourier transform taken to give g(y) where'(y) = g(e-Y) = g(X). 00 (y) = 1 f F(u) e-iyUdu (78) 2 " K(u) Basically then, the method requires two integrations and an evaluation of a complex I function. The details of the numerical integration will be given later, but it is clear that we cannot integrate numerically from - oo to o in Equations (76) and (78). [Equation (76) will be discussed here, but the remarks apply also to Equation (78)] Therefore, we must introduce the limits xo and -xo into Equation (76). These are the cut-off points of the integral, and naturally we would like to

-141have ix o approach co as closely as possible. Instead of Equation (76) we now have: I F(u) = f ex (x)eiUdx + E(x, u) (82) -Xo The cut-off at x = + x and the subsequent dropping of E(x, u) is equivalent to folding the transform of F(u) against the transform of the unit step function: O, x > x m(x) = 1, -x i x (83) O, x <-xo Consider the following sketch: f (x) xo x Figure. Sketch of (x) as a function of x Figure 47. Sketch of T(x) as a function of x

-142There is no difficulty in choosing a zero on the x axis so that the necessary conditions of symmetry are maintained. However, it is necessary that a finite xo exist such that the value of E(xo, u) is sufficiently small to render a good solution possible. Since f(x) is obtained from experimentally measured quantities, knowledge out at least as far as x is required. If this knowledge cannot conveniently be obtained experimentally, it must be obtained in some other way. More will be said about this later. The major difficulty in this method is now apparent. We are trying to simulate a curve with an abrupt tut-off at x by a sum of exponentials that extend to x = o. This process introduces high frequency Fourier components into F(u) which will tend to obscure the results. In fact, if E(xo, u) is large enough, it will be impossible to obtain a good solution. It is most unfortunate that once the frequency spectrum of f(x) has been warped by folding against the step function m(x), which may be represented by sin(ux)/ux, the warping cannot be removed from g(y) by any unfolding. It will help in the discussion of the cut-off error if the form of the results is examined. g(X) may be thought of as a aum of A functions. However, due to errors in the initial data, the integration process, the evaluation of the r function, and in the numerical calculations due to round-off and truncation of numbers, the best that can be hoped for, when cut-off errors are not considered, is a frequency spectrum similar to the one shown in the following sketch which

-143represents a hypothetical four component curve. X Figure 48. Sketch of g(X) as a function of X Each peak in the spectrum indicates a component. The abscissa value at the center of a peak is the decay constant, while the height of the peak gives the number of atoms of the component at time zero. An error analysis would come from the profile of the peak. The beauty of the method lies in the fact that the number of components auto

-144matically falls out of the analysis. How well the components may be resolved naturally depends on how good the initial data are, and also on the errors introduced during the numerical calculation. The region of the sketch around X2 and X3 illustrates the case where two half-lives are close together. A better resolution may be effected by starting at one end of the total curve, fitting a distribution curve to the first peak, and subtracting it from the total curve. In this way, each component may be "peeled" off of the gross curve. Returning to the discussion of the cut-off error, it is of interest to see how this factor will affect the overall results. It was noted that the cut-off at a finite value of x has the effect of adding into F(u) Fourier components which extend the range in u on which F(u) maintains appreciable value. In order that the integral in Equation (78) converge, F(u) in the quotient F(u)/K(u) must approach zero at infinity faster than K(u). Now, +(l +iu)j { 1 (84) T4T sinh(Ju) and so K(u) tends to diminish quite rapidly. Therefore, when F(u) is warped by the cut-off (or by other sources of error), for some value of u the quotient F(u)/K(u) will begin to grow without bound. The result of the extra Fourier components in F(u) is not only to smudge out the peaks but also to cause them to be riding on a function that is tending to increase in value. The sketch below illustrates this case for a hypothetical two component decay curve. Even this

-145poor solution (and it could be far worse) does yield some information. g (X) XI ~I Al X2 Figure 49. Sketch of g(X) showing effect of cut-off error The X's may be determined with some accuracy, and an estimate of g(Xl) and g(X2) may be obtained if the questionable procedure of trying to subtract off the contribution of the rising function is employed. If nothing is done to counteract the effect of the cut-off, at least something may be said about error formulae for the cut-off.

-14bThere are two cut-off integrals to be evaluated: 00 (a) feX?(-x)elLidx E (a) e-X (-x)e - =E, x o 0 and (85) (b) fIeX1(x)elxdx = E x o Perlis in Reference 115 states the necessary conditions for the behavior of f(x) at xo, and obtains the following results' (a) E- _ o(e (n)o) In - 11 - lu (86) (b) E+ ~ O(e(n+l)xo) in + 1i - lul where n L 1 in the case (a) and n L -1 in case (b). He obtains a more specific result for the cut-off in forming K(u). Range Error x too -o [(l + iu) + (1+ iu) (87) oto -rl (1+ iu) xo where 17 refers to the incomplete P function of x with integration limit j. As x~ - oo, the above errors -, 0. Thus, in both F(u) and K(u), the errors should be quite small at any reasonably large value of x. There are two principal approaches to the problem of cut-off errors. The first and less satisfactory approach is to introduce a convergence factor into Equation (78). Through a sacrifice in

-147"resolution' of g(y) a convergent numerical solution is available. au2 Two convenient factors given in Reference 115 are e and sin2au/(au)2. Therefore, instead of calculating g(y), we find the quantity \* 00 -iyu g (y) = 1 F(u) e T(u)du 2 ( 1 + iu) (88) 00 ='(y - s)t(s)ds -co where T(u) may be e-au for example. Then Equation (88) becomes g(y) = 1 g(y - s)e ( ds (89) -00 \/2a The solution is seen to be folded against the transform of T(u). The purpose of this is to have the transform of the convergence factor approximate the actions of a unit step function. That is, we would like to have the solution multiplied by a factor that was essentially one up to the point where the solution started to rise without bound. Then the convergence factor should cause the solution to fall to zero. How rapidly the convergence factor falls to zero may be controlled by the constant a. The greater the value in u at which F(u)/K(u) remains well behaved, the closer to zero the constant a may be made. The smaller the value of a, the narrower the spread by which each point in g*(y) will be weighted. This solution to the cut-off problem suggested in Reference 115, is undesirable from several respects. First of all, it introduces into'*(y) oscillations of decreasing amplitude which tend to obscure the results. Secondly, it can at best only try to subtract out the "background noise" on which the peaks are riding and keep the range of

-148the number size from becoming overly large. Actually, the subtraction might be done graphically. Finally, the g*(y) produced cannot be unfolded from the transform of T(u), and the g*(Xi) found are only estimates of the desired g(Xi) being in error by an unknown amount. A second, and it is believed better, approach to the cut-off problem is that of artificially extending the observed values of f(t) to the point where the cut-off becomes relatively insignificant. We shall consider two cases, in each of which it is assumed that the half-lives of the components are short enough so that it is convenient to follow the gross decay until only the longest lived component remains. Case 1. Here we require also that there is not a vast difference in half-lives between the shortest- and the longest-lived components. The numerical calculations are to be carried out on a computer, and so the maximum difference in half-lives will be a function of the size of the computer. The procedure in this case would be to follow the decay of the longest-lived component until its decay constant can be determined with sufficient accuracy using a least squares fit. Similarly, its ordinate intercept at time zero is determined. This information will then allow a simple extrapolation to give values of f(t) at appropriate intervals out to where the experimental counting rate becomes comparable with the background. A cut-off at this point will s'till introduce some error, but this error should be almost negligible and g(X) should be very well behaved over the range of interest. In the integration scheme to be presented later, equally spaced points are required. It is necessary that small enough in

-149tervals be taken so that no appreciable error is introduced when integrating under the shorter-lived components. This implies that if the interval is small enough to suit the shorter-lived components then a great many values of f(t) that define only the longest-lived component will be required. This is a failing of most numerical integration schemes, but may be circumvented to a degree by using an unequal interval integration scheme. It becomes apparent then that if the distribution of half-lives is sufficiently poor so that several hundred values of f(t) are required (obtained by interpolation), then a fairly large computer like the IBM 704 is necessary to handle the problem. It should be remembered that K(u) must be evaluated at the same points as F(u). Finally, g(X) will also be evaluated over a large set of values. Case 2. This is the case where the half-life of the longest-lived component differs greatly in magnitude from the next longest-lived activity. For example, a three-component system might have the following half-lives: 1 hour, 10 hours, and 200 days. In 60 hours the 10 hour half-life will have been reduced by a factor of 64 while the 200 day activity will still be above 99% of its value at time zero. Therefore, to a good approximation we may subtract off the 200 day activity as a constant factor from the gross experimental decay curve. The remaining two components may then be analyzed much more accurately than if the 200 day component had been extrapolated down to background and the entire system treated. If this system had been followed for 200 days and then cut off, even with a convergence factor the results would probably not have been very good.

e. The Numerical Integration Scheme. The numerical evaluation of an integral such as b I = f f(p)cos(xp)dp (90) a presents considerable difficulty due to the rapid oscillations of the function cos(xp). An ordinary integration scheme, such as Simpson's rule, would require excessively small intervals of integration in order to approach the accuracy of a method specifically designed to handle trigonometric functions. L. N. G. Filon (42) has developed an excellent method for the evaluation of such integrals that requires integration intervals only as small as would be necessary to evaluate the integral without the trigonometric factor. Since the method does not seem to be well known it will be developed very briefly here. Let the range of integration be divided into 2n equal parts with an interval h, so that b a + 2nb (91) The following notation will be used: xh =, a + sh = s (92) f(a + sh) = fs where s = an integer. It is assumed that over the range (ps - h, s + h), that is, (PS-1l PS+1) the function f(p) can be approximated sufficiently well by the quadrature formula, f(p) = A + B(p -ps) + C(p - ps)2 (93)

Differentiation leads to f' = (3fs+ + +4fS)/2h s+l1 (94)'1 = (4f s - 3fs-l)/2h If I = s+ f(p)cos(xp)dp (95) Ps -1 Then xIs = f(p) - 2Cx- sin(xp) + f'(p)x-lcos(xp) (96) s -1 This reduce to XOIs = (fs+l - fs-l)(cos9 - sing)sin(xps) + (fs+l + fs-1)(@sinQ - 2Q-lsing + 2cosQ)cos(xps) + 4f (Q-lsing - cosQ)cos(xp ) (97) If we write o3a = 92 + GsinQcosg - 2 sin2g 93 = 2l(1 + cos2) - 2singcosQ] 93y = 4sing - Gcos9] (98) Then Equation (97) can be written as Is h[c(fs+lSin(XPs+i) fslsin(xPs1)) + 1/2(fs+lcos(xps+l) + fs cos(xpl ))+ S (XPs (99) If Is is summed for S = 1, 3, 5,..., 2n-1, we have the expression a f(p)cos(xp)dp = haf()sin(Cb) - f(a)sin(a)} + C2s + 7C2s a'^^10

-152where C2s denotes the sum of all even ordinates of the curve y = f(p)cos(xp) between a and b inclusive minus one-half of the first and last ordinates, and C2s denotes the sum of all of the odd ordinates. The quantities C, 3, and y, given in terms of 9 = xh, are determined by the set of equations (98) but for reasonable accuracy at small values of 0 it is necessary to expand the trigonometric terms as follows a = 203 - 2g5 + 207 -... 45 315 4725 = 2 + 22 - 404 + 26... (101) 3 15 105 v6T y = 4 - 22 + 4 _ 6 +... 3 15 210 11,340 When -0 O, then - O- 0, -, 2/3, and y -> 4/3. This results in Equation (100) becoming b f f(p)dp = (h/)[2C2s + 4C2s-1] (102) a which is Simpson's Rule. Values of a, 3, and y have been tabulated for small angles in radians (Reference 141) and in degrees (Reference 42). In a similar way, the formula for the integral containing the sine term can be deduced. b f f(p)sin(x dp = h[-a f(b)cos(xb) - f(a)cos(xa) (103) a + PS2s + YS2s-1] where S2s and S2s_1 have the same meanings as C2s and C2s_1 except that they pertain to the curve y = f(p)sin(xp). An estimation of the error may be made by first evaluating the integral using an interval h, and then halving the interval and repeating the process. Clearly if halving

-153the intervals makes no appreciable difference in the results, the process cannot be very far from the limit. An actual error estimation based on this halving procedure is given in Reference 42. f. Details on the Numerical Solution. The size of this problem necessitates a fairly large computer, something in the range of an IBM 704. It is conceivable that an IBM 650 could be used, but this would at least require that the problem be broken up into several small sections. Even then the number of points would impose a serious problem. It is not the purpose of this section to present a complete computer program. It is rather to stress a few details in the numerical solution so that when a large enough computer becomes available the programming will be facilitated. To begin with, the original data will have to be treated to minimize the cut-off error as was mentioned in section d. If it only requires the subtraction of a constant term this might best be done by hand. If an extrapolation is needed, a separate program such as the weighted least squares routine described previously will give the necessary constants enabling the extrapolation to be performed or else an interpolation program (which is also needed later) may be used to interpolate on the "tail.' Preferably then the original countingrate data will have been corrected for background and cut-off before the principal program begins. In forming f(t) it is first necessary to integrate to form the growth curve of the daughter products. Simpson's rule will serve as the integration scheme. The limit of the growth curve may be taken as the ordinate at tax. F(t) is then formed by subtracting the growth

-154curve from its limit and plotting the results against t. All these preliminary calculations are carried out automatically in the computer. We make the change of variable t = eX, and f(t) = f(eX) = f(x). The function f(t) ranges from 0 to tmax while f(x) now ranges from - o to in tmax = xmax. The point at t = - o is dropped. The remaining points are predominantly larger than x = 0. Therefore, since we have to integrate between - Xmax, it would be best to choose a new zero so that the experimental points are equally distributed on either side, primarily to make use of the symmetric properties of the sines and cosines used later. Each f(x) must be multiplied by ex and then the following integral may be set up F(u) 1 f max [f*(x) + f*(-x)]e dx (104) or more conveniently, F(u) 1 fmax [f*(x) + f*(-x)lcos(ux) + i[f*(x) - f*(-x)sin(ux)Jdx 42itw^~~~~~~~ ( ~(105) where f*(x) = eXf(x). Essentially then Equation (105) and K(u) must be evaluated over the same range of u values. The Ffunction K(l + iu) - (Z) may be evaluated from a formula such as 1 ==Ze 1 + Z)e - (106) pT ~~n-il n J

-155where y is Euler's constant (y = 0.57721566...). It is also tabulated in the form of lnr(x + iy) = U + iV in Reference 100 for the ranges x: 0 to 10 y: 0 to 10 Therefore F(Z) = e(cosV + i sinV) (107) The range of the table may be doubled by use of a duplication formula given in Reference 100 and further extended by means of extrapolation. By using Equation (106) to calculate first r(iu) and then obtaining r(1 + iu) from (i + z) = z() (108) the computation is facilitated. F(u) and K(u) are complex numbers. The real part of F(u), termed F comes from the cosine terms in the integration, and the imaginary part Fs comes from the sine terms. Similarly, K(u) is composed of Kc and Ks. The quotient of F(u) and K(u) ist G(u) = F(u) = F + iF = (F + iFs)(Kc - iK) (109) Ku) K + iK K2 +K Since co 2tg(y) = f F(u) e-Udu, (110) we get If F(u) e-iydu = +iFK + leK, s+ FsKs (cos yu - i sin yu)du I K + K J(11) (III)

-156Making use of the properties of harmonic functions, Real Imaginary G(u) = G(-u) G(u) = -G(-u) all odd terms vanish in the integration (86), so that t(y) = T cF + FKss cos(yu)+ sKc - KcKs sin(yu du tK ~ 2 2 2 Kc c +Kc s (112) It should be noticed that in Equation (112) all imaginary quantities have vanished. Actually the integral in Equation (112) has the limits 0 to u, and the integration takes place over the same values of u used to max determine F(u) and K(u) or an equidistant subset of u's. Finally then, a set of values of g(y) are obtained over the range 0 to Ymax From these, the plot of g(X) vs. X is easily obtained. Some prior knowledge about the probable range of X's would be, of course, of use in determining what the limits on y should be. g. Summary and Conclusions. A method for the analysis of decay curve type problems, based on the Fourier transform has been proposed. The results appear as a frequency spectrum of g(X) vs. X. Fundamental to the entire analysis is the assumption that, while in theory g(X) could be represented by a discontinuous set of A functions, it is impossible to measure g(X) exactly by an experimental means and so, due to the inherent error in the data, g(X) lapses into a continuous function. However, it is also true that the resolution of the peaks in g(X) may be made as fine as desired merely by extracting an increasing number of values of f(x) from the initial data, and refining the integration scheme accordingly.

-157An error analysis is available for the integration scheme and also for the cut-off error. The error in the original data plus the error in the numerical calculations due to rounding off numbers, etc., will produce peak profiles showing a normal distribution as long as the cut-off error is very small and the integration scheme employs small enough intervals. Therefore, a standard deviation may be obtained from the peak itself. Under favorable conditions, the standard deviation should not change greatly from peak to peak in a multicomponent system, and hence may be used as a test to see if a peak might contain more than one component for X's extremely close together. That the number of components present may be determined by counting the number of peaks in the results is one of the best features of this method. A fairly large computer will be required to handle the problem if the maximum resolution inherent in the method is to be approached. However, the computer program itself should not be logically complex. Most of the subroutines needed, for example the determination of sine, cosine, exponential, logarithm, and routines for simple integration and interpolation, are available at most computer installations. Only the subroutines for the determination of the complex function, and Filon's integration scheme would then be needed. Besides the fact that the number of components falls out of the analysis, it is felt that this method has certain other advantages over previous methods. For example, the restriction to equal intervals is merely a requirement of the particular integration scheme used, and

-158is not basic to the procedure as in Cornell's, Householder's, and Prony's methods. Human judgment in the resolution process has been reduced to roughly estimating the spread in half-lives so that a minimum of computational time may be used. This is an advantage over all previous methods. The fact that errors entering the program at different places may be estimated individually may prove advantageous. Finally, the numerical computation itself is not as complicated as in Prony's or Householder's method, and the occurence of two half-lives very close together does not endanger the entire calculation as it does in the former methods.

-159CHAPTER VI NUCLEAR DECAY SCHEME STUDIES In the following sections the results of the investigations carried out on the isotopes Co62, y92 Ir196 and Sc47 are given. Where possible the results are discussed from the point of view of nuclear shell structure. Figure 50 gives a key to the charts of isotopes. A. The Isotope Co62 Relatively few reports have appeared in the literature on the decay of Co62 and the energy levels in Ni62. Pariley, et al.(ll4)in 1949 found that the half-life of the Co62 ground state decay was 13.9 minutes. They measured a maximum n-energy of 2.3 Mev by absorption techniques, and noted the presence of 7 rays of about 1.3 Mev. A much shorter lived activity of about 1.6 minutes that emitted both { rays and 7 rays also was found, and assigned to an isomeric state of Co62. In 1954, Nussbaum, et al.107)measured a maximum 3-ray energy of 2.8 + 0.2 Mev for this isotope in addition to a number of y rays (see Table XIV), and in 1955, Kraushaar, et al.(79)found in the decay of Cu62 a group of 7 rays also ascribed to Co62 (Table XIV). Since over 98% of the Cu62 decay is bo the ground 62 state of Ni62 the results obtained for the 7 rays are somewhat tentative. A search of the literature up to March, 1957, revealed no further work on this isotope. A chart of the isotopes of cobalt, nickel, and copper is given in Figure 51, to indicate the possible isotopes that might be produced by low energy deuteron bombardment of nickel. 1. Chemical Separation The Co62 was produced by the (d,a) reaction on a nickel oxide

-160Cu 60 Element, Mass Number 25 m Half Life I/+.r Mode of Decay Radioactive |Ni60 | Element, Mass Number 26.2 Percent Abundance Stable ft negative beta particle ^+ ~positive beta particle Y gamma-ray K electron capture IT isomeric transition s second m minute h hour d day y year Figure 50. Key to the charts of the isotopes

Cu57 Cu58 Cu59 Cu60 Cu61 Cu62 Cu63 Cu64 Cu65 Cu66 Cu67 3s 81s 25m 3.3h 9.7 rn 69.0.12.8h 31.0 5.2 m 59h _______ r __ /3 y f\YKy 12y 8Kr __ e6,r Ni56 Ni57 Ni58 Ni 59 Ni 6Q Ni 6 N6 2 N63 N64 Ni65 Ni66 6.4d 36 h 67.9 8x104y 26.2 1.2 3.7 85y 1.0 2.56h 56 h K27 f K, y K _ _ f 3, y 1 3 Co55 C056 0057 C058 C059 Co60 Co61 Co62 Co63 Co64 Co 18h 77d 270d 9h 72d 100 11 15.2y 1.7h 14m 1.6m 5m RNK,y /3:K y /39., r iT 3,+Ky 1 > _______________ Figure 51. Chart of isotopes of cobalt, nickel and copper

-162TABLE XIV. RADIATIONS EMITTED DURING THE DECAY OF Co AND Ni62 7 Rays Reference Energy (Mev) Relative Abundance 114 1.3 7/3-1 107 1.0 40 1.17 + 0.03 100 1.5 5 1.7 10 2.0 15 2.5 -3 79 0.66 + 0.03'1 0.86: 0.03'1 1.18 + 0.02 "1 1.36 + 0.02 1.46 + 0.03 1.55 + 0.02 1.67 + 0.02 1.98 + 0.03 2.24 + 0.03 1 Present Worlk 1.17 + 0.01 100 1.17 + 0.01 82 1.47 + 0.02 18 1.74 + 0.03 18 2.03 + 0.03 7 2.5 + 0.2 <2 3 Rays 114 2.3 + 0.1 107 2.8 + 0.2 Present Work* 0.88 + 0.04 25 + 3 2.88 + 0.03 75 + 3 Coincidence Results 7 Rays "in coincidence" Radiation (Mev) (Mev) Present Work 1", 1.2 1.17 7, 1.17 1.17, 1.47, and higher * Stated errors are standard deviations, estimated for the B results.

-163target enriched in Ni. Table XV lists the isotopic percentages, along with the spectrographic analysis for the target material as given by the Isotopes Division of the Atomic Energy Commission. TABLE XV. ANALYSIS OF ENRICHED NICIKEL ISOTOPES Mass Analysis Isotope Atomic % Precision 58 1.99 + 0.10 60 1.20 + 0.05 61 0.14 + 0.01 62 0.77 + 0.03 64 95.90 + 0.16 Spectrographic Analysis Element. Cu 0.05 Fe 0.02 Mn 0.01 Eight bombardments were required to complete this work. Usually about 50 mgs. of the oxide was enclosed in a 1.5-mil aluminum envelope for the bombardment in the 7.8-Mev deuteron beam of the University of Michigan cyclotron. A chemical separation was required to isolate the cobalt in a carrier-free form from the copper and nickel activities produced by the (d,n) and(d,p) reactions in the target element, as well as from the products of these reactions on the impurity elements. Two approaches to this separation problem were tried, both using ion-exchange as the basis for the separation (77). With Dowex-2 resin, nickel will not adsorb onto the column at any concentration of HC1. Cobalt and copper, however, both show appreciable adsorption in

-164the range from 6 to 12 N HC1, with cobalt being the more strongly adsorbed element. In 6 N HC1, cobalt and copper are adsorbed about equally. In the range from 6 to 3 N HC1, cobalt rapidly loses its tendency to be adsorbed while copper remains fairly constant. Using two different sized resin beds, one 8 nmml x 140 Immu and the other 4 I mm x 10 imm, it was found that trying to remove the copper before the cobalt by eluting with solutions in the range of 8 to 10 N in HC1 required mlore time than removing these elements in the reverse order. While the smaller column decreased the tim;e of the separation, the degree of separation was much poorer at elution rates of 1-2 drop/sec. It was therefore decided to perform the separation as indicated in Table XVI using the larger column. Since it was necessary to use the target material in several bombardmlents and to return it to the Atomic Energy Commission in purified form, a procedure for reprocessing the target material was required. This was done as follows: The 8 N HC1 solution containing Ni(II) was first evaporated to dryness, and then the residue was dissolved in water. Next Ni(OH)2 was precipitated with NaOH. The hydroxide was repeatedly washed in a centrifuge tube with distilled water until the wash water was neutral. Finally, the hydroxide was heated until the oxide was obtained. The cobalt fraction was mounted directly by evaporation onto a thin cover glass for 7-ray measurements. For the 3-ray measurements the cobalt fraction was first evaporated to dryness, then taken up in water, and finally mounted on the aluminum coated 1/4-mil Mylar film previously described.

-165TABLE XVI. COBALT SEPARATION CHEMISTRY Element separated: Cobalt Procedure by: Gardner Target Material: Nickel Oxide Time for sep'n: 1/2 hr. Type of bbdt: 7.8 Mev deuterons Equipment required: Ion exchange column of AG 2-x8, 200-400 mesh resin obtained from Bio-Rad Labs. Resin bed was 8 iim x 140 im. Degree of purification: Estimated, 105 from Ni and Cu. Advantages: Carrier-free separation with high decontamination. Procedure: 1. Dissolve NiO in 5 ml of 10 N. HCI. Evaporate to small volume, place in column previously washed with conc. HC1. 2. Remove Ni from column by elution with 15 ml of 8 N HC1 at rate of 1 or 2 drops/sec. Both Co and Cu remain on column. 3. Remove Co by elution with 6 - 8 ml. of 4 N HCl. Last ml. or so possibly contaminated with small amount of Cu. 4. Mount for counting. Remarks: 1. General reference: Hicks, H. G., et al., "The Qualitative Anionic Behavior of Resin,'Dowex 2' ", Livermore Research Laboratory Report, LRL-65, Dec. 1953.

-1662. Experimental Me thods y-ray spectra were obtained using both the photographic technique, and also by automatic plotting of the spectrum from the pulse height analyzer on the strip chart recorder. The P-ray spectrum was measured using the 3.6-Mev hollow plastic scintillator. The 7-7 coincidence studies were carried out using two 1" x 1-1/2" NaI (Tl) crystals. A 1.6 gm/cm2 aluminum absorber was placed between the sample and each crystal to absorb the 0-rays. For the 1-7 coincidence studies, one NaI (T1) crystal was replaced by the plastic scintillator previously mentioned. No absorber was used between the source and the plastic scintillator. Since the half-life of Co62 is so short, it was only possible to obtain coincidence information through the use of the photographic coincidence method described in Chapter III. 3. Results The half-life of Co2 was determined to be 13.909 + 0.013 minutes, the error being the standard deviation, by following the decay of the 7 rays in a scintillation well-type counter. The decay curve is shown in Figure 52. A small amount of Ni65 is seen to be present. No other cobalt activity could be seen. The resolution was obtained using the Weighted Least-Squares program. Figure 53 shows the Co 7-ray spectrum taken from the strip chart recorder. Measurements of this type along with the oscilliscope pictures indicated that the primary 7-ray transition was 1.17 + 0.01 Mev in energy. Four higher 7 rays of lower intensity were also observed. These were 1.47 + 0.02, 1.74 + 0.03, 2.03 + 0.03, and 2.5 + 0.2 Mev.

-167DECAY CURVE OF Co FRACTION 104 H DD e 0 HL0 2: c-J o100 Co62 0 1 2 5 4 5 HOURS Figure 52. Decay curve of cobalt fraction

COUNTS PER MINUTE (,O ca ~~~~O 3m H C'] 0~ 0 )LD'~~C) ror',,o0M \ C:~ rr^^ 4 c+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~', C")^V ^ ^^ ^/ ^ ^^ ^ 0~~~, _r>^ ^\~~~~~~~~ ^o) ^ c^_C __. - ^^^ ^^ ^~~~~~~~~~(,

-169Integration under the Speedomax curves, correcting for the NaI(T1) crystal efficiency (92)yielded intensities of 11, 11%, 4%, and less than 1%, respectively, for these 7 rays relative to 100% for the 1.17-Mev 7 ray. A tentative decay scheme based on the 7-ray information allowed for the possibility of three,-ray groups: a 0.8, a 2.3 to 2.5, and a 2.8-Mev transition. Due to the short half-life, a complete:-ray spectrum could not be obtained on a single sample with enough points or statistical accuracy. Therefore, it was necessary to combine the results of two bombardmnents. Figure 54 shows the results of this combination. Each point on the curve was corrected for decay during the time of the count which ranged between 3 to 8 minutes. Then each point was corrected for decay back to the time of the start of the first count. Correction was also made for the 7-ray background in the plastic scintillator. Since the upper portion of the P-ray spectrum yielded a Kurie plot straight down to 1.5 Mev, the possible 2.3-2.5 Mev transition was assumed to be absent, and the two portions of the Kurie plot obtained on different days could be fitted together easily. The resulting total Kurie plot, when corrected for instrumlent resolution, could be resolved into two components: a 0.88 + 0.04 and a 2.88 + 0.03 Mev transition. Relative abundances are 25 + 3% and 75 + 3%, respectively. Errors in the energies and intensities are estimated probable errors. Log (ft) values are 4.5 and 5.8, respectively, indicating both transitions are allowed. This is borne out, in the case of the 2.88-Mev n-ray, by the shape of the Kurie plot. While the shape of the lower energy group also appears to be allowed, it is subject to the usual subtraction errors.

o KURIE PLOT OF o Co62 0 O-SECOND BOMBARDMENT 0 X-THIRD BOMBARDMENT @-LOWER ENERGY COMPOr- ^0 NENT (RESOLVED) F 0 0 0.4 0.8 1.2 1.6 2.0 2.4 2.8 ENERGY (MEV) Figure 54. Kurie plot of Co62

-171Assuming the 2.88-Mev f group populated the 1.17-Mev 7 level, the 7-ray intensities calculated from the observed B-ray intensities differed from the observed 7-ray intensities by a factor of 2. Furthermore, the total disintegration energy for Co62 predicted from ( systematics is about 5 Mev, whereas the 2.88-Mev 1 and the 1.17-Mev 7 transitions add up to only about 4 Mev. By assuming a cascade doublet following the higher B ray, the difficulties were resolved. Both 7 rays would have an energy of 1.17 Mev, and differ in energy by less than 1-2%. The corrected percentage abundances for the 1.17 (doublet), 1.47, 1.74, 2.03, and 2.5 Mev 7-rays are then, respectively, 100% and 82% for the doublet, 18%, 18%, 7% and less than 2%. Coincidence information was obtained by photographing the y rays in coincidence with all rays above 1.2 Mev. This showed only the 1.17 Mev 7-ray. 7-7 coincidence information was obtained by setting the interdicting channel on the 1.17-Mev 7-ray. The results show strong coincidences with a 7 ray at 1.17 Mev, and also with higher gammas. 4. Conclusions The above information has been correlated into the level diagram shown in Figure 55. The P rays are shown with their energies in Mev's first, followed by the % abundance and the log (ft) value. The 7 rays are shown similarly with their energies in Mev's first and then the % abundance. Tn the nuclide Co6 the ground state configuration of the protons is (if 7/2)7/2, where the superscript -1 indicates one nucleon less than a closed shell and the subscript 7/2 indicates the resultant proton spin. The ground state configuration of the neutrons is somewhat

-1721.6m (5+) 14.0m Co62 0.85 (25%o,4.5) (4+) 1.73(18%) 2.88 (75%,5.8) \ 2.03(7%) (3+)______ 1.47(18%) 1.17(8 2%) I(2) ______ _.64 (< 2 %) (2+). 17( 100%) (0+4 Stable N i 62 Figure 55. Decay scheme of Co62

-173less uniquely predictable. Here both the (2p 3/2) and the (if 5/2) levels are relatively close together in energy (91). A consideration of the other nuclides in this neutron range indicates that the probable configuration is (2p 3/2)372 (if 5/2)o-2 According to Nordheim's "weak" rule (103) the proton and neutron spins will tend to add, giving a resultant spin of 4 or 5 and even parity. 62 In the case of 28Ni34, the (if 7/2) proton shell is clearly filled with a resultant spin of 0. Again the neutron configuration is unclear, but is probably (2p 3/2)02 (if 5/2)02 so that the ground state 62 of Ni would have 0 spin and even parity. The first excited state of 62 Ni is very probably 2+, that is spin 2 and even parity. The second state at 2.34 Mev and the fourth state at 4.37 Mev are both required to have high spins and even parity since the ground state of the parent is 4+ or 5+ and both transitions are allowed. A spin of 4+ for the second state is reasonable from the shell model. The 2.5 + 0.2 Mev 7-ray transition cannot be a cross-over from the second excited state to the ground state, and therefore, must arise from the third excited level. Hence the order of the 1.74-and 1.47-Mev y rays given in Figure 55 is the most likely. This would require a spin of 2 or 3 for the third level. The fact that the third level is not populated by 3 decay, together with the ratio of the intensities of the 1.75 and 2.03-Mev 7 rays suggest the following level assignments: spin 3+ for the third level, spin 4+ for the fifth level, and spin 5+ for the ground state of Co62. The 1.6-min. isomer of Co62 found by Parrlley (108) was not found here due to the length of time required by the chemical separation.

-174However, it is quite reasonable to expect that the first excited state of Co62 would have a proton configuration of (If 7/2)/2, and a neutron configuration of (lf 5/2) 3. These would couple to give a spin of 1+, pro5/2 ducing an isomeric state with an energy not much above the ground state. This first excited state would have a high probability of decaying directly to the ground state of Ni62. For the sake of completeness this hypothetical excited level has been included in Figure 55. B. The Isotope y92 Several representations of the decay scheme of 92 have been reported in the literature, the most recent being found in the compilation of K. Way, et al. (144). This last was based primarily on the work of Cassatt and Meinke (20,21), and of Ames and his co-workers (4). Furthermore the investigation of Nb92 which decays by electron capture 92 (54,66) has contributed to the knowledge of the energy levels in Zr2. The half-life of 92 has been determined to be approximately 3.6 hours (4,20,21) Maximum P-ray energies of 3.4, 3.5, and 3.6 Mev have been reported (4,20,21,138). The presence of three P rays of maximum energies of 1.3, 2.7, and 3.60 Mev was postulated by Ames, et al. (4). y rays have also been measured (4,20,21,54,66). Table XVII lists the available information concerning the 3 rays emitted by y92 and the 7 rays from Zr92 Figure 56 shows the decay scheme for y92 - Zr92 _ Nb92 as given in reference (144). The present work was undertaken primarily because the placement of the 0.45-and the 0.21-44ev rays is not consistent with the available data, and the assignment of a spin of 2+ or 3+ to the ground state of Nb92 is unlikely from shell structure theory. A ground state spin of 2+ or 3+ would necessitate a spin of 6 or 7 for the

-175TABLE XVII. RADIATIONS EMITTED DURING THE DECAY OF y92 AND ZR92 y Rays Reference Energy (Mev) Relative Abundance 20, 21 0.20 10 0.48 (double?) 11 0.94 (double?) 19 1.45 10 1.9 1 2.4 0.2 4 0.45 0.55 0.93 1.39 1.83 54,66 0.900 0.930 1.83 2.35 Present Work * 40.03 (postulated) 6 - 12 0.07 + 0.004 1.9 + 0.3 0.14 - (?) weak 0.21 + 0.005 3.1 + 0.7 0.29 + 0.02 3.1 + 0.7 0.47 + 0.02 (double) 7.2 + 0.5 o.54 + 0.02 0.5 + 0.2 0.90 + 0.01 3.9 + 0.5 0.93 + 0.02 11.2 + 1.5 1.44 + 0.04 4.2 + 0.5 1.86 + 0.05 0.5 + 0.2 2.40 t 0.05 0.08 + 0.05 B Rays 20, 21 3.6 (spectrum complex) 4 1.3 11 2.68 12 3.60 77 Present Work * 1.26 + 0.06 9 1.75 + 0.15 (complex?) 3 3.60 + 0.03 88 * Errors are standard deviations, estimated for the B results.

18 3S Kr...... 16 ~~~~~~~ Kr92 ~ ~ -—'~ 36 56 Rb91+p 14 tu 14 -12.48 7.35 Mo91+n 12- ~~~~Sr9i+n AN- 8S Short 92 Sr91+P 9.. y91'n~~~~~~~b.6 1 ~,-~,;.....X9.453 lived37 ~ 6.61 Y sn A/ 92 Eo (ev).......A5''m'mc 8 ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~~~N I" \9A5 43 49 O+ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ +......~Z 2s N9+ n........'~ 4.1 F 6.305. 2.7h Sr92 I 5- 38 54 (4.) 7.71 C 8.66 1.93'v'O.55 3893!T /00I% - 4..3 1.. 4 2, 22 9 ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~~~~9 36th y0. o 3 39 53, 0.4 _ 13 h 1.3 11%/0,6.7 r2A2 2:35,5 0.55 23+ 2- 2.68 9pj 9 A 12% o 8.0 10 d4Nb 0.5~0.2 O 1.4 1 1 5 1: Q9'P. 5 Stable M050 b (d,p). 3.60''.8 2A 09 0 /o I77.83.20 3.5%6<6.6 21.2 42 50 t (n,2n) U (V,n) <1.33,1 un ~0.93.96.5%e,~ ~;.o!A,B,C Comments b 6%6' 60ML Moss Unk 0- + i II I Stable oZr 92 40 52 Figure 56. Decay Scheme for y92 Zr92 NbW92 as Given by K. Way, et al.

-17792 isomeric state in Nb, whereas a low spin and odd parity would be expected from shell theory. If this isomeric state actually populates the 2.35-Mev level in Zr92, as shown in Figure 56, which in turn may decay directly to the ground state, then the spin of the isomeric Nb92 level could not be more than 4 or 5. Finally if the ground state spin of Nb92 is actually 5+ or 6+, as would be expected, then the 0.93-Mev level shown in Figure 56_with a spin of 2+ would not be fed by a transition from the Nb9 ground state. These and other considerations indicated that a reinvestigation of this decay scheme was warranted. A search of the literature up to March, 1957 revealed-no further work on this isotope. A chart of the nuclides of yttrium, zirconium, and niobium is given in Figure 57, to indicate the possible isotopes that might be produced by low-energy deuteron bombardment of zirconium. 1. Chemical Separation The yttrium activity was produced by bombarding high-purity zirconium metal foil with 7.8 - Mev dieuterons in the University of Michigan cyclotron. The foil was obtained from the Foote Mineral Company, Philadelphia, Pennsylvania, and the chemical impurities are listed in Table XVIII. A chemical separation was necessary to separate the yttrium from the zirconium and niobium activities produced by (d,n) and (d,p) reactions in the target element, as well as from products of these reactions on the impurity elements. Seven cyclotron bombardments were required to complete this study. When the yttrium chemical separation was being planned, a possible study of the 17-minute isotope Y94 was considered. This study would have required a sample of zirconium enriched in the isotope Zr96. Although

Nb90 Nb91 Nb92 Nb93 Nb94 Nb95 Nb96 Nb97 Nb98 24S 15h 64d long 13h 10.d 5.6m5x 90h 35d 23.35h 60s 72.1m 30m(?) IT IT KKly IT IT lT Py y IT /Ir _ - Zr88 Zr89 Zr90 Zr91 Zr92 Zr 93 Zr94 Zr95 Zr 96 Zr97 85d 4.4m1793hl 51.46 11.23 17.11 9.5x105y 17.40 65d 2.80 17.h O Y87 Y88 Y 89 Y 90 Y91 Y92 Y 93 Y 94 Y95 14h 800h 104d ~1100 64h 51Dm 58d 3.60h 10.Oh 16.5m 10.5m IT T'K {4si lr K 7y rIT,- IsTIp f y r,,r Ir Figure 57. Chart of the isotopes of yttrium, zirconium, and niobium

-179it was later decided not to attempt this study due to the cost of the enriched zirconium isotopes, the chemical separation was developed with the short half-life of Y9 in mind. TABIE XVIII. CHEMICAL IMPURITIES IN TYPICAL ZIRCONIUM TARGET FOIL AS DETERMINED BY SPECTROGRAPHIC ANALYSIS Element % Element % Al 0.04 - 0.08 Mo 0.001 Ca 0.003 N 0.004 - 0. Cr 0.001 Pb 0.001 Cu 0.007 Si 0.02 Fe 0.01 - 0.20 Sn 0.001 Hf 2.5 - 3.0 Ti 0.04 Mg 0.003 W 0.001 Mn 0.001

-180As in the case of the cobalt separation previously described, an ion-exchange separation was decided upon (78). It seemed desirable to introduce a rapid, fairly thorough separation of yttrium from zirconium and niobium before the ion-exchange stage. This initial step would then 94 constitute the entire separation if the Y isotope were to be studied. However, this first step would be used in the more thorough separation, and hence should not introduce any inactive yttrium carrier. The zirconium metal foil used in the Y2 bombardments could best be dissolved in dilute HF. Niobium and zirconium are soluble in HFj while yttrium would form the more insoluble YF3. Therefore, if an element which forms an insoluble fluoride were to be added to the HF solution of the zirconium target foil, the radioactive yttrium would be carried down on the insoluble fluoride precipitate. The alkaline and rare earth elements would be suitable, but they would present difficulties in designing a rapid separation to remove them at a later stage. For these reasons it was decided to use lead as the carrier for yttrium. Other advantages of lead are that its fluoride can easily be dissolved in HC1 solution, and that it can be rapidly separated from yttrium on an ion exchange column. One disadvantage is that the fluoride tends to dissolve slightly when washed with water. The principal place in the separation where some yttrium is most likely to be lost occurs in the first step where the target foil is dissolved in dilute HP. It was noted that an appreciable amount of the yttrium activity tends to adhere to the walls of the Lusteroid centrifuge tube. When the PbF2 is precipitated, and in later washings, the precipitate should be thoroughly stirred in order to remove as much yttrium from the walls of the tube as possible. The actual chemical separation that was developed is given in Table XIX.

-18iTABIE XIX. YTTRIUM SEPARATION CHEMISTRY Element separated: Yttrium Procedure by: Gardner Target Material: Zirconium Time for sep'n:;1 hr. Type of bbdt: -7.8 Meev Equipment required: Lusteroid centrifuge tubes, two ion exchange columns of deuterons AG 2-X8, 200-400 mesh resin obtained from Bio-Rad. Lab. Resin beds were 8 mm x 140 mm. Yield: ^90% Degree of purification: Estimated >106 from Nb and Zr Advantages: Essentially carrier-free separation with high decontamination. Can be shortened to' 5 min. if carrier can be tolerated. Procedure: (1) Place target foil in a Lusteroid centrifuge tube containing 5 ml of 6 N HF, and 10 mg of Nb holdback carrier. (Sample dissolved immediately. ) ++ (2) Add 10 mgs of Pb, stir precipitate, and centrifuge. (See remark 1). (3) Wash precipitate twice with'"5 ml of ~1 N HF, centrifuging in between. Wash once with a minimum amount of H20. (4) Transfer precipitate to a glass beaker. Dissolve in several ml. of conc. HC1 and evaporate to dryness with a few drops of HNO3. Take up in conc. HC1 and evaporate to dryness again. (5) Take up in several drops of conc. HC1 and place on column previously washed with conc. HC1. (6) Elute with conc. HC1. at rate of'1-2 drops/sec. The Y will begin to come off after 3-4 mls of eluent; most will appear in another 4-6mls. (See remarks 2 and 3). (7) If desired, the Y may be separated from the lead as follows. Concentrate the eluent from step 6 to <1 ml, dilute with H20 untilv%5-7 N in HC1, and place on column prepared by washing with 6 N IHC. Elute with 6 N HC1. (See remark 4).

-182TABIE XIX. YTTRIUM SEPARATION CHEMISTRY (Cont d) (8) If the Pb carrier can be tolerated, along with an almost negligible amount of Nb and Zr, the precedure may be shortened to ~ 5 minutes if the sample is mounted after step 3. Filtering, using non-glass equipment, will reduce the time even further. Remarks: (1) Essentially all of the Y follows the Pb. (2) A large part of the Pb will also appear in Y fraction. (3) Any small amounts of Zr and Nb possibly present will remain on column. (4) Pb will remain on column. (5) General reference for this type procedure: Hicks, H. G., et al., "The Qualitative Anionic Behavior of a Number of Metals with an Ion Exchange Resin,'Dowex 2' ", Livermore Research Laboratory Report, LRL-65, December, 1953.

-1832. Experimental Methods 92 The 7 rays from 92 were studied in the usual manner. The 3 rays were measured using the 3.6-Mev hollow plastic scintillator. In obtaining the P-ray spectrum of y92, corrections were made for the contribution of Y and Y which were produced along with the Y92 The y - y coincidence studies were carried out using two 1" x 1-1/2" NaI(Ti) crystals, with 1.5 gm/cm2 Be absorbers between the source and the detectors. For 3 -7 coincidence studies one 7-ray detector was replaced by the 3.6-Mev P scintillator. Since the half-life of y92 is relatively long, manual coincidence spectra could be obtained as well as photographic coincidence results. 3. Results 92 The half-life of Y was found to be 3.66 + 0.66 hours by a least-squares analysis of the decay curve of a portion of the P-ray spectrum at 3 3.0 Mev. This is beyond the endpoint of the Y9 P-ray spectrum. Figure 58 shows the decay curve and the least-squares line fit. Other decay curves taken with a well-type y-scintillation counter showed the presence of only three components, a 3.6-hour, a 65-hour, and a 105-day activity 92 90 88 due to y, y and y88 respectively. The 7 rays found, and their relative intensities are listed in Table XVII. Figure 59 shows the 7-ray spectrum taken from the Speedomax strip chart recorder, while Figure 60A shows an oscilliscope picture of the same spectrum. These results reveal the presence of four intense y rays with energies of 0.21, 0.47, 0.93, and 1.45 Mev along with several weaker transitions. The 7 intensities were obtained by integration under the Speedomax curves, correcting for the NaI(Tl) crystal efficiency (92).

6.6 DECAY CURVE FOR y92 6.4 D62 0 HALF- LIFE 3.66 HOURS > 0 6.0 50 100 150 200 250 MINUTES Figure 58. Decay curve of y92 with least-squares line fit

-185[ I l lIl.I 0.07 0.21 92 10:5 103 0 0.47. 0.29 H- r( D L I l1l, 0.93 LU 054 - cn 0.90 0.2 0.4 0.6 0.8 1:0 ENERGY (MEV) Figure 59. Gamma-ray spectrum of Y92

(A) (B) (C) Figure 60. (A) Total y G ania-ray Spectrum. (B) Gammaray Spectrum in Coincidence with 0.47-Mev Gamma ray. (C) Gamma-ray Spectrum in Coincidence with 0.21-Mev Gamma ray.

-187Figure 61 gives the P-ray spectrum for y92, and has been corrected for decay as well as for the Y90 and the y88 components. The conversion peak at ^ 0.45 Mev should be noted. After a decay scheme had been postulated the total conversion coefficient could be estimated, and was found to be aT 0.20 + 0.12. Using the tables of Rose, et al. (119) it was found that the measured conversion coefficient fell into the range of the following transitions: M3, M4, and E4. This indicated that the 0.47-Mev 7-ray arose from a spin change of 3 or 4 units, and originated at a meta stable state whose half-life was on the order ofc 2 milliseconds to 5 minutes depending on the transition involved. Figure 62 illustrates the Kurie plot of the y92 spectrum. Here the major transition has a maximum energy of 3.60 + 0.03-Mev and a relative abundance of~ 88%. A small amount, roughly 3%, of an intermediate transition with a maximum energy 1.75 + 0.15-Mev appeared to be present, along with a o 9% of a 1.26 + 0.06-Mev transition. The 3.60-Mev transition definitely had a "unique" first forbidden shape, but the shapes of the lower energy transitions could not be determined with accuracy. No trace of the 2.68-Mev transition found by Ames, et al. could be seen; an upper limit for this transition, consistant with the present data, would be 1 or 2%. o-yr, -3 and 7-y coincidence data were obtained both by the photographic method and also by point-by-point manual sweeps. These are included in Table XX. Figure 60B is an oscilliscope photograph of the pulses in coincidence with the 0.47-Mev 7-ray. One of the more interesting features of the coincidence data is that the 0.93-Mev 7-ray does not seem to be in strong coincidence with either 1 transitions or other y transitions.

10,000 BETA SPECTRUM OF Y92 8,000 LLJ D: 6,000 H 4,000 0 () 2,000 0.5 1.0 1.5 2.0 2.5 3.0 ENERGY (MEV) Figure 61. Beta-ray spectrum of y92

o KURIE PLOT OF Y92 0 0 0 0 ^^^ o o F 000 0.5 1.0 1.5 2.0 2.5 3.0 5.5 ENERGY (MEV) Figure 62. Kurie plot of Y92. Ordinate scale for top curve differs from scale for lower curves

-190TABLE XX. COINCIDENCE RESULTS FOR ISOTOPE y92 Type of Energy of radiation Energy of radiation Energy of radiation Coincidence (Mev) in "strong" coinci- in "weak" or undence (Mev) certain coincidence (Mev) - 7 0.21 0.29, 0.47 0.14, 0.54, 0.67, 0.93 0.47 ~0.21 0.07, 0.93, 0.55 0.14, 0.4-0.9 0.93 0.07, 0.47 0.54, 0.9, 1.44 -7 ~o0.9 -y0.2 7 rays below 0.21 and at least up to 0.93 1.6 ~0.2 0.9 3.0 none none 7 —03 0.45 1.2 0.93 1.6 This raises the possibility that the 0.93-Mev level is fed primarily by a transition from astate that has a lifetime in the order of 1 2 second, since the resolving time of the coincidence analyzer was 2 pA seconds.

-1914. Conclusions The above information has been correlated into the level diagram shown in Figure 63. In the nuclide Y9 the ground state configuration 39 53 of the protons is (2p )1/2 where the subscript indicates the resultant proton spin and the supercript indicates the number of nucleons in the shell. The neutron configuration is (ld5/2)5/2. The resultant spin is 2-, that is, spin 2 and odd parity. This assignment is verified by the "unique" forbidden shape of the 3.60-Mev 3 transition. The ground state configuration 92 1 2 2 of 40Zr52 would then be (2 ) (ld5/2) for the protons and neutrons respectively. The resultant spin is then 0 +. Now the ground state configuration of 41Nb51 should be (lg9/2)9/2 (ld5/2)/2. Here the proton and neutron spins will tend to add to give a possible spin of 5+, 6+, or 7+. 92 A possible configuration for the first excited state of Nb would be 2p )1/2 (lg9/2)] [(ld5/2)5/ yielding a spin of 2-, although a P-2)1/2 0 5/2 lower spin might also be produced depending upon how the (lg9/2) protons are coupled. Assuming a ground state spin of 5+ or 6+ and a first excited 92 state of 1- or 2-, the isomerism of Nb is understandable. However, since the first excited state of Zr9 almost necessarily has a spin of 2+, the 92 electron capture from the ground state of Nb cannot proceed to this state. 92 Two possibilities then arise. Either the first excited state in Zr is very low in energy (on the order of a few tens of kilovolts), or else it is quite close in energy to the 0.93-Mev level. The first possibility is unlikely from the systematics of the first excited states of even-even nuclei, and also because the 3.60-Mev f transition from Y9 clearly has the unique shape indicating a spin change of two units and a change in

(2-)3.66hr Y92 (,2-) 13hr Nb92 MEV 1.20 (9%,6.7), 2.40 (56 I)Od ^^ \ \\ I I ^~ {~~~~~~5f6+) 10d 0.54 0.47 I 1.67-1.74(3%77.3) +); _ _1.93 +, I 1o)o 7 1: 1 0.21 1.86 ^^ ~ ~ __ \H^^ ^7^ 1___ 0 ^ (516-) I______ ___ o0 14 _1^2 1.72 0. 90 1.44 0.29 <0.41.S,6.)n 3.60(88%, 7.9) (7,84) 0 1.2 1.43<0.4(3.50/o,<6.6) 1.86 0.47 \ (4 I) _ t 0.9 <1.3 (96.5 %-6.0) \ 2 _+) _ _0.03 _ _ 0.93 F>2x10^ S I HALF-I LIVES 2.40 I1x1 I \,S 0.93 ((0+) V f 0 STABLE Zr92 Figure 63. Decay scheme of Y92 - Zr92 - Nb92

-193parity. Further, Griffith (48) has found the first excited state to be 0.926 Mev using inelastically scattered 4-Mev neutrons. Therefore, the 92 second excited state ii Zr was postulated to lie just above the 0.93Mev level. Accurate 7 ray spectra could not be obtained in the kilovolt region, but a reasonable order of magnitude for the difference between the first and second excited states might be AQ30 Kev or less. If the second excited Zr9 state has a spin of 4+, which could then be populated by electror capture from the ground state of Nb92, the ^ 30 kev transition might have a life time of the order of a few p seconds as indicated by the coincidence data. While this transition should be converted to a large extent, the conversion electron peak would be too low in energy to be seen on the n-ray spectrometer that was used. The 2.38-Mev level, being populated by transitions from both Y92 and the isomeric Nb92level would then have a spin of 2- or 3-. This is substantiated by the presence of a 2.40-Mev 7 ray. Since the I.86-Mev level can decay directly to the ground state, a spin of 3+ is reasonable. A negative parity is unlikely because the Y 3 transitions to this level are in such low abundance. Therefore, the weaker electron capture from the Nb92 ground state cannot populate this level, but more likely reaches the 1.93-Mev level for which a spin of 4+ is likely. The 1.72-Mev level is probably not populated by either electron capture or negative { emission, indicating a spin of 5- or 6-. A high spin at this level is likely because it must lead to the 1.43-Mev level having a still higher spin. The 1.43Mev level was postulated to be a meta-stable state leading to the converted 0.47-Mev 7 cransition. Therefore, a spin of 7+ or 8+ is required. The

-194isomeric level seems best placed at this point because of the coincidences found between the 0.21-Mev and the 0.29-Mev 7 transitions. It was not possible to determine which level the 1.44-Mev 7 ray populated. Either the first or second excited states or both are possible. The decay scheme presented in Figure 63 appears to fit all of the coincidence data together with the data on 3- and 7-ray intensities. A more concrete assignment of spins must await angular correlation and/or better conversion electron data. 196 C. The Isotope Ir9 196 It had been noted in 1953 that an isotope of iridium, Ir96, was listed in isotope tables and compilations of nuclear data (60,101) as having a half-life of 9 days and was supposed to decay by the emission of 196 an ~ 0.08-Mev 3 ray. A half-life of that order was surprising for Ir since the isotope Ir 94 has a half-life of only 19 hours. Furthermore, the half-life and n-ray energy indicated that the transition was allowed, an unlikely occurrence for this section of the isotopic table. In January, 1954 Butement and Poe (19), who had contributed the original information on this isotope, published further information substantiating their previous findings and listed 7 rays at energies of 0.58, 0.76, ando 1.0 Mev. A search of the literature was made in 1954 and has been extended periodically through March, 1957. No further work on this isotope has been found to date. A chart of the isotopes of iridium, platinum, and gold is given in Figure 64 1. Chemical Separation The radioactive iridium was to be produced by the (d,a) reaction from deuteron bombardment of metallic platinum foils 1.05 mils in thickness. Therefore, a chemical separation was needed to obtain the iridium free from

Au'90 Au191 Au192 Au'93 Au194 Au195 Au196 Au197 Au198 Au199 Au200 3 h 4.8d 4s 17h 39h 30s80d14h 5.6 d 2.7d 3.2d 48m _____K, y,K,y IT K,y K IT KT IK y ~, /~y pt 89 Pt190 Pt 191 Pt 192 Pt193 Pt194 Pt95 p96 Pt 197 Pt98 Pt199 0Od3 3d 3 d 1.4h 19h 0.012 0.78 jiong 32.8 33.7 25.4 7.2 30 K,_ K, IT iT _', Ir 188 Ir 189 Ir 190 r 191' Ir 192 r93 r 194 Ir Ir 196 Ir 198 40h 10d 3h ^10d 1.4m 74d 19 h 2.3 h short 7m 50s ^K./ K ^ K l 38.5 161.5 3KFy K. CKt, pKy a, l Figure 64. Chart of the isotopes of iridium, platinum and gold

-196contamination by platinum and gold activities produced by (d, p) and (d, n) reactions on the platinum target, and also from activities produced by reactions on the impurity elements in the target. Three bombardments were obtained on "commercial grade" platinum foil, 99.5% pure, obtained from Baker and Co., Inc., Newark, New Jersey. Table XXI shows the spectrographic analysis for a typical foil. No numerical estimates of the impurities were available from Baker and Co. TABIE XXI. SPECTROGRAPHIC ANAIYSIS OF TYPICAL PLATINUM FOIL Element Amount Code Pt M M - major Ir VST + Os VST S - strong Pd ST Rh ST VST - very strong trace Ru T Au T ST - strong trace Ag T Cu T T- trace Fe ST - Ni T 0 - not detected Pb T Sn T Zn T Si T Mg T Ca ST Mn T Al ST A final bombardment was obtained on a sample of platinum enriched in the 198 isotope Pt

-197The chemical separation that was finally developed is based to a certain extent on the information given in Noyes and Bray (104), and also in Meinke's compilation (93,94). Table XXII shows the separation. TABLE XXII. IRIDIUM SEPARATION CHEMISTRY Element separated: Iridium Procedure by: Gardner Target Material: Platinum Time for sep'n: 8 hours metal Type of bbdt: 20-Mev Equipment required: No deuterons special equipment Yield: 70% 6 Degree of purification: - 1 Advantages: Separates Ir from Pt, Au, Cu, Ni, and Zn with high decontamination factor. Procedure: (1) The Pt target was dissolved in boiling aqua regia and 3 mg Ir carrier and 10 mg each of the following carriers added; Cu, Ni, Zn and Au. (2) The solution was evaporated to incipient dryness, diluted to 10 ml with water and 10 drops conc. HC1 added. (3) Au was extracted 6-8 times with ethyl acetate. (See remark 1) (4) The aqueous phase (yellow-brown) was treated with 1-2 drops H2NNH2 to destroy NO0. (Solution turns pale). )See remark 2) (5) Pt was reduced with SnC1 (solution turns deep red) and extracted 6-8 times with ethyl acetate. (See remark 3) (6) The aqueous phase was evaporated to dryness with 2-3 ml aqua regia and 5 mg Pt carrier added. The NO was removed by addition of 20 drops conc. HC1 and evaporation to dryness. (7) The residue was dissolved in 2 ml H20 and 4 drops 6 N HC1.

-198TABIE XXII. IRIDIUM SEPARATION CHEMISTRY (Cont'd) (8) The solution was saturated with solid NH4Cl, warmed to dissolve any excess and cooled in ice for 0.5 hr. The red precipitate of Pt and Ir was washed several times with saturated NH4Cl solution. (9) The precipitate was dissolved in hot water, and NH4, removed by evaporation to dryness with 2-3 ml aqua rogia. (10) The residue was dissolved in 2-3 ml H20plus 2-5 drops conc. HC1 and evaporated to dryness again. (11) The residue was dissolved in 6-8 ml H20 and made basic to litmus with a few drops of saturated Na2CO3 solution (solution turns from brown to yellow). (12) The solution was heated to boiling, 4-6 ml NaOBr solution (0.5 ml 1M Na2CO 1 ml saturated Br2 solution) added and heating continued until the solution turned greenish-blue. (See remark 4) (13) 1-2 drops 6 N HC1 was added to the still warm solution and the solution digested until the IrO2 coagulated. (14) The precipitate was washed several times with H20, dissolved in conc. HBr and evaporated to dryness. The residue was dissolved in H20 and mounted for counting. Remarks: (1) For both the Au and Pt extractions, aqueous phase 2-4 organic phase (2) Metallic Pt may precipitate from hot solution. (3) If a red precipitate forms, conc. HC1 is added until it dissolves. In the presence of a large amount of Pt a series of partial reduction and extractions to remove all the Pt is preferred. (4) The addition of a drop of 1 M Na CO solution and/or saturated Br2 solution may expedite the formation of the blue-green color.

-199It was later found in the bombardment work that it was necessary to perform the chemical separation only once. If the final iridium activity was reprocessed using the same separation, no noticeable improvement in the purity could be detected. In the final enriched isotopes bombardment it was necessary to return the target material in a purified form. This was done as follows: The ethyl acetate containing the platinum was evaporated to dryness, and then strongly heated to drive off organic materials. In the process the Pt(II) was reduced to metallic platinum. The residue was washed repeatedly with hot concentrated and dilute HC1. The metallic platinum was dissolved in aqua regia, evaporated to dryness,and again reduced to metallic platinum by heating. After repeated washings, the metallic platinum was dissolved in aqua regia, and PtC14 formed by evaporation to dryness with HC1. The PtC14 was dissolved in water slightly acid with HC1. Next, metallic platinum was precipitated from solution by the addition of metallic zinc. The zinc was removed by dissolving in dilute HC1. The remaining platinum was washed in HC1 and water, and finally dried. 2. Experimental Procedure The three bombardments on natural platinum foil were obtained using the Argonne National Laboratory cyclotron. The cyclotron at the University of Michigan could not be used since a beam energy higher than 7.8-Mev is required to obtain useable (d, a) reaction yields at atomic numbers in the rare earth region and above. The Argonne cyclotron supplies deuterons at an energy of 20.4 Mev. Using stacks of 2 and 3 platinum foils separated by copper absorbers of appropriate thicknesses (calculated using the curves of Aron et. al., reference 5), a total of 7 platinum targets were bombarded at various energies ranging from 9.6 to 20.4 Mev.

-200Figure 65 shows the foil arrangements for each of the three bombardments. All platinum foils are 1.05 mils thick, while the thickness of the copper absorbers are indicated. The energy of the deuteron beam striking each platinum foil is also shown. The range and rate of energy loss of deuterons in platinum was obtained from the information for deuterons in lead, as given in Reference 5. To convert the range of deuterons in lead to the range in platinum the following formula was used. Ro =R Zpb Apt Pt Pb A Zpt (113) Here R is the range in mg/cm2, Z is the atomic number, and A is the atomic weight associated with the elements represented by the subscripts. In the first bombardment, the energy of the deuteron beam striking the second platinum foil is unknown since both the first platinum foil and the 15.6 mil copper absorber showed signs of being burned.

-201Pt Cu Pt --- —.-19 Mev —... 6 to 10 Mev~ b —----------—.._bt / 1.05 mil 15.6 mil 1.05 mil / / Cu Pt C~i Pt Deuteron B 81 \ 6.7 mil 1.05 mil 3.7 mil 1.05 mil Pt, Pt --—.t.- 19.2 ev — 4 17.9 Mev — 49.6 Mev _-3wl bdt 1.05 mil 1.05 mil 9.4 mil 1.05 mil Figure 65. Foil arrangement for the three natural platinum bombardments After the chemical separation, samples were mounted on 1/4 mil Teflon film. Decay curves were taken on 4ic n-ray counters, and on 7-ray counters. Since this work was performed before the hollow scintillator P-ray spectrometer was built, data relating to P-ray spectra were obtained in two ways. The first involved the use of the 1800 magnetic spectrometer described by Meinke, et. al. (97), and the second was the aluminum absorption curve method. Gamma-ray data were obtained primarily through the use of oscilloscope photographs, although Speedomax sweeps were also taken. Coincidence data were obtained manually, since the photographic method had not as yet been developed when this work was being done.

-2023. Results Obtained From Natural Platinum Bombardments The decay curves of eight separate samples obtained at various bombarding energies show the presence of only three components with halflives of 19 hours, -8 days and 75 days. The 19-hour activity is due to Ir9, while the 75-day activity comes from Ir192. The latter two isotope assignments were substantiated by n-and 7 -ray data presented later. A half-life of 8.3 + 0.5 days was obtained for the intermediate lived activity by averaging the values obtained from eight decay curves. This may be compared with the value of ^ 9 days found by Butement and Poe. Most of the above decay curves were begun'"16 hours after bombardment, and were followed for 3 months. While activities with halflives longer than 75 days might have been present, any appreciable amount of a half-life longer than about 5 hours and less than 75 days would have been detected. The presence of only three half-lives indicates that the chemical separation was adequate, and that the 8.3-day activity was indeed an isotope of iridium as claimed by Butement. This is further confirmed by the fact that repeating the chemical separation did not change the ratio of the three activities. Finally, the relative cross-section values for the 8.3-day activity obtained on four different samples at bombardment energies of 19.2 and 20.4 Mev were all equal to within experimental error. 194 Rough cross section values for the (d,y) reaction on Pt9 and Pt196 to produce Ir192 (75 day) and r196 (19 hour) were obtained using Equation 114.

-203(Tl/2 hr) (Co/m) (Mwt) a = 1.117 x 10-13 x (abund. ) (Q) (Chem Yld.) (Aliquote) (114) where a = The cross-section in barns (10-24cm2) T1/2 hr = The half-life in hlurs Co/m = The activity in counts per minute at "zero" time. Mwt = Mass of target nuclei Abund. = Fractional abundance of target nuclei Q = Deuteron beam current in A amp. hours Assuming for the moment that the 8.3-day activity does belong to 196 Ir 96 it will be shown later that this is not true), the calculated crosssections appear in Table XXIII. The errors are estimated standard deviations. The errors in beam energy were estimated from the variations in foil thicknesses and from the initial beam distributions. The chemical yield were obtained by weighing as IrBr4 an aliquote of the final iridium fraction. No corrections were made for the self absorption of the D rays in their sources. Within the stated errors, self-absorption should only affect the "Ir 196 cross-sections, and it will be shown later that these values must be discarded anyway since they do not actually refer to Ir96.

-204TABIE XXIII. CROSS-SECTIONS FOR THE (d, a) REACTION ON PLATINUM Target I_____ Energy of Deuterons in Mev larger Nuclei 9.6 * 1.3 11.3 T 1.6 15.9 - 1.5 19.2; 1.6 20.4 t 0.8 3-4'1-0 x10-4 Pt 194 l. 2+0.4xl05 l.lto.5xlO'4 1.7+.8x1O- 3.8t1.2xl04 2.1~0.9x104 2.0t0.6xl0'4 196 +4~ + 4 + 4 7.2~3.8x10I. Pt6 4-o1.6xoi.6xO' 331x 4.4+1.5xlO 2. lt.9x 3.31 2...O. 9xlO 1.9.0.7x10'4 The following information was obtained using the magnetic 3 -ray spectrometer: TABIE XXIV. BETA-RAY SPECTROMETER DATA FROM NATURAL PLATINUM BOMBARDMENTS ays After Bombardment Beta-Ray Information 3 No + Maximum O energy t 2 Mev. Lower 3 -component energy, 0.7 Mev. 8 Very little 2-Mev 3 component left. Maximum 3 energy of principle component x 0.7 Mev. Indication of low energy 1 component s0,~07 - 0.1 Mev. 10-15 | -0.08 Mev 3 component decays faster than the 0.7-Mev component

-205Figure 66 shows the n-ray spectrum obtained on a typical sample 8 days after bombardment. The above information is in accord with the literature values of Emax for Ir92 and Ir94 which emit rays with maximum energies of 0.72 and 2.2 Mev respectively, and with Butement's observationof an " 0.08-Mev ray 196 ascribed to Ir. Aluminum absorption curve data confirmed the above spectrometer data. By following the decay of a portion of the P-ray spectrum at about 0.05 Mev it was found that the lowest energy group was decaying with a half-life of l0 + 2 days. The decay was not followed long enough to define the half-life more exactly. Gamma-ray decay curves were obtained which showed the presence of the 8.3-day activity. Oscilloscope photographs and Speedomax sweeps showed the presence of the following 7 rays which could not be ascribed to either Ir92 or r194 TABIE XXV. GAMMA RAYS FROM 8.3-DAY IRIDIUM ACTIVITY'' m - Mev 0.064 0.378 0.097. 407 0.1-0.29 0.448 0.341 0.548 ---— J 3 _ _ -: -: - Z --- - S -.r -I::. - -: -- - - -. - - S Y ~:: -. S or |_ - - _ The energy range from 0.1 - 0.29 Mev was obscured by the Compton smear from higher energy 7 rays, and no accurate data could be obtained at these energies. The relative intensities of the above r rays could not be accurately estimated due to the intense background of similar energy y rays arising from Ir9. Most of the coincidence data-that were obtained proved inconclusive again due to the interference from Ir 92. All that can be said is that the 0.34-Mev 7 ray is in coincidence with 7 rays in the range 0.40-0.45 Mev,

r m I I I I I I I I I 1 600- BETA-RAY SPECTRUM Hi OF Ir ACTIVITIES 500 LLz -400 (I) 0 o200 100 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 ENERGY (MEV) Figure 66. Beta-ray spectrum of iridium activities eight days after bombardment

-207and that the 0.064-Mev transition (which corresponds in energy to a Pt x-ray) is in coincidence with itself. To summarize, the above results correlate well with those obtained 196 by Butement and Poe who ascribed the 8.3-day activity to Ir 96 Even the rough cross-sections seemed to indicate that their assignment was correct. However, before a decay scheme could be formulated, better coincidence information was needed. This could not be obtained in the presence of a 192 significant amount of Ir9. Therefore, it was decided to obtain a sample of platinum enriched in the isotope Pt, the parent of Ir9 in the (d,C ) reaction. 4. Results Obtained From the Enriched Platinum Bombardment Table XXVI contains the analysis of the enriched platinum target. 198 194 With this target the ratio of Pt 9/Pt was increased by a factor of 100 over the natural ratio. Fifty mgs. of the powdered platinum metal were bombarded with 23.5-Mev deuterons at the University of California cyclotron. Six days after bombardment the sample was worked up and a decay curve started. TABIE XXVI. ANAIYSIS OF ENRICHED PLATINUM ISOTOPES Mass Analysis Normal Isotope Atomic % Precision Atomic % 190 <0.01 0.012 192 0.042 0.013 0.78 194 3.57 0.10 32.8 195 8.97 0.28 33.7 196 26.47 0.78 25.4 198 60.95 0.72 7.2

-208TABIE XXVI. ANAIYSIS OF ENRICHED PLATINUM ISOTOPES (Cont'd) Spectrographic Analysis Element Element % Al <0.05 Mo <0.02 Ba < 0.02 Na <0.02 Be < 0.001 Nb <0.1 Ca <0.2 Ni <0.05 Co <0.05 Pb <0.1 Cr <0.05 Pd <0.04 Cu <0.05 Rh <0.04 Fe <0.02 Ru <0.1 Ir <0.1 Si <0.05 K 0.01 Ti <0.02 Li <0.01 V'0.02 Mg 0.02 Zr <0.1 Mn <0.02 Only a very small amount of activity was found in the iridium fraction after bombardment. Gamma-ray data showed that the primary activity was still Ir 92, with the addition of a small amount of 7 radiation in the 0.7-0.9 Mev range. A large amount of 7-7'coincidence information was obtained, but all of it could be explained by the presence of 192 192 194 Ir. An excellent discussion of the decay of Ir 9 and Ir9 is given by Johns and Nablo in Reference 68. The decay curve for this bombardment was followed for 50 days. 192 When resolved it showed only two activities to be present, the 75-day Ir and a small amount of a 3-4 day activity which probably was a mixture of 196 199 Au96 and Au 99 No indication of an 8.3-day activity was found at all. The atom yield of the 3-4 day activity was less than the atom yield for Ir by a factor of 6. In a normal platinum bombardment where the Ir

-209activity would be increased by a factor of 10, the small amount of the 3-4 day activity would have escaped attention. 5. Conclusions While the natural platinum bombardments strongly supported Butement and Poe in their assignment of the 8.3-day activity to Ir 196 the enriched isotope bombardment completely disqualified this assignment. The results of this work, however, indicate that the 8.3-day activity does belong to some isotope of iridium. The following reinterpretation of the data was made. The presence of a relatively large amount of platinum x-rays shows that the 8.3-day activity decays primarily by electron capture, with little or no positron emission as indicated by the n-ray spectrometer data. This assumption allows the assignment of the ~0.08-Mev 3 particles to fluorescent electrons, or conversion electrons of a low energy Y ray. This last is a logical assignment that removes the difficulty of demanding that the 3 decay be allowed, as would be necessary if Emax = 0.08 Mev and the half-life were 8 days. Decay by electron capture requires that the iridium isotope be 189 190 on the neutron deficient side of stability. Here Ir and Ir appear as possibilities, since each has an isomer with a half-life of the right order of magnitude. Neither of these isotopes has been well characterized, but recently Aten and co-workers (6) have presented some information on Ir9, and Smith and Hollander (134) have reported on Ir 9. Table XXVII lists their findings.

-210TABIE XXVII. GAMMA-RAY ENERGIES IN MEV FROM Ir189 AND Ir190' —. _,'''_- -. -- -' -- -' -: —: — _*" - L-.' i189 190 Ir Ir x-ray x-ray 0.135 0.186 0.245 0.360 0.400 0.540 0.580 0.800 weak 1.33 weak The y rays found in this work permit the assumption that the observed 8.3-day activity is the result of the presence of both Ir89 an 190 189 Ir. Ir could not be obtained by the (d.,) reaction on platinum, and 190 Ir would have to come from Pt192 with a natural isotopic abundance of 190 0.78%. Furthermore Ir has a 3-hour isomer which would further reduce the amount of 8.3-day activity. In the enriched platinum target, Pt192 was reduced to 0.042%, almost a factor of 100 less than Pt194 which produces the 75-day Ir192 activity. It is quite understandable that in the enriched isotope bombardment no 8.3-day activity was found. To account for the 8.3-day activity found in the natural platinum bombardments, it is suggested that a small amount of osmium impurity was in the platinum targets. This would produce iridium activities by the 189 190 favorable (d,n) reaction. Hence both Ir and Ir could be produced from Os and Os with natural abundances of 13.3 and 16.1% respectively.

-211Os1 (26.4%) and Os (41.0%) both yield stable iridium isotopes. Osl86(1.6%) 187 would produce a small amount of the 12-hour Ir which would not have been resolved from the more intense 19-hour Ir. Os 7 (1.6%) would yield 188 189 some 41-hour Ir, which would be masked to a certain extent by the Ir 190 and the Ir activities. A re-examination of the decay curves did indicate that a small amount of the 41-hour Ir might have been present. Since the (d, n) reaction is more favorable than the (d, a) reaction by a factor of 2 3 10 -103, the amount of osmium impurity needed to account for the results would only be 0.01 to 0.1%. The spectrographic analysis of the natural platinum indicates that this assumption is not unreasonable. Furthermore, Chu (25) 189 190 has produced Ir and Ir by deuteron bombardment of osmium. Concerning the true Irl96 activity, the present work sugges-t that an upper limit for the half-life would be about 5 hours. Actually, a much shorter half-life would be anticipated, perhaps in the range of minutes. 196 That Ir does not decay with a half-life of 8.3 days eliminates a situation that would have been difficult to explain using present theories of nuclear structure. 47 D. The Isotope Sc In 1953 several reports (24,27,87) appeared concerning the decay 47 of the 3.4-day isotope Sc. In each report two P-ray transitions and one 7 ray transition were found. However, the energies of these transitions varied from report to report. In Table XXVIII these values are listed. It therefore was of interest to make a careful study of this isotope, perhaps with the aid of enriched isotopes, to better determine its decay characteristics. While this work was in a preliminary stage, an excellent article

-212TABLE XXVIII. RADIATIONS EITTED DURING THE DECAY OF Sc47 Reference P-ray.-ray 7-ray y-ray Energy(Mev) Abund.(%) Energy(Mev) Abund(%) 24 0.435 O 66 0.185 100 0.622 / 34 27 0.64 0.160 87 0.28 28 fd 0.22 0.49 72 82 0.450 74 0.160 100 0.610 26 by Lidofsky and Fisher (82) appeared describing the decay of Sc4. Since it did not appear likely that their results could be improved upon with the equipment that was available in this laboratory, the project was dropped. The chemical separation that was developed in the present work, and the results of four bombardments on TiO2 targets are given below. 1. Chemical Separation The radioactive scandium was produced by the bombardment of TiO2 targets using the (d,a) reaction. The chemical separation must therefore isolate scand-ium in a carrier-free form from titanium and vanadium activities also produced during the bombardment. The first problem to be encountered was to find a way to dissolve the TiO2. Attempts to dissolve TiO2 in concentrated H2S04 plus H202 were not successful, since several days were usually required to obtain a solution. It was decided then to try a fusion process.

-213Using K2S207, obtained by heating K2S208 until So3 fumes ceased to evolve, the fusion process required only 10-15 minutes. After the flux cooled below e-2000C, H2S04 was added. The flux was then reheated.until all solids dissolved in the H2S04. Upon cooling, the flux would remain liquid which facilitated its removal from the crucible. The first attempts at a chemical separation were along the lines that Hall (51) was following in the development of his scandium separation chemistry. In this case the scandium activity was to be separated in a carrier-free form by filtering out the scandium radiocolloid on filter paper from a solution of pH 8.5. The titanium would be kept in solution by forming the complex with H202. Very little success was found with this method. The titanium was difficult to keep in solution and was constantly appearing as a dark orange-brown compound on the filter paper. Furthermore, the scandium colloid had a tendency to stick to the walls of the glassware, and also was difficult to remove from the filter paper. It seemed advisable to try another approach. Ion-exchange techniques were to be the basis of the new separation, since Ti (IV) and V (V) will adsorb on Dowex-2 resin in 11-12 N HC1 while scandium will not (78). The difficulty at this stage is that V (IV) does not adsorb on the column, and V (V) adsorbed on the resin will tend to reduce to V (IV) in concentrated HC1 solutions. After trying several approaches it was found that a small amount of KC103 in the concentrated HC1 solutions will keep V(V) from reducing. In fact, the V (V) adsorbs so well that both scandium and titanium may be removed from the column with 11 N HC1 and the V removed later with 6 N. HC1, providing a clean separation for vanadium.

-2l1The details of the scandium chemical separation are given in Table XXIX. It will be noted that in steps 6-10 a rough separation of scandium from titanium and vanadium is made prior to the use of the ionexchange procedure. This was done primarily to insure that no titanium contamination would be found in the scandium fraction. The final scandium fraction may contain a small amount of NaCl. If necessary the scandium may be separated from the salt by slow elution on an exchange column using concentrated. HC1. The salt will come off first. TABLE XXIX. SCANDIUM SEPARATION CHEMISTRY Element separated: Scandium Procedure by: Gardner Target Material: Titanium Oxide Time for sep'n:' 2 hours Type of bbdt: 7.8 Mev Equipment required: Ion deuterons exchange column of AG 2-X8, 200-400 mesh resin obtained from Bio-Rad Laboratories. Resin bed was 8 mm x 140 mm. Yield: > 60c 6 5 Degree of purification: >10 from V; 10 from Ti Advantages: Carrier-free separation with high decontamination Procedure: (1) Prepare column by washing several times with conc. HCU, then with conc. HC1 + KCIO (See remark 1). (2) Place TiO2 in a small porcelain crucible and add roughly 1-2 gms K2S O per 100 mgs of TiO2. Cover and bring to dull read heat. Continue heating until TiO2 dissolves in flux (10-15 min) (See remark 2). (5) Cool below 2000C. and cautiously add 1-2 ml of conc. H2S04.

-215TABIE XXIX. SCANDIUM SEPARATION CHEMISTRY (Cont'd) (4) Heat until all solids dissolve in the H2SOL. Upon cooling, contents of crucible should remain liquid. If solidification occurs, add more conc. H2SO4 and repeat heating. (5) Pour liquid into a beaker using a small amount of H20 as a wash. (6) Add solid NaOH until solution is basic to litmus and TiO2 precipitates out as a white hydrous oxide. Stir thoroughly for a few minutes to insure that all Sc will be carried down on the TiO2. (See remark 3) (7) Centrifuge out the TiO2 and salt. (Most of the V activity remains in the supernate). (8) Wash the TiOp and salt several times with 1 N NaOH until all of the salt has dissolved. (See remark 4). (9) Dissolve TiO2 in a few drops of cone. HC1 and reprecipitate with saturated Na2CO3 solution. (Much of the Sc goes into the supernate in this step.) (10) Remove supernate and repeat step 9 twice more, saving the supernatents at each step and combining them in the end. (See remark 5). (11) Total supernate is acidified with conc. HC1 and evaporated down to less than 20 ml. Solution is then cooled in an ice bath and saturated with HC1 gas. (Saturation is assumed complete when no more salt comes out of solution.) Centrifuge and discard salt. (12) Continue evaporation almost to dryness. Extract residue with 5 ml of 11 M HC1 three times. Evaporate HC1 solution to 5-10 ml, add - 5 mgs of KC103, cool in ice, and saturate with H(31 gas. Centrifuge to remove any salt and then place on column prepared in step 1. (13) Force solution through column at rate of 1-2 drops/sec., stopping just before liquid reaches top of resin. (Do not allow top of resin to ever become dry ) Some Sc will appear after first 5 ml. Elute with 11 M HC1 + KC103 at rate of 1-2 drops/sec. (See remark 1) (14) Most of Sc comes off in first 10-15 mls of 11 M HC1. Ti begins to come off after 25-30 mls of 11 M HC1. (See remark 6) (15) A very small amount of salt may follow the Sc. This may be removed by evaporating eluent containing the Sc down to 5-10 ml, saturating it with HC1 gas, and then running it through the column washed with just conc. HC1 - at a much slower rate. The salt comes off first.

-216TABIE XXIX. SCANDIUM SEPARATION CHEMISTRY (Cont'd) Remarks: (1) All HC1 + KCLO3 solutions contain just enough KC10O to impart a faint yellow color to the solution. Usually this was 0.5 mg KC103. ml (2) K2S207 was made by heating K2S208 until SO3 fumes ceased to evolve. (3) Some water may be added along with the solid NaOH, but the volume of solution should be kept small, usually 20 ml or less (4) All of the washings in step 8 will usually contain much less than 10% of the total Sc activity. (5) Of the total Sc activity in the TiO2 precipitate over 60% is removed in the first reprecipitation and over 20% is removed in the second. (6) The eluent (including original solution) is usually taken in separate units of 3-5 ml. Aliquots of these are counted to give the Sc and Ti peaks, then the center 60-70% of Sc peak is taken for analysis. (7) General Reference for this tyee procedure: Hicks, H. G., et al., "The Qualitative Anionic Behavior of a Number of Metals with an Ion Exchange Resin,'Dowex 2' ", Livermore Research Laboratory Report, LRL-65, Dec. 1953. 2. Bombardment Results. Four bombardments in all were obtained using natural TiO2 powder targets. The TiO2 was listed as grade C.P. by J. T. Baker Chemical Company, Phillipsburg, New Jersey. The bombardments were obtained to produce radioactive sources of titanium, vanadium, and scandium to be used in the development of the chemical procedure. In the fourth bombardment the target was inadvertantly destroyed in the cyclotron.

-217Decay curves taken on the scandium fraction from the third bombardment were followed for three months. These showed the presence 46 of the 85-day Sc and also a shorter-lived component of r3 days. This latter activity was due to a combination of the following activities: Sc44 (2.4 day), Sc47 (5.43 day), and Sc48 (44 hour). No trace of any activities that could be assigned to vanadium or titanium was found. Oscilliscope photographs showed the presence of five y rays. These are listed in Table XXX with their identification. TABIE XXX. GAMMA RAYS FROM SCANDIUM FRACTION 7-ray energy Isotope Literature value (Mev)' of 7-ray energy. Mev. 0.16 Sc 0.16 47 Sc 0.16 0.28 Sc4m 0.27 44 0 52 Sc annihilation radiation 0.511 48 0.9-0.97 Sc 0.98 46 Sc 0.89 1.1-1.3 Sc44 1.16 46 Sc 1.12 Sc48 0.99, 1.04, 1.33

-218No 3-ray spectral data were obtained from any bombardments. In summary, it may be said that a chemical separation was developed that will separate scandium and vanadium in carrier-free form with good decontamination from a titanium target. A 7-ray with an energy of 0.16 47 Mev was found that could be assigned to the isotope Sc. The 7-ray r47 results are in agreement with the latest evidence on the isotope Sc

APPENDIX A Fabrication of Hollow Plastic Scintillator The first attempts at constructing a hollow n-ray detector involved drilling and grinding out a cone-like cavity in the center of a solid plastic scintillator. The grinding was done using dentist gem stones. After the depression had been roughed-out, the walls of the phosphor were polished with crocus cloth and rouge on felt. The phosphor was then wrapped with 1.5-mil aluminum foil. Two detectors were made using the above process. One of these was in the form of a right circular cylinder, 3/4" in diameter and 1/2" high, while the other detector was a truncated cone 1" in diameter and 3/4" high. Both of these detectors proved unsatisfactory, primarily due to cracks, which formed from strains set up in the plastic during the hollowing-out process. These cracks reduced the resolving power below that of a solid detector. Besides, it proved almost impossible to produce a symmetrical, conical cavity with a reasonably smooth bottom within the plastic and yet retain a small entrance hole in the top. Polishing the interior surfaces was also difficult. All attempts at solvent polishing, using kerosene, acetone, alcohol, or benzene, failed. The usual result produced by these solvents is an opaque, soft, plasticized layer formed on the surface. The following procedure was devised to produce satisfactory hollow detectors: Each detector was constructed in two parts which consist of a flat, solid, right circular cylinder used as the base, and a top piece consisting of a hollowed-out cylinder or truncated cone. -219

-220The cross-section views for the three detectors were given in Figure 13. As previously mentioned, the three detectors are to be used with maximum p energies of < 1.5 Mev, <2.2 Mev, and c 3.6 Mev. The pieces of phosphor were machined on a lathe from a solid piece of stock rod. The plastic melts at 110 C., so the machining should be done slowly while a stream of air is played over the piece. Every effort must be made not to introduce strains into the plastic. A sharp, finely pointed tool should be used as the plastic has a tendency to chip. It is not necessary to anneal the pieces if care is taken when machining. No cracks have appeared in the pieces machined in the air stream. Since each detector is constructed in two pieces, the polishing of the pieces is greatly facilitated. Particular care should be given the polishing of the top and bottom surfaces of the base, and the inside and bottom surfaces of the hollow piece. The point in polishing the interior surfaces of the detector is that light generated in the phosphor will escape from a rough surface more easily than from a polished surface. The polishing is accomplished using increasingly fine grades of sandpaper, ending with crocus cloth. The final finish is produced using Bon Ami on felt or soft paper tissue. If the initial polishing is done on a lathe, care should be taken not to burn or melt the plastic. Wet sandpaper may be used if desired when the polishing.is done by hand. The detectors are assembled and mounted in the following way: 1. Wipe each piece to remrove any dirt or grease. 2. Apply a small amount of Dow-Corning 200 Silicone fluid to the bottom surface of the hollow section. Wipe off excess fluid

-221 3. Place the hollow section onto the base. Press together until excess fluid is squeezed from between the pieces. Make sure that there are no air bubbles in the junction. Just enough fluid should be used so that only a very small amount is squeezed out. Wipe off excess fluid from outside of detector. 4. Remove excess fluid from inside of piece with tissue paper wrapped around a thin dowel rod. Be careful not to scratch plastic. The fluid must be carefully removed since the ( particles will loose energy passing through it. that will not appear as light. 5. Cut a circular piece of Reynold's Wrap slightly larger than the top surface. Cut another piece of aluminum foil to approximate the shape of the outside surface of the phosphor. With the shiney side of the foil toward the phosphor, wrap the detector with the foil leaving the bottom face uncovered. 6. Cut away the portion of foil covering the entrance hole in the detector. 7. Tape the foil-enclosed phosphor to the aluminum collar shown in Fig. 16. 8. Mount detector on phototube with Dow-Corning 200 Silicone fluid. 9. Place phototube in preamplifier housing.

APPENDIX B Energy Units In the Beta Decay section of Chapter II reference was made to D energies measured in relativistic (rest mass) units. This unit is somewhat more convenient than Mev (million electron volt) units when spectra are to be corrected for forbiddenness. The above units are related by the following equation. W= - +1 (8) mc where W is energy in relativistic units, E is energy in Mev, and mc is the rest mass of the electron (0.5109 Mev). Using Figure 67 the two units may be interchanged easily. Relativistic Mev Units 01 o.5- -2 1.0 —3 1.5- -4 2.0- -5 2.5- -6 3.0- -7 3.5- -8 4.0- -9 Figure 67. Relation Between Mev and Relativistic Units -222

BIBLIOGRAPHY 1. Aebersold, P. C., International Conference on Peaceful Uses of Atomic Energy, Geneva, Paper 308 (USA) (1955). Importance of isotopes in technology and industry. 2. Aebersold, P. C., Rupp, A. F., Nucleonics 10, No. 1, 24-7 (1952). Radioisotope production in the United States. 3. Ageno, M., Nuovo cimento (9) 1, 33-40 (1943). Radioattivita provocata nello zirconio da bombardamento di neutroni veloci. 4. Ames, D. P., Bunker, M. W., Starner, J. W., U. S. Atomic Energy Commission, reported in Los Alamos Scientific Laboratory Classified Report WASH-75 (January 1952). 5. Aron, W. A., Hoffman, B. G., Williams, F. C., U.S. Atomic Energy Commission, Rept. AECU-663 (May 1951). Range-energy curves. 6. Aten, Jr., A. H. W., DeFeyfer, G. D., Sterk M. J., and Wapstra, A. H., Physica 21, 740 (1955). The isomers of Osl0O and Irl90. 7. Beaufait, Jr., L. J., Lukens, Jr., H. R., U.S. Atomic Energy Commission, Rept. NP-5056 (May 1953). Handbook of radiochemical analysis. Radiochemical techniques. 8. Behrens, D. J., Brit. Atomic Energy Research Establishment, Rept. AERE-T/R-629 (January 1951). Some notes on the nature of experiments, on the statistics of counting, and on the fitting of exponential decay curves to the results of counting experiments. 9. Bernard, R., Davoine, F., Hirtz, J., Compt. rend. 232, 1826-8 (1951). Etudes de couches minces adsorbees au moyen d'el6ments traceurs radioactifs. 10. Bernstein, W., Chase, R. L., Schardt, A. W., Rev. Sci. Instr. 24, 437-44 (1953). Gray wedge pulse-height analysis. 11. Birks, J. B., "Scintillation Counters," McGraw-Hill, New York, 1953. 12. Blatt, J. M., Weisskopf, V. F., "Theoretical Nuclear Physics," Wiley, New York 1952. 13. Bleuler, E., Goldsmith, G. J., "Experimental Nucleonics," Rinehart, New York, 1952. 14. Bohr, N., Nature 137, 344-8 (1936). Neutron capture and nuclear constitution. 15. Bohr, N., Wheeler, J. A., Phys. Rev. 56, 426-50 (1939). The mechanism of nuclear fission. -223

-22416. Born, H. J., Seelman-Eggebert, W., Naturwissenschaften 31, 201-2 (1943). Die Maximalenergien der Beta-Strahlen einiger bei der Uranspaltung entstehender Isotope. 17. Bracewell, R. N., J. Opt. Soc. Amer. 45, 873-6 (1955). Simple graphical method for correcting for instrumental broadening. 18. Brookhaven National Laboratory, U.S. Atomic Energy Commission, Rept. BNL-364 (AS-9) (July 1955). Annual report. 19. Butement, F. D. S., Poe, A. J., Phil. Mag. 45, 31 (1954). Radioactivity of neutron excess iridium isotopes. 20. Cassatt, W. A., Ph.D. Thesis, University of Michigan, 1954. Decay scheme and characterization study of radioactive products of deuteron reactions. (Also AECU-2958) 21. Cassatt, W. A., Meinke, W. W., Phys, Revs. 99, 760 (1955). Decay Scheme of Y92. 22. Chase, R. L., U.S. Atomic Energy Commission, Rept. BNL-263 (T-42) (April 1954). Coincidence gray wedge pulse height analyzer BNL Model EH1-501. 23. Christian, D., Mitchell, R., Martin, Jr., D., Phys. Revs. 86, 946 (1952). Radioactivities of Pt and Ir from photonuclear reactions in platinum. 24. Chenz, et al., Phys. Rev. 90, 886 (1953). Tle disintegration of Sc47. 25. Chu, T., Phys. Rev. 79, 582 (1950). New radioactive isotopes of iridium. 26. Cook, Jr., T. B., U.S. Atomic Energy Commission, Rept. AECU-1598 (July 1951). The scintillation spectrometer as a means of measuring conversion coefficients. 27. Cork, J. M., et al., Phys. Rev. 92, 367 (1953). The radioactive decay of calcium 47. 28. Cornell, R. G., U.S. Atomic Energy Commission, Rept. ORNL-2120 (eptember 1956). A new estimation procedure for linear combinations of exponentials. 29. Coryell, C. D., Sugarman, N., ed., "Radiochemical Studies: The Fission Products," National Nuclear Energy Series IV, Vol. 9, McGrawHill, New York, 1951. 30. Ibid., p. 660-2. 31. Ibid., p. 1410-16.

-22532. Curran, S. C., "Luminescence and the Scintillation Counter," Academic Press, New York, 1953. 33. Demichelis, F., Malvano, R., Nuovo cimento (9) 11, 49-52 (1954). Discriminatore differenziale oscillografico e sua applicazione nei circuiti a coincidenza. 34. Deming, W. E., "Statistical Adjustment of Data," Wiley, New York, 1946. 35. DerMateosian, E., Goldhaber, M., Phys. Rev. 82, 115-6 (1951). Measurement of isomeric transition energies with a scintillation spectrometer. 36. Dhar, S., Indian J. Phys. 29, 329-51 (1955). On the energy response and resolution of a scintillation counter. 37. Dunworth, J. V., Rev. Sci. Instr. I1, 167 (1940). The application of the method of coincidence counting to experiments in nuclear physics. 38. Faltings, V., Angew. Chem. 64, 605-6 (1952). Eine Anwendung der natuirlichen 14C-Aktivitat in der Lebensmittelchemie. 39. Feenberg, E., Hammack, K. C., Phys. Rev. 75, 332 (1949). An interpretation of nuclear shell structure. 40. Feenberg, E., Hammack, K. C., Phys. Rev. 75, 1877-1893 (1949). Nuclear shell structure. 41. Feenberg, E., Trigg, G., Revs. Modern Phys. 22, 399-406 (1950). The interpretation of comparative half-lives in the Fermi theory of betadecay. 42. Filon, L. N. G., Proc. Roy. Soc. Edin. 49, 38 (1928-29). On a quadrature formula for trigometric integrals. 43. Finston, H. G., "Radiochemical Separations Compilation," Atomic Energy Commission, National Research Council, and Brookhaven National Laboratory, 1954. 44. Flanary, J. R., Culler, F. L., International Conference on Peaceful Uses of Atomic Energy, Geneva, Paper 539 (USA) (1955). Solvent extraction separation of uranium and plutonium from fission products. 45. Freedman, M. S. et al., Rev. Sci. Instr. 27, 716 (1956). Correction for phosphor backscattering in electron scintillation spectrometry. 46. Friedlander, G., Kennedy, J. W., "Introduction to Radiochemistry," Wiley, New York, 1949. 47. G6tte, H., Naturwissenschaften 29, 496-500 (1941). Neue bei der Uranspaltung auftretende Strontium- und Yittrium-Isotope.

-22648. Griffith, G. L.4 Phys. Rev. 103, 643-5 (1956) First excited states in Zr92 and Zr 9. 49. Grodzins, L., Rev. Sci. Instr. 26, 1208-9 (1955). Coincidence sorter for scintillation spectrometers. 50. Gardner, D. G.,Meinke, W. W., "Decay Scheme of Co62" to be published. 51. Hall, K. L., Ph.D. Thesis, University of Michigan, 1955. Determination of (d, alpha) reaction cross sections. (Also AECTZ-e126). 52. Halliday, D., "Introductory Nuclear Physics," Wiley, New York, 1955. 53. Hartree, D. R., "Numerical Analysis," Clarendon Press, Oxford, 1952. 54. Hayward, R. W., Hoppes, D. D., Ernst, H., Phys. Rev. 98, 231 (1955). Disintegration of niobium-91 and -92. 55. Hendricks, S. B., Dean, L. A., Ann. Rev. Nuclear Sci. 1, 597-610 (1952). Radioisotopes in soils research and plant nutrition. 56. Heusinkveld, M., Freier, G., Phys. Rev. 85, 80-4 (1952). The production of polarized protons and the inversion of energy levels of the P 1/2 - P 3/2 doublet in Li5. 57. Hicks, H. G., Gilbert, R. S., Stevenson, P. C., Hutchin, W. H., U.S. Atomic Energy Commission, Rept. LRL-65 (December 1953). The qualitative anionic behavior of a number of metals with an ion exchange resin, "Dowex 2." 58. Hildebrand, F. B., "Introduction to Numerical Analysis," McGraw-Hill, New York, 1956. 59. Hill, D. L., Wheeler, J. A., Phys. Rev. 89, 1102-45 (1953). Nuclear constitution and the interpretation of fission phenomena. 60. Hollander, J. M., Perlman, I., Seaborg, G. T., Revs. Modern Phys. 25, 469-73 (1953). Table of isotopes. 61. Householder, A. S., "Analyzing Exponential Decay Curves," Seminar in Scientific Computation, Nov. 1949 62. Householder, A. S., U.S. Atomic Energy Commission, Rept. ORNL-455 (February 1950). On Prony's method of fitting exponential decay curves and multiple-hit survival curves. 63. Hunt, W. A., Zaffarano, D. J., U.S. Atomic Energy Commission, Rept. ISC-469 (March 1954). Radioactive disintegration spectra of some short-lived nuclides. 64. Irvine, Jr,, J. W., Nucleonics 3, No. 2, 5-17 (1948). Production of radionuclides.

-22765. Jahnke, E., Emde, F., "Tables of Functions with Formulae and Curves," Dover Publications, New York, 1943. 66. James, R. A., Phys. Rev. 93, 288-90 (1954). Excitation functions of proton-induced reactions of Nb93*. 67. Johansson, S.A.E., U.S. Atomic Energy Commission, Rept. ISC-431 (December 1953). A coincidence scintillation spectrometer. 68. Johns, M. W., and Nablo, W. V., Phys. Revs. 96, 1599 (1954). Disintegration of Ir192 and Ir194. 69. Jordan, W. H., Ann. Rev. Nuclear Sci. 1, 207-44 (1952). Detection of nuclear particles. 70. Katcoff, S., Miskel, J. A., Stanley, C. W., Phys. Rev. 74, 631-38 (1948). Ranges in air and mass identification of plutonium fission fragments. 71. Ketelle, B. H., Phys. Rev. 80, 758 (1950). The Ca45 beta-distribution obtained in a split crystal scintillation spectrometer. 72. Kleinberg, J., et al., U.S. Atomic Energy Commission, Rept. LA-1566 (February 1953).' Collected radiochemical procedures. 73. Kleinberg, J., et al., U.S. Atomic Energy Commission, Rept. LA-1721 (September 1954). Collected radiochemical procedures. 74. Kohl, J., Chem. Eng. Progr. 48, 611-6 (1952). Radioisotopes in progress instrumentation. 75. Konopinski, E. J., Langer, L. M., Ann. Rev. Nuclear Sci. 2, 261-304 (1953). The experimental clarification of the theory of beta-decay. 76. Kraus, K. A., Nelson, F., International Conference on Peaceful Uses of Atomic Energy, Geneva, Paper 837 (USA) (1955). Anion exchange of studies of fission products. 77. Kraus, K. A.,Moore, G. E., J. Am. Chem. Soc., 75, 1460 (1953). Ion exchange studies Vi. The Divalent transition elements manganese to zinc in hydrochloric acid. 78. Kraus, A. K., Nelson, F., Smith, G. W., J. Phys. Chem. 58, 11 (1954). Ion exchange studies IX. Absorbability of a number of metals in hydrochloric acid solutions. 79. Kraushaar, Brun, Meyerhof, private communication to K. Way in June, 1955. See reference 144, page 80. 80. Kurie, F. N. D., Richardson, J. R., Paxton, H. C., Phys. Rev. 49, 368-81 (1936). The radiations emitted from artifically produced radioactive substances. The upper limits and shapes of the beta-ray spectra from several elements.

81. Lidofsky, L. J., et al. Bull. Am. Phys. Soc. 30, #3 (Washington), 1953. 82. Lidofsky, L. J., et al., Phys. Rev. 104, 759 (1956). Isobaric level structure for A = 47. 83. MacGregor, M. H., Ph.D. Thesis, University of Michigan, 1953. The decay of rubidium 87. 84. Maeder, D., Helv. Phys. Acta 20, 139-52 (1947). Elektronenr5hrenSpektrograph zur Messung von Ionisationskammerimpulsen. 85. Mandeville, C. E., Schert, M. V., Nucleonics, 3, 2 (1948). Nuclear disintegration schemes and the coincidence method. 86. Margenau, H., Murphy, G. M., "The Mathematics of Physics and Chemistry," Van Nostrand, New York, 19566. 87. Marguez, L., Phys. Rev. 92, 1511 (1953). Disintegration of Ca47, Sc47, V, and Rh100. 88. Mayer, M. G., Phys. Rev. 74, 235-39 (1948). On closed shells in nuclei. 89. Mayer, M. G., Phys. Rev. 75, 1969-70 (1949). On closed shells in nuclei. 90. Mayer, M. G., Phys. Rev. 78, 16-23 (1950). Nuclear configurations in the spin-orbit coupling model, empirical evidence. 91. Mayer, M. G., Jensen, J. H. D., "Elementary Theory of Nuclear Shell Structure," Wiley, New York, 1955. 92. McLaughlin, P. W., O'Kelly, G. D., U.S. Atomic Energy Commission, Rept. MTA-40 (September 1953). Counting efficiencies for sodium iodide crystals. 93. Meinke, W. W., U.S. Atomic Energy Commission, Rept. AECD-2738, UCRL-432 (August 1949). Chemical procedures used in bombardment work at Berkeley. 94. Meinke, W. W., U.S. Atomic Energy Commission, Rept. AECD-2750 (December 1949). Chemical procedures used in bombardment work at Berkeley. 95. Meinke, W. W., Anal. Chem. 28, 736-56 (1956). Nucleonics. 96. Meinke, W. W., Anderson, R. E., Anal. Chem. 25, 778-83 (1953). Activation analysis using low level neutron sources.

-22997. Meinke, W. W., Cassatt, W. A., Hall, K. L.,. S. Atomic Energy Commission, Rept. AECU-2944 (June 1954). A variable field survey beta-ray spectrometer. 98. Mitchell, A. G. C., Rev. Mod. Phys. 20, 296 (1948). Use of coincidence counting methods in determining nuclear disintegration schemes. 99. National Bureau of Standards, Applied Mathematics Series 13 (1952). Tables for the analysis of beta spectra. 100. National Bureau of Standards, Applied Mathematics Series 34 (1954). Tables of the gamma function for complex arguments. 101. National Bureau of Standards Circular No. 499, U.S. Government Printing Office, Wash., D.C. 1950. Nuclear data including supplements and cards bringing table up to date. 102. Nordheim, L. W., Phys. Rev. 75, 1894-1901 (1949) On spins, moments, and shells in nuclei. 103. Nordheim, L. W., Revs. Modern Phys. 23, 322-7 (1951). Nuclear shell structure and beta-decay. Even A nuclei. 104. Noyes, A. A., Bray, W. C., "A System of Qualitative Analysis for the Rare Elements," The MacMillan Company, New York, 1948. 105. Nuclear Chemical Research, Progress Report 4, University of Michigan, 1955. 106. Nucleonics 8, No. 5, 5-11 (1951). Industrial utilization of fission products. 107. Nussbaum, R. H., Wapstra, A. H., VanLieshout, R., NiJgh, G. J. rnstein, L. Th. M., Physica 20, 571-2 (1954). On the decay of 6Zn, 2Cu and 2Co. 108. Oak Ridge National Laboratory, U.S. Atomic Energy Commission, Rept. ORNL-940 (December 1950). Physics division quarterly progress report. 109. Oak Ridge National Laboratory, U.S. Atomic Energy Commission, Rept. ORNL-1092 (June 1951, Physics division quarterly progress report. 110. Oak Ridge National Laboratory, U.S. Atomic Energy Commission, Rept. ORNL-1159 (July L951). Instrumental research and development quarterly progress report. 111. O'Kelly, G. D., Lazar, N. H. Eic ler, E., Phys. Rev. 101, 1059-63 (1956). Decay chain Ca49-Sc9-Ti 9. 112. Owen, G. E., and Primakoff, H., Phys. Rev. 74, 1406 (1948). Relation between apparent shapes of monoenergetic conversion lines and of continuous 5-spectra in a magnetic spectrometer.

-230113. Palmer, J. P., Laslett, L. J., U.S. Atomic Energy Commission, Rept. AECU-1220 (March 1951). Beta-ray spectrometry with an anthracene scintillation spectrometer. 114. Parmley, T. J., Moyer, B. J., Lilly, R. C., Phys. Rev. 75, 619-23 (1949). The radioactivities of some high mass isotopes of cobalt. 115. Perlis, A. J., U.S. Atomic Energy Commission Rept. NP-786 (September 16, 1948). Numerical solution of linear integral equations. 116. Putnam, J. L., International Conference on Peaceful Uses of Atomic Energy, Geneva, Paper 463 (UK) (1955). Development in thickness gauges and allied instruments. 117. Rapp, K., Wallack, S., Moos, A. M., U.S. Atomic Energy Commission, Rept. NYO-685, Pt. I and II (July 1950). 118. Rose, M. E., Science 118, 635-42 (1953). Angular correlation of nuclear radiations. 119. Rose, M. E., Goertzel, G. H., Perry, C. L., U.S. Atomic Energy Commission, Rept. ORNL-1023 (June 1951). K-shell internal conversion coefficients; revised tables. 120. Sagane, R., Kojima, S., Miyamoto, G., Ikawa, M., Phys. Rev. 57, 1179-80 (1940). Artificial radioactivity induced in zirconium and molybdenum. 121. Scarborough, J. B., "Numerical Mathematical Analysis," John Hopkins Press, Baltimore, 1950, p. 460. 122. Schardt, A. W., U.S. Atomic Energy Commission, Rept. BNL-237 (T-37) (January 1954). Gray wedge pulse height analysis. 123. Scharff-Goldhaber, G., Phys. Rev. 90, 587-602 (1953). Excited states of even-even nuclei. 124. Schott, G. L., Miipke, W. W., Phys. Rev. 89, 1156-7 (1953). Assignment of y92 and Y94. 125. Scott, E. J., "Transform Calculus with an Introduction to Complex Variables," Harper and Brothers, New York, 1955. 126. Seelmann-Eggebert,. Naturwissenschaften 31, 510 (1943). Uber einige aktive Yttrium-Isotope. 127. Shuttleworth, R., Brit. J. Appl. Phys. 4, 326-9 (1953). The use of radioactive isotopes in metallurgy. 128. Siegbahn, K., ed., "Beta- and Gamma-Ray Spectroscopy," Interscience, New York, 1-23 (1955).

129. Ibid., p. 132-64. 130. Ibid., p. 273-91. 131. Ibid., p. 314-56. 132. Ibid., p. 433-67. 133. Siri, W. E., "Isotopic Tracers and Nuclear Radiations; with Applications to Biology and Medicine," McGraw-Hill, New York, 1949. 134. Smith, W. G., and Hollander, J. M., Phys. Rev. 98, 1258 (1955). Radiochemical Study of neutron-deficient chains in the noble metal region. 135. Snedden, I. N., "Fourier Transforms," McGraw-Hill, New York, 1951. 136. Stevenson, P. C., Hicks, H. G., Ann. Rev. Nuclear Sci. 3, 221-234 (1953). Separation techniques used in radiochemistry. 137. Stotler, H. H., U.S. Atomic Energy Commission, Rept. AECU-99 (1946). Fitting in a straight line. 138. Sugarman, N., Phys. Rev. 86, 604 (1952). Meson fission of mercury. 139. Titchmarsh, E. C., "Introduction to the Theory of Fourier Integrals," Clarendon Press, Oxford (1937). 140. Titchmarsh, E. C., "The Theory of Functions," Oxford University Press, London, 1939. 141. Tranter, C. J., "Integral Transforms in Mathematical Physics," Wiley, New York, 1956. 142. Von Kgrman, T., Biot, M. A., "Mathematical Methods in Engineering," McGraw-Hill, New York, 1940. 143. Wahl, A. C., Bonner, N. A., "Radioactivity Applied to Chemistry," Wiley, New York, 1951. 144. Way, K., King, R. W., McGinnis, C. L., VanLieshout, R., U.S. Atomic Energy Commission, Rept. TID-5300 (June 1955). Nuclear level schemes A^40 - A^92. 145. Weiner, N., "The Fourier Integral and Certain of it Applications," Dover Publications, New York, 1951. 146. Whittaker, E. T., Watson, G. N., "A Course in Modern Analysis," Macmillan, New York, 1943.

-232147. Wiedenbeck, M. L., and Chu, K. Y., Phys. Rev. 72, 1164 (1947). Coincidence measurements. Part I, beta spectra; Part II, internal conversion. 148. Widder, D. V., "The Laplace Transform," Princeton University Press, Princeton, 1946. 149. Wilson, H. A., Phys. Rev. 77, 516-20 (1950). The spherical shell nuclear model. 150. Worthing, A. G., Geffner, J., "Treatment of Experimental Data," Wiley, New York, 1943. 151. Youden, W. J., "Statistical Methods for Chemists," Wiley, New York, 1951.