Calculating DPA Cheat Sheet

Displacement cross section integrals

For simple, hard sphere elastic interactions the maximum energy transfer, T^, is given as,

T^=γEi

where,

γ=4M1M2(M1+M2)2

See Chapter 1.1.1 of Was, 2nd Edition for derivations. We then need to compare T^ with the displacement energy, Ed and the electronic cutoff energy, Ec. Based on where T^ lands within these ranges, we can determine the displacement cross section, σd(Ei). For simple calculations assuming isotropic, elastic scattering using the Kinchin-Pease model we can use simple “if-then” logic to determine the displacement cross section:

If T^ < Ed, then

σd(Ei)=σs(Ei)γEiTˇT^υ(T) dT=σs(Ei)γEiTˇT^0 dT=0

 σD(Ei)=0

If  Ed<T^ <2Ed, then

σD(Ei)=σs(Ei)γEi(0Edυ(T) dT+EdT^υ(T) dT)=σs(Ei)γEi(0+EdT^1 dT)

σD(Ei)=σs(Ei)γEi(T^Ed)=σs(Ei)γEi(γEiEd)=σs(Ei)(1EdγEi)

 σD(Ei)=σs(Ei)(1EdγEi)

If  2Ed<T^ <2Ec, then

σD(Ei)=σs(Ei)γEi(0Edυ(T) dT+Ed2Edυ(T) dT+2EdT^υ(T) dT)

σD(Ei)=σs(Ei)γEi(0Ed0 dT+Ed2Ed1 dT+2EdT^T2EddT)

σD(Ei)=σs(Ei)γEi(0+Ed+T^2(2Ed)24Ed)=σs(Ei)T^(T^24Ed)

σD(Ei)=σs(Ei)(γEi4Ed)

If T^ >2Ec, then

σD(Ei)=σs(Ei)γEi(0Edυ(T) dT+Ed2Edυ(T) dT+2EdEcυ(T) dT+EcT^υ(T) dT)

σD(Ei)=σs(Ei)γEi(0Ed0 dT+Ed2Ed1 dT+2EdEcT2EddT+EcT^Ec2EddT)

σD(Ei)=σs(Ei)γEi(0+Ed+Ec2(2Ed)24Ed+Ec2Ed(T^Ec))

σD(Ei)=σs(Ei)γEi(0+Ed+12EdEc2(2Ed)22+Ec2Ed(γEiEc))

σD(Ei)=σs(Ei)(Ec2Ed)(1Ec2γEi)

The Kinchin-Pease model for displacements

The Kinchin-Pease (K-P) model is a simplified estimate of the number of displaced atoms (Nd) per PKA[^2] by means of the simple set of proposed relationships:

Nd={0for 0<E<Ed1for Ed<E<2EdE2Edfor 2Ed<E<EcEc2Edfor Ec<E 

Where E is the energy of the PKA and Ed is the threshold energy. At energies above Ec, the recoils lose energy only by electron excitation, while below Ec the energy transfer and loss is controlled by hardsphere elastic interactions (e.g. scattering). Note, the simplified model does not account for the effects due to the crystal lattice (e.g. crystal structure effects) for either single interactions or more complex damage cascades.