Open access peer-reviewed chapter

New Methods in Doppler Broadening Function Calculation

By Daniel Artur P. Palma, Alessandro da C. Gonçalves, Aquilino Senra Martinez and Amir Zacarias Mesquita

Submitted: May 16th 2012Reviewed: August 17th 2012Published: February 6th 2013

DOI: 10.5772/52464

Downloaded: 1532

1. Introduction

In all nuclear reactors some neutrons can be absorbed in the resonance region and, in the design of these reactors, an accurate treatment of the resonant absorptions is essential. Apart from that, the resonant absorption varies with fuel temperature, due to the Doppler broadening of the resonances (Stacey, 2001). The thermal agitation movement of the reactor core is adequately represented in microscopic cross-section of the neutron-core interaction through the Doppler Broadening function. This function is calculated numerically in modern systems for the calculation of macro-group constants, necessary to determine the power distribution in a nuclear reactor. This function has also been used for the approximate calculations of the resonance integrals in heterogeneous fuel cells (Campos and Martinez, 1989). It can also be applied to the calculation of self-shielding factors to correct the measurements of the microscopic cross-sections through the activation technique (Shcherbakov and Harada, 2002). In these types of application we can point out the need to develop precise analytical approximations for the Doppler broadening function to be used in the codes that calculates the values of this function. Tables generated from such codes are not convenient for some applications and experimental data processing.

This chapter will present a brief retrospective look at the calculation methodologies for the Doppler broadening function as well as the recent advances in the development of simple and precise analytical expressions based on the approximations of Beth-Plackzec according to the formalism of Briet-Wigner.


2. The Doppler broadening function

Let us consider a medium with a temperature where the target nuclei are in thermal movement. In a state of thermal equilibrium for a temperature T, the velocities are distributed according to Maxwell-Boltzmann distribution (Duderstadt and Hamilton, 1976),


where Nis the total number of nucleus, Mis the mass of the nucleus and kis Boltzmann’s constant.

Considering the neutrons as an ideal gas in thermal equilibrium, it is possible to write the average cross-section for neutron-nucleus interaction taking into consideration the movement of the neutrons and of the nucleus as:


where f(V)is the distribution function of Maxwell-Boltzmann as given by equation (1) and V=VΩ^is the velocity of the target nuclei. Denoting vr=vVthe relative velocity between the movement of the neutron and the movement of the target nucleus and considering the isotropic case, that is, with no privileged direction, it is possible to separate the integration contained in equation (2) in the double integral:


It is possible to see clearly in equation (3) that the cross-section depends of the relative velocity between the neutrons and the target nuclei. As the nuclei are in thermal movement, the relative velocity can increase or decrease. This difference between relative velocities rises to the Doppler deviation effect in cross-section behaviour. After integrating equation (3) in relation to the azimuthal angle (ϕ)the average cross-section for neutron-nucleus interaction can be written thus:


Denoting μ=cosθso that dμ=sinθdθ, equation (4) takes the form of:


From the definition of the relative velocity one has the relation,


and, as a result,


With the aid of a simple substitution, using relations (6) and (27), equation (5) is thus written as:


In equation (8), the limits of integration are always positive due to the presence of the module. As a result, one should separate the integral found in equation (8) into two separate integrals, as follows,


It is possible to modify the limits of integration for equation (9) taking into account that the mass of the target nucleus is much larger than the mass of the incident neutron. In terms of relative velocity, equation (9) can be written as:


In replacing the expression of the Boltzmann distribution function, equation (1), in equation (10) one has:


where it was defined β2M2kT. Introducing the variables for reduced velocities ϖr=βvrand ϖ=βv, equation (11) is written by:


Integrating equation (12) in relation to V one gets to the expression:


For resonances (that is, for the energy levels of the composed nucleus) it is possible to describe the energy dependence of the absorption cross-section by a simple formula, valid for T=0K, known as Breit-Wigner formula for resonant capture, expressed in function of the energy of the centre-of-mass by,


where E0is the energy where the resonance occurs and ECMis the energy of the centre-of-mass of the neutron–nucleus system. Apart from that, we find in equation (14) the term σ0, that is the value of the total cross-section σtotal(E)in resonance energy E0that can be written in terms of the reduced wavelength ƛ0by:


where the statistical spin factor gis given by the expression:


where Iis the nuclear spin and Jis the total spin (Bell and Glasstone, 1970).

In replacing the expression (14) in equation (13) one finds an exact expression for the average cross-section, valid for any temperature:


In a system with two bodies it is possible to write the kinetic energy in the centre-of-mass system, by


where MR=mMm+Mis the reduced mass of the system.

For the problem at hand, of a neutron that is incident in a thermal equilibrium system with a temperature T, it is a good approximation to assume that vvr. Thus, the ratio between the kinetic energy of the incident neutron and the kinetic energy of the centre-of-mass system is thus written


where A is the atomic mass of the target core. Resulting:


and denoting β2=12vth2one finally obtains the expression for the cross-section of radioactive capture near any isolated resonance with an energy peak E0, as written by:




where vris the module for relative neutron-nucleus velocity, vis the module for neutron velocity, and


The Doppler width for resonance ΓDis expressed by:


All the other nuclear parameters listed below are well established in the literature,

  • A= mass number;

  • T= absolute temperature;

  • E= energy of incident neutron;

  • ECM= energy of centre-of-mass;

  • E0= energy where the resonance occurs;

  • Γ= total width of the resonance as measured in the lab coordinates;

  • ΓD=(4E0kT/A)1/2= Doppler width of resonance;

  • v= neutron velocity module;

  • vr=|vV|= module of the relative velocity between neutron movement and nucleus movement;

  • vth=2kTM= module of the velocity for each target nucleus.

3. The Bethe and Placzek approximations

The expression proposed by Bethe and Placzek for the Doppler broadening function ψ(x,ξ)is obtained from some approximations, as follows:

  1. one neglects the second exponential in equation (22), given that it decreases exponentially and is negligible in relation to first integral in equation (22) given that (v+vr)2>>(vvr)2.

  2. it is a good approximation to extend the lower limit for integration down to in equation (22), given that the ratio between the energy of neutron incidence and the practical width is big.

  3. being ECMthe energy of the system in the centre-of-mass system and Ethe energy of the incident neutron, the following relation is always met:


where it was denoted that η=ECMEE. Equation (25) can be expanded in a Taylor series and, to the first order, is written by


In terms of the masses and velocities, equation (26) is written as follows:


where MRis the reduced mass of the system. For heavy nucleus MRmand equation (27) can be written as:


so that,


In replacing approximation equation (29) in the remaining exponential of equation (22) one finally obtains the Doppler broadening function that will be approached in this chapter,


The approximations made in this section apply in almost all the practical cases, and are not applicable only in situations of low resonance energies (E<1eV) and very high temperatures.

4. Properties of the Doppler broadening function ψ(x,ξ)

The function ψ(x,ξ)as proposed by the approximation of Bethe and Placzek has an even parity, is strictly positive and undergoes a broadening as that variable ξdiminishes, that is, varies inversely with the absolute temperature of the medium. For low temperatures, that is, when temperature in the medium tend to zero, the Doppler broadening function can be represented as shown below:


Equation (31) is known as an asymptotic approximation of the Doppler broadening function. For high temperatures, that is, when the temperature of the medium tends to infinite, the Doppler broadening function can be represented through the Gaussian Function, given that:


Figure 1.

The Doppler broadening function forξ=0.05,0.15and0.5.

The area over the curve of the Doppler Broadening function is written as below and, as it consists of separable and known integers it is possible to write:


From equation (33) once concludes that the area over the curve of the Doppler Broadening function is constant for the intervals of temperature and energy of interest in thermal reactors. This property is valid even for broadened resonances as shown in Figure 1, considering the different values of variable ξ.

5. Analytical approximations for the Doppler broadening function

This section describes the main approximation methods for the Doppler broadening function, according to the approximation of Bethe and Placzek, equation (30).

5.1. Asymptotic expansion

A practical choice to calculate the Doppler broadening function is its asymptotic expression resulting from the expansion of the term 11+y2in equation (30) in a Taylor series around y=x.


In replacing equation (34) in equation (30) and integrating term by term, one obtains the following the asymptotic expansion:


Despite equation (35) being valid only for |x.ξ|>6, it is quite useful to determine the behaviour of the Doppler Broadening function in specific conditions. For high values of x, it is possible to observe that function ψ(x,ξ)presents the following asymptotic form:


5.2. Method of Beynon and Grant

Beynon and Grant (Beynon and Grant, 1963) proposed a calculation method for the Doppler broadening function that consists of expanding the exponential part of the integrand of the Doppler broadening function ψ(x,ξ)in the Chebyshev polynomials and integrate, term by term, the resulting expression, which allows writing:

ψ(a,b)=1a{πcos(ab)[1- E2(a)]ea2+J(a,b)}e-14b2,E37

where a=12ξand b=ξxand still,



E2(a)=2π0aey2 dy.E39

For values where the condition |x.ξ|>6was met (Beynon and Grant, 1963) is recommended the use of the asymptotic expression of the function ψ(x,ξ), equation (36). It should be pointed that the results obtained by this method have become a reference in several works on the Doppler broadening function.

5.3. Method of Campos and Martinez

The core idea of the method proposed by (Campos and Martinez, 1987) is to transform the Doppler broadening function from its integral form into a differential partial equation subjected to the initial conditions. Differentiating equation (30) in relation to xone obtains:


Acknowledging in equation (40) the very Doppler broadening function and the term of interference as defined by the integral:


it is possible to write:


Deriving equation (42) again in relation to x, after expliciting function χ(x,ξ)in the same equation, one has:


The right side of equation (43) can be written in another way, given that y21+y2=111+y2,


In replacing the result obtained in (44) in equation (43) one obtains the differential equation where Campos and Martinez based themselves to obtain an analytical approximation for the broadening function ψ(x,ξ):


subjected to the initial conditions:


Admitting that function ψ(x,ξ)may be expanded in series,

ψ(x,ξ)=n=0cn(ξ) xnE47

and in replacing-se equation (47) in the differential equation as given by equation (45), one obtains after some algebraic manipulation the following polynomial equation:





and all the other terms are calculated from the following relation of recurrence:


The representation in series for the Doppler broadening function, as given by equation (47), is valid only for |x.ξ|<6. For the cases where |x.ξ|>6, (Campos and Martinez, 1987) used the asymptotic form as given by equation (35), as well as proposed by Beynon and Grant.

5.4. Four order method of Padé

The Padé approximation is one of the most frequently used approximations for the calculation of the Doppler broadening function and its applications and can efficiently represent functions, through a rational approximation, that is, a ratio between polynomials. For the four-order Padé approximation (Keshavamurthy& Harish, 1993) they proposed the following polynomial ratio:


whose coefficients are given in Tables 1 and 2.


Table 1.

Coefficients pand qof the four-order Padé Approximation


Table 2.

Coefficients h, aand bof the four-order Padé Approximation

From the coefficients of Tables 1 and 2, and of equation (49), one obtains in the end the following analytical approximation for function ψ(x,ξ), according to the four-order Padé approximation:


where η(x,ξ)and ω(x,ξ)are the following polynomials:




5.5. Frobenius method

In this method the homogeneous part of the differential equation that rules the Doppler broadening function, equation (45), is solved using the Frobenius Method (Palma et. al., 2005) that consists fundamentally of seeking a solution of the differential equation in the form of series around the point x=x0, with a free parameter, that is, as follows:

ψ(x,ξ)=xsn=0cn(ξ) xn=n=0cn(ξ) xn+s,E53

with c00and where s is the parameter that grants the method flexibility.

Deriving equation (53) and replacing it in the homogeneous equation associated to equation (45) one obtains, after grouping the similar terms:

n=0cn(n+s)(n+s1) xn+s2+n=0cnξ2[(n+s)+ξ2+24] xn+s+ξ24n=0cn xn+s2=0.E54

The initial equation of the problem, obtained when n=0, remembering that c00is


From equation (55), as c00, one obtains that s=0or s=1. Using first s=0and c00one obtains the following relations of recurrence:

cn=ξ2(4n+ξ26)4n(n+1)cn2, valid forn=2orn=3E56
cn=ξ2[cn2(4n+ξ26)+cn4ξ2]4n(n+1),valid for  n4.         E57

Considering the case where s=1, one obtains the other series linearly independent with the first term, not null, denoted by c˜0:

cn=ξ2(4n+ξ22)4n(n+1)cn2, valid for  n4E58
cn=ξ2[cn2(4n+ξ22)+cn4ξ2]4n(n+1), valid for  n4E59

With this the homogeneous solution assumes the following form:


where the coefficients are all known from the relations of recurrence, equations (56) to (59). In writing function ψh(x,ξ)as:


it is possible to determine all the coefficients Anequalling, term by term, equations (60) and (61) so to write:


Acknowledging the expansion of the cosine and sine functions, one obtains an analytical form to solve the homogeneous part of the differential equations that rule the Doppler broadening function:


In order to obtain the particular solutions of equation (45), and consequently its general solution, it is possible to apply the method of parameter variation from the linearly independent solutions:


Supposing a solution thus,


where functions u1(x)and u2(x)are determined after the imposition of the initial conditions expressed by equations (46a) and (46b) and of the imposition of the nullity of the expression:


That, along with the condition,


Which results from the very equation (45), form a linear system whose solution is given by the equations:


Integrating equations (69) and (70),


it is possible to write the solution particular of equation (45) as follows:


As the general solution of differential equation (45) is the sum of the solution of the homogeneous and particular equations, the initial conditions are imposed, as expressed by equations (46a) and (46b), to determine the constants:


Finally, according to the Frobenius Method, the Doppler Broadening function can be written thus:


where ϕ(x,ξ)=erf(iξxξ2).

5.6. Fourier transform method

In doing the transformation of variables u=ξ2(xy)in the full representation of the Doppler broadening function, equation (30), one obtains the expression

ψ(ξ,x)=1π+eu2du1+(x2uξ)2 ,E77

that can be mathematically interpreted as the convolution of the Lorentzian function with a gaussian function, as exemplified by the equation below:

ψ(ξ,x)=f*g+g(u) f(x-u) du ,E78

where f(x-u)=11+(x2uξ)2is the lorentzian function and g(u)=1πeu2the gaussian function. Function f(x-u)admits a full representation through the Fourier cosine transform (Polyanin and Manzhirov, 1998), as being

f(x-u)=0ewcos[(x2uξ)w] dw.E79

In replacing-se equation (79) in the integer of convolution, as given by equation (78), applying the properties of the integrals of convolution one gets to the following expression:

ψ(ξ,x)=f*g 0e-w+g(u) cos[(x-2uξ )w] du dw =0e-w I(w) dw ,E80


I(w)  1π+eu2 cos[(x-2uξ) w] du=1πcos(xw)+eu2cos[(2uξ) w] du=e-wξ22cos(xw).E81

In replacing equation (81) in the equation (80), one obtains a new full representation of the Doppler broadening function, interpreted as a Fourier cosine transform (Gonçalves et. al., 2008):

ψ(ξ,x) = 0e-w2ξ2wcos(wx) dw= 12[ 0ew2ξ22wadw +  0ew2ξ22wbdw ],E82

where a(1ix)2and (1+ix)2.

The integrals on the right side of equations (3.25) and (3.26) are known as complementary error functions, in which case one can conclude that:

0ew2ξ22 wadw =ξπ2e(xi-1)24ξ2erfc (ξiξx2)E83
0ew2ξ22 wbdw =ξπ2e(xi+1)24ξ2erfc (ξ+iξx2).E84

In replacing equations (83) and (84) in equation (82) it is possible to write the following expression for the Doppler broadening function:

ψ(x,ξ)= ξπ4e(xi1)2ξ42[1+erf (iξxξ2)]+ξπ4e(xi+1)2ξ42[1erf (iξx+ξ2)].E85

With some algebraic manipulation it is easy to prove that the Fourier transform method and the Frobenius method, equations (85) and (76) respectively, provide identical results.

5.7. Fourier series method

From the representation of the Doppler broadening function in a Fourier cosine transform, equation (82), it is possible to write

ψ(ξ,x) = 0e-w2ξ2wcos(wx) dw=0G(w)ewcos(wx) dw,E86

where function G(w)=e-w2ξ2is even and can be expanded into a Fourier series in cosines:



an=ξπ2Le(nπξ2L)2[erf (2L+nπξ2i2ξL)+erf (2Lnπξ2i2ξL)].E89

In replacing equation (87) in equation (86) and integrand, it is possible to write the following expression for the Doppler broadening function in the form of Fourier series:

ψ(x,ξ)=ξπ2L(1+x2)erf(Lξ)+ξπLn=1Fn(x,ξ,L) Re [Z(ξ,L)],E90


Fn(x,ξ,L)[(nπ)2+L2(1+x2)]e (n π ξ2 L)2L2(1+x2)2+(nπ)2(22x2+(nπ/L)2),E91
Z(n,ξ,L)=erf (nπξ2i+2L22ξL).E92

5.8. Representation of function ψ(x,ξ)using Salzer expansions

Although the formulations obtained for function ψ(x,ξ)from the Frobenius method, Fourier transform and Fourier series methods only contain functions that are well-known in literature, it can be inconvenient to work with error functions that contain an imaginary argument. One of the ways to overcome this situation is to calculate the real and imaginary parts of function ϕ(x,ξ)using the expansions proposed by Salzer (Palma and Martinez, 2009)



Reϕ(x,ξ)erf(ξ2)+exp(ξ24)×{1πξ[cos(ξ2x2)1]+2πn=1nmaxexp(n2/4)n2+ξ2 fn(x,ξ)}E94
Imϕ(x,ξ)exp(ξ24){1πξsin(ξ2x2)+2πn=1nmaxexp(n2/4)n2+ξ2 gn(x,ξ)},E95

where auxiliary functions fn(x,ξ)and gn(x,ξ)are written by:


5.9. The Mamedov method

Mamedov (Mamedov, 2009) put forward an analytical formulation to calculate function ψ(x,ξ), based on its representation in the form of a Fourier transform, equation (82). Using the expansions in series of the exponential and cosine functions,


and the well-known binomial expansion


Mamedov proposed the following expressions for the Doppler broadening function:

for ξ>1and x>1


for ξ1and x15


for ξ1and x=0


where Γ(x,ξ), γ(x,ξ)and Γ(x)are the well-know incomplete Gamma functions and Fm(n)are binomials coefficients defined by:

Fm(n)={n(n1)...(nm+1)m!, for integer n(1)mΓ(mn)m!Γ(n), for non-integer nE104

6. Numerical calculation of function ψ(x,ξ)

The numerical calculation of the Doppler broadening function consists of calculating a defined integral. There are many methods in the literature for this calculation, but in this chapter we will describe a numerical reference method based on the Gauss-Legendre quadrature. In basic terms, the Gauss-Legendre quadrature method consists of approximating a defined integer through the following expression:


where Nis the order of the quadrature, ηiis the point of the quadrature and withe weight corresponding to the point of quadrature. The points of the Gauss-Legendre quadrature are the roots of the polynomials of Legendre (Arfken, 1985) in the interval [1,1], as generated from the Rodrigues‘ formula,


for an isotope at a given temperature, that is, for a fixed value for variable ξ, the function ψ(x,ξ)decreases rapidly and a very high value is not necessary for what we will consider our numerical infinite. This fact can be evidenced at Figure 1.

For that an adequate numerical infinite (x=5000)was considered, as well as a high-order quadrature (N=15), whose points of Legendre and respective weights are found in Table 3. The results obtained with this method, whose handicap is the high computing cost, can be seen in Table 4.


Table 3.

Points of Legendre ηiand respective wiweights.


Table 4.

Reference values for Doppler Broadening Function ψ(x,ξ).

7. Conclusion

A brief retrospective look at the calculation methodologies for the Doppler broadening function considering the approximations of Beth-Plackzec according to the formalism ofBriet-Wigner was presented in this chapter.


This research project is supported by Instituto Alberto Luiz Coimbra de Pós-graduação e PesquisaemEngenharia/Universidade Federal do Rio de Janeiro (COPPE/UFRJ), ComissãoNacional de Energia Nuclear (CNEN), Centro de Desenvolvimento de Tecnologia Nuclear (CDTN). It was financially supported by through National Institute of Science and Technology of Innovative Nuclear Reactors,Brazilian Council for Scientific and Technological Development (CNPq,) and Research Support Foundation of the State of Minas Gerais(FAPEMIG).

© 2013 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms of the Creative Commons Attribution 3.0 License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

How to cite and reference

Link to this chapter Copy to clipboard

Cite this chapter Copy to clipboard

Daniel Artur P. Palma, Alessandro da C. Gonçalves, Aquilino Senra Martinez and Amir Zacarias Mesquita (February 6th 2013). New Methods in Doppler Broadening Function Calculation, Current Research in Nuclear Reactor Technology in Brazil and Worldwide, Amir Zacarias Mesquita, IntechOpen, DOI: 10.5772/52464. Available from:

chapter statistics

1532total chapter downloads

More statistics for editors and authors

Login to your personal dashboard for more detailed statistics on your publications.

Access personal reporting

Related Content

This Book

Current Research in Nuclear Reactor Technology in Brazil and Worldwide

Edited by Amir Mesquita

Next chapter

Isothermal Phase Transformation of U-Zr-Nb Alloys for Advanced Nuclear Fuels

By Rafael Witter Dias Pais, Ana Maria Matildes dos Santos, Fernando Soares Lameiras and Wilmar Barbosa Ferraz

Related Book

First chapter

Experimental Investigation of Thermal Hydraulics in the IPR-R1 TRIGA Nuclear Reactor

By Amir Zacarias Mesquita, Daniel Artur P. Palma, Antonella Lombardi Costa, Cláubia Pereira, Maria Auxiliadora F. Veloso and Patrícia Amélia L. Reis

We are IntechOpen, the world's leading publisher of Open Access books. Built by scientists, for scientists. Our readership spans scientists, professors, researchers, librarians, and students, as well as business professionals. We share our knowledge and peer-reveiwed research papers with libraries, scientific and engineering societies, and also work with corporate R&D departments and government entities.

More About Us