Submit manuscript...
eISSN: 2576-4543

Physics & Astronomy International Journal

Review Article Volume 2 Issue 6

A Gamma-ray in a uniform compton scatterer

Vladimir Uchaikin

Professor, Uljanovsk State University, Russia

Correspondence: Prof Vladimir Uchaikin, Uljanovsk State University, L.Tolstoi str. 42, 432000 Russia, Tel 79876875971

Received: February 01, 2018 | Published: December 6, 2018

Citation: Uchaikin V. A Gamma-ray in a uniform compton scatterer. Phys Astron Int J. 2018;2(6):573-576. DOI: 10.15406/paij.2018.02.00144

Download PDF

Abstract

The structure of a narrow ray of gamma radiation propagating through a medium filled with free electrons (called Compton scatterer) is investigated. The dual representation is used including primary (basic) and adjoint forms of kinetic equations. In the ray propagation problem primary form possesses cylindrical symmetry whereas the adjoint form (in case of isotropic detector in an unbounded homogeneous medium) has the spherical symmetry. The analysis performed on the base of adjoint function singularities shows, in particular, that in a vicinity of the narrow gamma-ray single-scattered radiation predominates over all other components. Analytical representation of the field in the vicinity of a cylindrical primary ray has been found. The result can be important in gamma astronomy processing.

Keywords: Compton Scatterer, neumann series, boltzmann equation

Introduction

The concept of adjoint function (or importance function) can be unlocked as follows. Let the problem under consideration requires finding some linear functional
J[f()]=(P,f)P(x)f(x)dxJ[f()]=(P,f)P(x)f(x)dx (1)
from a solution  of the equation

Lf(x)=Q(x),Lf(x)=Q(x),  (2)

Describing the transport of particles in a medium from a source with phase density Q(x)Q(x)  immersed in a medium, interaction of which with the particles is described by the operator L. In the problem under consideration,Q(x)Q(x)  may denote the flux of photons, or their concentration, or some other local characteristics of the photon field at a pointx=(r,p)x=(r,p) of the phase space. For the sake of convenience, we will represent the photon momentum pp through the pair of variables Ω=p/pΩ=p/p(the direction unit vector) and E=cpE=cp (photon’s energy)

The direct way of finding the value of the functional (1), expressing the reading of a photon detector, lies through solving of the kinetic equation (2) and computing integral J=f+(x)Q(x)dxJ=f+(x)Q(x)dx over the photon field with a weighting function P(x)P(x) expressing the contribution of a single photon placed at point X into resulting reading of the detector. But this is not a unique way to reach this result. Another way is based on the other equivalent representation of the detector reading,1−3 using so-called adjoint (or importance) function f+(x)f+(x) being a solution of the adjoint (in Lagrange’s sense) equation
L+f+(x)=P(x).L+f+(x)=P(x).  (3)

By definition, the primary LL and adjoint L+L+ operators are connected via relation

(f+,Lf)=(L+f+,f).(f+,Lf)=(L+f+,f).

Inserting now Eq.(2) into the LHS of the latter equation and Eq.(3) into the RHS of it, and using Eq.(1), we understand that numerically f+(x0)f+(x0)is the detector reading under condition that the source is localized at the point x0x0 and has a unite power: Q(x)=δ(xx0)Q(x)=δ(xx0). A more deep exposition of the concept of the importance function is given in books.13

Integral equations and neumann series

The mathematical description of penetration of gamma radiation (X-rays) through a scattering matter is given in the well-known monograph.4,5 Primary transport equation (2) is written as a linearized kinetic Boltzmann equation, for time-independent problem having the form:

ΩI(r,Ω,E)+μ(r,E)I(r,Ω,E)=Q(r,Ω,E)+dE'μs(r,E')I(r,ΩE')W(ΩE'Ω,E).ΩI(r,Ω,E)+μ(r,E)I(r,Ω,E)=Q(r,Ω,E)+dEμs(r,E)I(r,ΩE)W(ΩEΩ,E). (4)

Here I=I(r,Ω,E)I=I(r,Ω,E) denotes the differential intensity (so, IdSdΩdEIdSdΩdE is a mean number of particles, crossing over a unit elementary area dSΩdSΩ and belonging to intervals ×dE×dE), μμ and μsμs are linear attenuation (total and scattering) coefficients Q(r,Ω,E)Q(r,Ω,E) is the source phase density, WW the scattering indicatrix, normalized to 1:

4πE'0dEW(ΩE'Ω,E)=1.4πE0dEW(ΩEΩ,E)=1.

The correspondent adjoint function is of the form:

ΩI+(r,Ω,E)+μ(r,E)I+(r,Ω,E)=P(r,Ω,E)+dΩdEμs(r,E)W(Ω,EΩ',E)I+(r,Ω',E).ΩI+(r,Ω,E)+μ(r,E)I+(r,Ω,E)=P(r,Ω,E)+dΩdEμs(r,E)W(Ω,EΩ',E)I+(r,Ω',E).(5)

Both the equation can be transformed to the equivalent integral forms:67

I(r,Ω,E)=0eτ(r,rΩξ,E)˜Q(rξΩ,Ω,E)dξI(r,Ω,E)=0eτ(r,rΩξ,E)˜Q(rξΩ,Ω,E)dξ  (6)

with

˜Q(r,Ω,E)=Q(r,Ω,E)+4πdΩ'EEdEI(r,Ω',E)μs(r,E)W(Ω',EΩ,E),  (7)

and

I+(r,Ω,E)=0eτ(r,r+Ωξ;E)˜P(r+Ωξ,Ω;E)dξ,  (8)

with

˜P(r,Ω,E)=P(r,Ω,E)+μs(r,E)4πdΩ'E0dEI+(r,Ω,E)W(Ω,EΩ',E).  (9)

In these equations, τ(r,r±ΩξE) denotes the optical distance between the points for γ-quanta with energy E indicated in the argument:

τ(r,r±Ωξ;E)=ξ0μ(r±Ωξ)dξ.  (10)

Inserting (7) into (6) and (9) into (8), we will see initial terms of the Neumann series representing solution of the integral equations. The terms with n=0 relate to non-scattered radiation, whereas the other terms describe contribution of scattered quanta (the integer n=1,2,3,...  indicate scattering multiplicity. In particular,

I0(r,Ω,E)=0eτ(r,rΩξ;E)Q(rΩξ,Ω,E)dξ,  (11)

I+0(r,Ω,E)=0eτ(r,r+Ωξ;E)P(r+Ωξ,Ω,E)dξ,  (12)

I(r,Ω,E)=n=0In(r,Ω,E),I+(r,Ω,E)=n=0I+n(r,Ω,E),

and next terms are computed by induction:

In(r,Ω,E)=0dξeτ(r,rΩξ;E)μs(rΩξ,E)×[4πdΩ'EdEW(Ω',EΩ,E)In1(rΩξ,Ω,E)],n=1,2,  (13)

I+n(r,Ω,E)=0dξeτ(r,r+Ωξ;E)μs(r+Ωξ,E)×[4πdΩ'E0dEW(Ω,EΩ',E)I+n1(r+Ωξ,Ω',E)],n=1,2,  (14)

A point mono-directional source

Now we consider the field created by a point source of gamma-quanta placed at the origin O of the Cartesian coordinate and directed along z-axis, and let the field be measured by a point isotropic detector showing the integral

J=4πdΩ0dEP(E)I(r1,Ω,E)

with  standing for the detector response to a single photon with energy E. Thus, the free term in the adjoint equation (5) takes the form

P(r,Ω,E)=P(E)δ(rr1),  (15)

Where as the primary equation (4) has it in the form

Q(r,Ω,E)=δ(r)δ(ΩΩ0)Q0(E)  (16)

Q0(E)  denotes the source energy spectrum). The medium is supposed to be homogeneous and isotropic, so a suitable choice of coordinate system is dictated by specifics of the source-detector arrangement. An important role in this is given to symmetries, allowing for the reduction of the number of independent variables in the problem and thereby facilitating its solution (Figure 1). As we saw above, the same physical problem is expressed in two different forms differing in the degree of symmetry: Eq.(15) exhibits spherical symmetry whereas Eq.(16) reveals cylindrical one. So, it is naturally to choose the spherical symmetry version, possessing higher degree of symmetry and requiring only one spatial variable|rr1|. In what follows, we put r1=0, i.e. combine the origin of coordinates with the detector position, use the vector r for the initial position of the quantum, and introduce the angle θ via relation:

cosθ=Ωrr.

On substitution (15) into Eq.(12), we obtain non-scattered term in adjoint function decomposition:

I+0(r,Ω,E)=12πr2eτ(r,r';E)P(E)δ(1cosθ).  (18)

For a homogeneous medium τ(r,r+Ωξ;E)=μ(E)ξ , so

I+0(r,Ω,E)=12πr2eμ(E)ξP(E)δ(1cosθ).  (19)

Figure 1Geometry of narrow ray. Q is the source place, P is a counter place.

Computing single scattered intensity, we take into account that the scattering indicatrix includes the delta-function expressing the Compton interrelation between initial quantum energy E, scattering angle θ'  and scattered quantum energy

E1=E1+(1cosθ)E/mc2.

Denoting r+Ωξ by r' and the regular multiplier in the indicatrix by V(cosθ) , we represent the latter as

W(ΩΩ')=V(cosθ)δ(EE1(cosθ,E)),

where

cosθ=Ωr'r  (20)

and

V(cosθ)=1σdσdΩ

is the normalized differential cross-section of the Compton scattering. In an explicit form,

dσdΩ=r202(E1E)2(EE1+E1Esin2θ)

(Klein-Nishina-Tamm formula) and

σε=34σ0{1+εε3[2ε(1+ε)1+2εln(1+2ε)]}+34σ0{12εln(1+2ε)1+3ε(1+2ε)2},ε=Emec2.  (21)

Inserting n=1 into Eq.(14), using Eq.(19) and property of delta-function, we arrive at the following expression for single scattering contribution into the sought quantity

J1=I+(r,Ω,E)=0eμ(E)ξμ(E1)rμs(E)V(cosθ)P(E1)dξr2.

With fixed variables r and Ω, the variable ξ is uniquely linked with the scattering angle θ'(Figure 1). Applying the sinus-theorem yields interrelations

r=rsinθsinθξ=rsin(θθ)sinθ,dξ=rsinθsin2θdθ.

Using these formulae for passage to other variables leads us to expression

J=μs(E)rsinθθsθexp{rsinθ[μ(E)sin(θθ)+μ(E1)sinθ]}V(cosθ)P(E1)dθ.  (22)

In case of an infinite homogeneous medium (coefficients μ and μs are independent of spatial variables), the upper limit of integration θs=, but this formula can be used and on the plain boundary with an absorbing medium, when the source is in a homogeneous semi-space and the point of measurement is on a plane surface perpendicular to direction Ω . In this case the lower limit of integrating is zero upper limit θs should be taken π/2

The formulae for single-scattered quanta obtained above, includes multiplier [rsinθ]1, which grows unboundedly with decreasing 0. Consequently, at small θ or r the rest terms (with n>1) may be ignored.

Singularity of the adjoint function at great depths

Let us pass to the cylindrical coordinate system each point of which is characterized by longitudinal (z=rcosθ)  coordinate (depth) and transverse one (ρ=rsinθ) (radius). Then

I+1(r,Ω,E)=μs(E)ρθsθexp{μ(E)ρctgθμ(E)ρcosecθ}V(cosθ)P(E1)dθ.  (23)

Expand the exponential function in series and denote ith term of this series by φi:i=0,1,2,... :

φ0=μs(E)ρeμ(E)zθsθV(cosθ)P(E1)dθ, (24)

φ1=μs(E)eμ(E)zθsθ[μ(E)ctgθμ(E1)cosecθ]V(cosθ)P(E1)dθ,  (25)

φ2=12μs(E)ρeμ(E)zθsθ[μ(E)ctgθμ(E1)cosecθ]2V(cosθ)P(E1)dθ,  (26)

and so on.

Let the counter be far enough from the source and from the boundary (i.e. z>ρ and h<ρ). In this case, θρz0 and θsπρhπ , therefore,

φ0=μs(E)ρeμ(E)zπ0V(cosθ)P(E1)dθ.  (27)

The integrand in Eq.(25) has a singularity at θ=0 and θ=π. Let us break the integral into three pieces θsθ=ερ/z+πεε+πρ/hπε , where ε is a small positive number choosing in such a way that all terms in this sum are positive. Neglecting the change of μ,K and P in the ε-interval and leaving only singular terms, we obtain

φ1=μs(E)eμ(E)z[μ(E)+μ(Eπ)]V(cosπ)P(Ez)lnρ,  (28)

where

EπE1(cosπ,E).

Acting in a similar way, one can show that  and next terms in this sum at ρ0 do not have singularity.
If the counter is placed on the very boundary, then singularity characterizes only the first term (24) in the expansion of (23)

φ0=μs(E)ρeμ(E)zπ/20V(cosθ)P(E1)dθ.  (29)

Let us introduce the notations

Θ(ϕ,ψ;E)=2πP(E)ψϕV(cosθ)P(E1)dθ,

and Θ=Θ(0,ψ;E) , then for both above mentioned positions of the counter the expression

I+(r,Ω,E)I+0(r,Ω,E)+I+1(r,Ω,E)12πρ[δ(ρ)+μs(E)Θ]P(E)eμ(E)z,ρ0,  (30)

stay be valid, if put θs=π/2  in the boundary case, and π otherwise (more weak logarithmic singularity are omitted). As a result, we get:

I+(r,Ω,E;)=12πρ[δ(ρ)+μsΘ]P(E)eμz.

Let I0(x,y)  be the initial profile of the ray at z=0 and I(x,y,z) the ray profile at depth . According to the above result, the latter is composed from two terms: unscattered part repeating the shape of the incident ray, Insc(x,y,z)=I0(x,y)eμz, and scattered component

Isc(x,y,z)=μsΘeμzΣI0(x,y)dxdy(xx)2+(yy)2, (31)

whereΣ  stands for a cross section of the incident ray. The formula establishes the connection between the clear image of the ray given by the unscattered component, and the diffusion "halo," responsible for which is the scattering process. This connection is reminiscent of the relationship known in electrostatics between the charge distribution σ(x,y)=μsΘI0(x,y)eμz  on the plane z=const and the potential ϕ(x,y)=Isc(x,y) created by it in this plane.

In case of a uniform distribution of the initial flux over cross-section having the circular form with radius a,

σ(r)={σ0,r<a;0,r>0,

and

ϕ(r)=a02π0σ0rdrdϑr2+r22rrcosϑ.

Integration with respect to r' yields

ϕ(r)=2σ0π0(r2+a22racosϑr++rcosϑlnarcosϑ+r2+a22racosϑr(1cosϑ))dϑ.

Further, introducing

χ=πϑ2

and integrating the term with logarithm by parts, we obtain

ϕ(r)=2σ0[(a+r)E(k)+(ar)K(k)],k=2rar+a,  (32)

where K  and E stand for total elliptic integrals of the first and second kind respectively. At r=0, E=K=π2 and Eq. (32) becomes ϕ(0)=2πaσ0.

In the general case of an axially symmetrical initial profile, the haloϕ(r) Is expressed via integral

ϕ(r)=02π0σ(r)rdrdϑr2r22rrcosϑ. (33)

Expanding inverse distance with respect to Legendre polynomials,

1r2+r22rrcosϑ=1r>n=0(r<r>)nPn(cosϑ)

and inserting this into (33) yields:

ϕ(r)=n=0Cn[r0(rr)n+1σ(r)dr+r(rr)nσ(r)dr],

where coefficients

Cn=2π0Pn(cosϑ)=2πP2n(0)

are equal to 0 for odd n, and to 2π[12n(nn/2)]2 for even. Expression (33) can be integrated over the angle directly.

ϕ(r)=40K(k)σ(r)rdrr+r,k=2rrr+r.

In case of absence of axial symmetry the problem can be solved by numerical integration in expression (31).

Acknowledgments

This work is partially supported by the Russian Foundation of Basic Research (Projects 16-01-00556 and 18-51-53018).

Conflict of interest

Authors declare there is no conflicts of interest.

References

Creative Commons Attribution License

©2018 Uchaikin. This is an open access article distributed under the terms of the, which permits unrestricted use, distribution, and build upon your work non-commercially.