NEUTRON CROSS SECTIONS M. Ragheb 1/17/11 INTRODUCTION Neutron interactions with matter can be either scattering or absorption reactions. Scattering can result in a change in the energy and direction of motion of a neutron but cannot directly cause the disappearance of a free neutron. Absorption leads to the disappearance of free neutrons as a result of a nuclear reaction with fission or the formation of a new nucleus and another particle or particles such as protons, alpha particles and gamma photons. The probability of occurrence of these reactions is primarily dependent on the energy of the neutrons and on the properties of the nucleus with which it is interacting. We here consider the different reactions by which a neutron can interact with matter. Figure 1. Plutonium 9 actual total, elastic, inelastic, fission and capture cross sections.
Figure. Plutonium 9 group averaged total, elastic, inelastic, fission and capture cross sections. Neutron cross sections data are normally expressed in units of barns where -4 1[barn] = 1 [ cm ]. They are tabulated and plotted as a function of the kinetic energy of the neutrons in electron volts [ev] or million electron volts [MeV]. The actual measured cross section (Fig.1) displaying a large number of resonances or their evaluated group averages (Fig. ), are normally used in numerical calculations. Cross sections averaged over the appropriate neutrons energy range are needed for realistic computations. A common error in numerical applications is the use of cross sections that are not correctly defined over the appropriate energy range. MICROSCOPIC CROSS SECTION, REACTION RATE DENSITY We consider a thought or gedanken experiment to determine the reaction rate R [reactions/sec] that would occur in a small volume of a thin target material of area A [cm ] and thickness x [cm] when a neutron beam of neutrons moving in the x direction with a density n [neutrons/cm ] and velocity v [cm/sec] as shown in Fig.. If the density of the material is ρ [gm / cm ], and its atomic weight is M [amu] we can use a modified form of Avogadro s law to determine the nuclei density in the target as:
N ρ = Av (1) M Figure. Geometry for neutron reaction rate in thin target of volume V. The reaction rate R should be construed to be proportional to the area of the target A, its thickness x, the number density of the particles in the neutron beam, n, the velocity of the neutrons, v, and the number density of the nuclei in the target N. This can be expressed mathematically as: Reactions] R A.x.n.v.N [ () sec The proportionality symbol can be replaced by an equality sign provided we add a proportionality constant, leading to; Reactions R = σ.a.x.n.v.n [ ] sec () Since the volume of the target is: V = Ax., (4) we can define the reaction rate density as: R R Reactions R'= = = n.v.n.σ [ ] V A.x cm.sec (5) From this equation we can infer the units of the proportionality constant as:
' R [Reactions/(cm.sec)] σ = n.v. N [neutrons/cm ].[cm/sec].(nuclei/cm ) or:[cm ] (6) Thus the proportionality constant σ has units of area and physically represents the area that a nucleus in the target presents to the interacting neutrons in the impinging beam. The number of particles in the beam crossing a unit area per unit time is designated as the beam intensity or the beam current, and is given by: neutrons cm neutrons I = n.v. or: cm sec (cm.sec) (7) Thus the microscopic cross section can be written as: ' R [Reactions/(cm.sec)] σ = IN. [neutrons/(cm.sec)].(nuclei/cm ), or :[cm ] (8) which identifies it as the reaction rate density per unit beam intensity per nucleus in the target per unit volume, and the overall unit is a unit of area: cm. However, the characteristic size of a nucleus is in the range of 1-1 [cm], and accordingly, neutron cross sections data are more conveniently expressed in terms of the square of this characteristic distance in the barn unit, where: -4 1[barn]=1 [ cm ]. In this case reaction rate densities can be estimated as: ' R = INσ.. (9) This allows an interpretation of the microscopic cross section as the reaction rate density per unit beam intensity per nucleus per cubic centimeter of the target. The cross section is not in general equal to the actual area of the nucleus. For instance the radiative capture cross section for Au 197 at the peak of 4.9 ev resonance is x1 - cm, whereas the geometrical area of its nucleus is just 1.98x1-4 cm. The reaction cross section is much greater than the physical cross section of the nucleus, except at very high neutron energies where the cross section becomes of the same order of magnitude as the nucleus. This can be calculated from the knowledge about the empirically determined expression for the radius of the nucleus as: r = r A 1 =, A is the mass number 1 r 1.5 1 [cm] (1)
The cross sectional area of a nucleus from Eqn. 1 becomes: s = πr = πr A [ cm ] (11) For Au 197, the area of the nucleus becomes: 1 s = π (1.5 1 ) (197) = 4 1.98 1 [ cm ] Since the cross section has such a small magnitude, the unit of the barn was chosen jokingly to express such a small area with the name of a large structure. Incidentally, the characteristic size of the atom is given by the first Bohr s radius as: - α 7.975 1 α = = = = 7 4πR 4π 1.977 1 α is the fine structure constant = 7.975 1 R 7 1 is the Rydberg constant = 1.977 1 [ m ] 11 9 5.9177 1 [ m] 5.9177 1 [ cm] - Thus the nuclear to atomic radius ratio is: 1 r 1.5 1 = =.555 1 9 α 5.9177 1 MICROSCOPIC REACTION CROSS SECTIONS Each probable reaction that a neutron can undergo with a nucleus is associated with a specific cross section. The most important of them are: 5 σ se = elastic scattering cross section σ si = inelastic scattering cross section σ γ = radiative capture cross section σ f = fission cross section σ p = (n, p) reaction cross section σ T = (n, T) reaction cross section σ = (n, α) reaction cross section α The sum of the cross sections that can lead to the disappearance of the neutron is designated as the absorption cross section: σ = σ + σ + σ + σ + σ +... (1) a γ f p T α
The sum of the all reactions is designated as the total cross section: σ = σ + σ + σ + σ + σ + σ + σ +... t se si γ f p T = σ + σ + σ se si a α (1) The capture cross section incorporates all the cross sections that do not lead to fission: Thus, for a fissile nucleus: σ = σ + σ + σ + σ +... (14) c γ p T σa = σc + σ f (15) The microscopic cross sections are tabulated as a function of energy and are compiled in data bases such as the Evaluated Nuclear Data File at Brookhaven National Laboratory (BNL) as ENDF-B VI, the latest number being the version number. Cross sections plots for a light element, O 16 are shown in Fig. 4, and for a heavy element, Pb 8 are shown in Fig. 5. α Figure 4. Neutron reactions in a light element, O 16.
Figure 5. Neutron reactions in a heavy element, Pb 8. The cross sections for U 5 and U 8 are shown in Figs. 6 and 7. These graphs display a large number of resonances. Figure 6. Actual neutron reactions in U 5.
Figure 7. Actual neutron reactions in U 8. To calculate reaction rates of interest such as tritium production from Li 6 (Fig. 8) and from Li 7 (Fig. 9), and the threshold neutron multiplying reactions from U 8 (Fig. 1), libraries of reaction cross sections or response functions are generated. Figure 8. Reaction cross sections for tritium production in Li 6.
Figure 9. Reaction cross sections for tritium production in Li 7. Figure 1. Threshold reactions in U 8.
ESTIMATION OF NUMBER DENSITIES The cross section is normally determined experimentally, and the number density N of the target material, can be calculated from a modified form of Avogadro s law as: ρ atoms or nuclei N = A v M cm (16) which is a generalization of the conventional form of Avogadro s law: N ' g = Av [ atoms or molecules] (16) M where: ρ is the density of the target in [gm/cm ], g is the mass of the target in [gm], M is the atomic or molecular weight in atomic mass units in [gm/mole] or [amu], A v =.6x1 4 [nuclei, atoms, or molecules / mole], is Avogadro s number. Four practical situations are encountered in the estimation of number density of materials. CASE I: PURE SUBSTANCES AND SINGLE SPECIES. The modified form of Avogadro s can be applied in a straight forward way in the case of single species materials. EXAMPLE 1, Single species For Sodium (Na), ρ(na)=.97[gm/cm ], M(Na)=.99[gm/mole], ρ(na).97 M(Na).99 4 4 N(Na)= A v=.6 1 =.54 When we do not have a single species or a pure substance, other cases present themselves in practice. CASE II: MOLECULAR AND ISOTOPIC COMPOSITIONS. ATOMIC PERCENTAGES, (a/). Given: 1. ρ mixture of elements. Atomic weights of constituent elements. Mixture proportions as: i) Molecular compositions, ii) Atomic percentages (a/o).
In this case: ρ N =.A.f mixture element v a Mmixture Number of atoms of element f= a Atom or molecule of mixture (17) EXAMPLE, Molecular composition In the case of molecules such as water, ρ mixture =ρ(ho)=1. [gm/cm ], M =M(H O)=( 1.797+1 15.994)=18.15 mixture.6x1 18.15 4 4 N(HO)=1 [1]=.4 1 [molecules/cm ] 4 N(O) =N(HO) [1]=.4 N(H) =N(H 4 O) []=.6686 EXAMPLE, Isotopic composition Isotopic abundances are normally reported in atomic percentages (a/o). Consider the case of Boron. ρ =ρ(b)=.45 [gm/cm ] mixture 1 11 M(B )=1., M(B )=11. 1 11 a/o(b ) 19.8 %, a/o(b ) 8. % M =.198 1 +.8 11 = 1.8 mixture = = 4 1.6 1 4 N(B )=.45 [.198]=.71 1.8 4 11.6 1 4 N(B )=.45 [.8]=.11 1.8 EXAMPLE, Mixed Molecules and isotopic abundances We consider the case of boron carbide, B 4 C.
ρ mixture =ρ(b4c)=.5 [gm/cm ] M(B)=1.81, M(C)=1.1 1 11 a/o(b ) = 19.8 %, a/o(b ) = 8. % 1 1 a/o(c ) = 98.89 %, a/o(c ) = 1.11 % M =4 1.81 +1.1 = 55.5 mixture 4.6 1 4 N(B)=.5 [4]=.199 55.5 4.6 1 4 N(C)=.5 [1]=.75 55.5 1 4 N(B )=.198 N(B)=.18 11 4 N(B )=.8 N(B)=.881 1 4 N(C )=.9889 N(C)=.7 1 4 N(C )=.111 N(C)=. CASE III: COMPOSITION BY WEIGHT, ALLOYS AND ENRICHMENT, (w/o). Given: 1. ρ mixture. Atomic weights of constituent elements. Mixture proportions as weight percentages (w/o) In this case: mixture w N element =.A v, Melement w ρ.f f =weight percentage of element or: ρ (element) effective N element =.Av Melement (18) EXAMPLE 4, Enrichment given in w/o We consider the case uranium enriched to w/o in the U 5 isotope.
ρ mixture =ρ(u)=19.1 [gm/cm ] 4 5.6 1 1 N(U )=(. 19.1) =9.79 5.49 4.6x1 =.86 8.58 5 ρ effective (U )=. ρ(u)=. 19.1 [gm/cm ] 8 ρ effective (U )=.8 ρ(u)=.8 19.1 [gm/cm ] 8 N(U )=(.8 19.1) EXAMPLE 5, Alloy compositions We consider Zircaloy-4 used as cladding material in fission reactors, with a density: ρ(zircaloy-4)=6.74 [gm/cm ], and composition: 98.4 w/o Zr, M(Zr)=91..1 w/o Cr, M(Cr)=51.996.1 w/o Fe, M(Fe)=55.847 1.45 w/o Sn, M(Sn)=118.69 4.6x1 N(Zr)=(6.745.984) =4.7 91. 4.6x1 19 N(Cr)=(6.745.1) =7.81 51.996 4.6 1 N(Fe)=(6.745.1) =1.5 55.847 4.6 1 N(Sn)=(6.745.145) =4.96 118.69 EXAMPLE 6, Aqueous Solutions (Mixed Case I and Case II) We consider the case of Sulfuric Acid 1 w/ in solution in water, with a density: ρ(h SO +H O)=1.8 [gm/cm ] 4
ρ effective (HO) =.9 1.8 [gm/cm ] ρ effective (HSO 4)=.1 1.8 [gm/cm ] 4.6 1 N(HO) =(.9 1.8) =.5 18 4.6 1 N(HSO 4)=(.1 1.8) =6.68 98 N(H)=.N(H O)+.N(H SO )=.5 1 + 6.68 1 =6.68 4 N(O)=1.N(HO)+4.N(HSO 4)=1.5 1 +4 6.68 1 =.518 N(S)=1.N(HSO 4)=6.68 N(H )=.99985 N(H)=6.67 1 1 [atoms/cm ] 18 N(D )=.15 N(H)=9.956 16 N(O )=.99756 N(O)=.59 17 19 N(O )=.9 N(O)=1.7 N(S )=.95 N(S)=6.6 N(S )=.75 N(S)=4.978 1 18 [atoms/cm ] 4 19 N(S )=.4 N(S)=.788 6 16 N(S )=.15 N(S)=9.956 CASE IV: HETEROGENEOUS SYSTEMS. VOLUMETRIC MIXTURES, (v/o). Given: 1. Volume fractions (v/o) of either cases I or II.. Densities of components of volume mix. In this case: * N (element)=n(element in primary mix) (Volume fraction of primary mix in secondary mix) (19) EXAMPLE 7, Volumetric mixture Here we consider a secondary mix of 1 v/ Zircaloy-4 and 9 v/ water. The Zircaloy-4 and water can be considered as primary mixes and we use their results from examples 1 and 5:
* N (Zr)=.1 N(Zr) * N (Cr)=.1 N(Cr) * N (Sn)=.1 N(Sn) * N (Fe)=.1 N(Fe) * N (H)=.9 N(H) * N (O)=.9 N(O) EXAMPLE 8: Volumetric Homogenization. We consider the homogenization of a uranium dioxide fuel rod of.7 cm diameter covered with stainless steel, considered as Fe, cladding of.5 cm thickness. The densities of these materials are: ρ(uo )=1.[gm/cm ],ρ(fe)=7.8[gm/cm ] We can calculate the following number densities; 7.8 55.85 1. 4 N(UO )=.6 1 =. 1 [molecules/cm (UO )] (8+) 4 N(Fe)=.6 1 =8.41 N(U)=N(UO )=. 1 [atoms/cm (UO )] N(O)=N(UO )=4.46 1 [atoms/cm (UO )] π 4 π 4 V(UO )= (.7) =.849[cm ] V(Fe)= [(.8) -(.7) ]=.1178[cm ] The homogenized number densities become: V(UO ).849 V(UO )+V(Fe).849+.1178 * N (U)=N(U) =. 1 =1.71 * * N (O)=N (U) =.4 V(Fe).1178 V(UO )+V(Fe).849+.1178 * N (Fe)=N(Fe) =8.41 1 =1.97
EXAMPLE 9 We consider some uranium dioxide fuel of density ρ ( UO ) = 1.5[ gm / cm ], that has its uranium enriched to ε = w/o in the U 5 isotope. To calculate the number density of U 5 we use the formula: NU 5 ( ) ρ U = MU ( ) 5 effective ( ) 5 A v The effective density for U 5 is: 5 5 M ( U ) gm( UO ) gm( U ) gm( U ) ρeffective ( U ) = ρ( UO ).. ε M ( UO ) cm gm ( UO ) gm ( U ) We need to determine M(U) and M(UO ). To determine the atomic weight of U, we use: NU NU NU 5 8 ( ) = ( ) + ( ) ρ U ρ U ρ U MU ( ) MU ( ) MU ( ) 5 8 ( ) ( ) ( ) Av = A 5 v + A 8 v Since : ρu 5 ( ) =. ( ) ρ U 8 ( ) = (1 ). ( ) Then : ερu ε ρ U 1 ε (1 ε) = + MU MU MU 5 8 ( ) ( ) ( ) From the last equation we can calculate the atomic mass of the enriched uranium as: MU ( ) = ε MU MU + ε 5 8 ( ). ( ) 8 5. MU ( ) (1 ). MU ( ) 5.4 8.58 = (. x 8.58) + (.7 x 5.4) = 7.14 We can also evaluate the molecular mass of uranium dioxide as: MUO ( ) = MU ( ) + MO ( ) = 7.14 + ( 15.994)=69.14
From which the number density of U5 can finally be calculated as: N U 5 ( ) ρ = MU ( ) 5 effective ( U ) 5 Av MU ( ) ρ( ).. = MU ( ) UO ε M ( UO ) A 5 v 7.14 1.5. = 69.14.6 1 5.4 = atoms cm 1 7.11 1 [ / ] 4 MACROSCOPIC CROSS SECTION AND MEAN FREE PATH The sum of the microscopic cross sections of the individual nuclei in the target per unit volume is designated as the macroscopic cross section and is given by: cm Σ= N or cm cm 1. σ :[ ] () If a mixture of nuclear species exists in a unit volume, the overall macroscopic cross section becomes: 1 1. σ1. σ... n. σn[ ] Σ= N + N + + N cm () The macroscopic cross section can be conveniently estimated from the relation: ρ σ M 1 Σ=.6. [ cm ] () where the microscopic cross section σ is expressed in units of barns. EXAMPLE 1 is: The macroscopic scattering cross section for U 8 with a density of 19 [gm/cm ] Σ = ρ σ M 1 s.6. s[ cm ] 19 Σ s = = 8 1.6 8.9.46[ ] cm
EXAMPLE 11 Consider a uniform mixture of U 5 and graphite with moderator to fuel ratio of 1,. The macroscopic absorption cross section for carbon with a density of 1.6 gm/cm is: 1.6 Σ ac = = 1 1.6.4 1.7[ ] The number density of U 5 is a factor 1-4 that of the graphite, hence: 1.6 Σ au = = 1 cm 4 1 1.6 681.545[ ] The total macroscopic cross section for the moderator and fuel mixture becomes: cm Σ =Σ +Σ = + = 1 a ac au.7.545.57[ cm ] A simple expression for the reaction rate density is just the product of the beam intensity and the macroscopic cross section as: R reactions cm.sec ' = I. Σ (1) The inverse of the macroscopic cross section has units of length and has the physical meaning of representing the average distance traveled by a neutron between two interactions or mean free path : NEUTRON FLUX 1 mean free path: λ = [ cm] () Σ In the case of a radiation shield or a nuclear reactor neutrons acquire a random directional motion, which suggests the replacement of the velocity vector in the definition of beam intensity or current by the magnitude of the velocity or speed, which defines the neutron flux as: neutron φ = n v = nv cm.sec () which is now a scalar rather than a vectorial quantity. The neutron flux can be interpreted as the total distance traveled by the neutrons per unit volume per unit time [cm/(cm.sec)].
The neutron flux can also be considered as the number of neutrons per second entering an imaginary sphere of surface area of 4 cm and with a diametrical plane of area equal to 1 cm. The reaction rate density can still be expressed in term of the neutron flux as: R reactions φ cm.sec ' =. Σ (4) The reaction rate is now: reactions R= φ. Σ. V sec (5) where V is the volume under consideration in cm. In general, the reaction rate is an integral of the form: reactions R = ϕ( E, r ) Σ( E, r ) dedv sec (6) EV where the energy and spatial dependencies of both the neutron flux and the cross sections are accounted for. This integral is approximated as a summation in terms of the flux per unit energy interval as the form: reactions R= ϕi, jσi, j Ei Vj E V sec (7) i j For a constant cross section in each volume: reactions R = ( φiσi). V i sec where : φ = ϕ E i i, j i E i (8) The last equation is sometimes expressed in the inner product format: reactions R=< φ( E), Σ ( E) >. V sec (9) THERMAL NEUTRONS MAXWELLIAN VELOCITY DISTRIBUTION
When fast neutrons from the fission process collide with the reactor or shield material they loose they energy to the thermal equilibrium energy E = kt where they have as much probability of gaining energy as of loosing it through further collisions. The thermalized neutrons diffuse through the surrounding material until they are absorbed or leak from its surface. They have a population density or number of neutrons per cm that follows the Maxwellian distribution: n(v) = n 4π v ( π kt/m) / e mv kt () This distribution is normalized over all neutron velocities to n as: n mv 4π v kt e dv = n / ( π kt/m) where: k is the Boltmann constant=1.85x1 [ / ] 7 m is the rest mass of the neutron=1.675x1 [ ] n is the thermal neutrons density [n/cm ] o T = 7 + C,is the Kelvin temperature Joule K kg If we let: C = then we can write: 4π ( π kt/m) / n(v) = n Cv mv kt e (1) MOST PROBABLE VELOCITY The most probable velocity of thermal neutrons would occur at the maximum point of the Maxwellian distribution and can be obtained by setting the derivative of the neutron density distribution with respect to the velocity equal to zero. mv mv kt kt dn(v) mv = ncve ncv e = dv kt There follows that:
Thermal neutrons are also designated as, m/sec neutrons. MEAN NEUTRON VELOCITY The mean, average or mathematical expectation of the thermal neutrons velocity can be estimated from: v= = = = v n(v)dv n(v)dv v n n v n 4π v ( π kt/m) ( π kt/m) ( π kt/m) ( π kt/m) / 4π v / 4π v / 4π v n / e e e mv kt mv kt mv kt e mv kt where we applied the normalization condition for the Maxwellian distribution. To estimate the integral, let us make the change of variables: dv dv dv dv v u = du = vdv m c=- kt Thus we can write: v = = = 4π v ( π kt/m) π ( π kt/m) π ( π kt/m) / / / e mv kt ue ue cu cu dv du du
EXAMPLE 14 The absorption cross sections for the gold isotopes at.5 ev are: σ σ 197 Au a 198 Au a (9) = 98.8b (9) = 6,.b For a foil of gold irradiated in a reactor at a temperature of 5 o C: n + Au Au + γ 1 197 198 79 79 The Maxwellian and temperature corrected absorption cross sections of the two gold isotopes will be: σ σ 197 Au a π 9 ( T ) = 6, 7 + 5 198 197 Au Au a T σa 1/ 1/ π 9 = 98.8 = 98.8x.8867x.9548 = 8.94b 1/ π 9 ( ) = (9) T = 6,x.8867x.9548 = 1,946 b EXERCISE 1. Compare the area of the nucleus of U 8 to its neutron radiative capture at thermal (.5 ev) and fission (1.99 MeV) energies.. Calculate the macroscopic absorption cross section for natural uranium.. Estimate the mean free path for thermal neutrons scattering in beryllium. 4. A stainless steel composition is 69 w/o Fe, 17 w/o chromium, 1 w/o nickel and w/o molybdenum. Calculate its absorption cross section for thermal neutrons. 5. For, m/sec or thermal neutrons, calculate the following quantities: 1. Number densities,. Total macroscopic cross-sections,. Total mean free paths, In the following materials: 1. Uranium,. Beryllium,
. Carbon in the form of graphite. (Note that diamonds is a form of carbon with a high density). REFERENCES 1. M. Ragheb, Lecture Notes on Fission Reactors Design Theory, FSL-, University of Illinois, 198.. J. R. Lamarsh, Introduction to Nuclear Engineering, Addison-Wesley Publishing Company, 198.. Arthur R. Foster and Robert L. Wright, Jr., Basic Nuclear Engineering, Allyn and Bacon, 1978.